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THE APPLICATION OF SEDIMENT-TRANSPORT MECHANICS 
TO STABLE-CHANNEL DESIGN 


Emmett M. Laursen, A.M. ASCE 
(Proc. Paper 1034) 


SYNOPSIS 


The three requisites for a stable alluvial channel are explicitly stated and 
he role of sediment transport in each assessed. The similarity of the many 
sediment transport formulas is demonstrated and the general method of ap- 
lication to design illustrated. The use of these formulas as scaling rela- 
ionships between different channels is advocated. 


Not enough is known about the mechanics of sediment transportation for 
heory alone to be used with confidence in the design of an alluvial channel. 
Jowever, enough is known for transport theory to be an aid tothe designer 
n the exercise of engineering judgment. In part this aid consists of dispel- 
ing confusion by demonstrating that design methods and criteria which are 
seemingly different are instead basically equivalent. 

By definition, a stable sediment-bearing canal in erodible material must 
1ot only convey the required water discharge, but must also transport the 
;ediment load supplied, and must not scour or silt either the bed or the banks. 
[These three factors neatly delineate the designer’s problem, the solution to 
vyhich is almost single valued—“almost” only because there may be a narrow 
wut finite range of flow conditions above which the banks will not silt and be- 
ow which they will not erode. The designer’s problem is the selection of the 
‘correct” combination of slope, size, and shape of channel which will result 

1 a stable canal for a particular combination of water discharge and sedi- 
nent load (or range thereof) and the particular materials through which the 
‘anal must pass. The selection must be based on personal experience guided 
y such relationships as have been deduced from general experience with 
xisting channels and with laboratory flumes. The difficulty in selecting the 
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correct combination stems from an always limited personal experience and 
from a lack of confidence in the guiding relationships. 

In order to design a stable channel, three relationships are needed, one 
for each of the factors in the previous definition—a flow equation, a sedimen 
transport equation, and a bank-erosion criterion. Although the flow condi- 
tions in the stable canal must satisfy all three of these relationships simul- 
taneously, for clarity in exposition they will be considered separately in the 
following discussion. 


The Flow Equation 


Of the three relationships needed for the design of a stable alluvial chan- 
nel, the best understood is the flow equation. Basically this relationship is 
the Chezy equation 


: 
: 
' 


Q=CAVRS (1) 


In the use of this equation the designer must assign, arbitrarily, a numerical 
value to the discharge coefficient C . In order to aid his judgment, he will 
usually prefer to use some empirical formula such as that of Manning: 


Q: 449 ares o/2 (2) 


This form of the flow equation was widely accepted because it was generally 
believed that the coefficient n was primarily if not entirely a function of the 
boundary material. Although it is now realized that n can be affected by 
several factors, a wealth of information has been accumulated over the years 
on the n values for existing channels covering a wide range of size, shape, 
alignment, bank and bed material, and condition. Because of this store of in- 
formation, n values for channels similar to that being designed can usually 
be found. The importance of this fact is self evident, but it is noteworthy 
here because the same sort of aid to judgment is needed for the other two 
relationships. 

Sediment transportation will affect the flow, and it will hence be reflected 
in the flow equation insofar as it affects the discharge coefficient. Dunes are 
almost always associated with the transport of sediment, and the resistance 
coefficient, which is the inverse of discharge coefficient, is affected by the 
resulting increase in boundary roughness. That sediment transport may in- 
fluence the resistance coefficient in other ways has been claimed(1) but has 
not been proved unequivocally. No matter what the causal relation, however, 
the rate of sediment transport is one of the several criteria which should be 
considered in appraising the similarity of channels, ; 

Two attempts to relate the resistance coefficient to sediment transport — 
should be mentioned. Einstein and Barbarossa(2) divided the bed resistance 
into two parts, that due to the sand roughness and that due to bar, or dune, 
roughness. To evaluate the sand roughness they used Strickler’s formula 


n=VKe /29.3 (3) 
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ere K, is the diameter (measured in feet) of the particle size coarser 
han 65% of the sediment. After eliminating bank effects in narrow channels, 
e effect of the sand roughness was subtracted from the total resistance and 
his remainder was related to a modification of Einstein’s y function for the 
ransportation of bed load. 

Barton and Lin(3) have proposed a relationship between the Chezy coef- 
icient, the Reynolds number of the flow, and the ratio between the sediment 
ize and the hydraulic radius. Although sediment transport does not enter the 
elationship explicitly, the data were obtained from flows transporting sedi- 
ent. 

Both of these proposals are so new that there has not been sufficient time 
0 evaluate their usefulness or their weaknesses. Both of them, however, by- 
ass the problem of predicting dune size, which is probably the major effect 
sediment transport on resistance to flow. Anderson(4) recently obtained 
solution for the length of dune as a function of the Froude number and pre- 
ented a limited amount of data corroborating his mathematical solution. His 
olution had a limit of a plane bed where the dune length became infinite. In 

ddition he presented one example of a plane bed at a Froude number of 1.6. 

measurements at the Iowa Institute, however, a plane bed (i.e., without 

lunes) has been obtained at a Froude number as low as 0.85. The sand in this 
xperiment was 0.1 mm in diameter, finer than the sand of Anderson’s exam- 
ple. Of course, it is possible for dunes to disappear because the dune height 
educes to zero as well as because the dune length increases in infinity. The 
eight of the dune, which was not and cannot be a part of Anderson’s analysis, 
is fully as important as the length. Very likely it is more a function of the 
article size and the Reynolds number of the flow than of the Froude number. 
The empirical equations describing stable channels which have been de- 
| in India(5) include relationships which are in essence evaluations of 
the Chezy C , just as the Manning equation is. These flow relationships 
turally are subject to the same criticism as is often made of the other 

dian criteria—that since they are based on data from Indian canals they are 

ilored to fit Indian conditions. This criticism, although valid, does not 
mean that due consideration should not be given to modifying the relation- 
ships to fit the conditions in other areas of the world. 
| Whatever flow equations and whatever numerical value of aresistance co- 
efficient the designer chooses, he will arbitrarily, though to the best of his 
judgment, have placed a restriction on the physical character of the channel 
eventually to be selected. Before the final decision two more restrictions in- 
volving the same type of evaluation must be explicitly stated. 


: The Sediment-Transport Equation 
| 


- If the flow entering the canal is not clear, but carries with it a sediment 
oad, the capacity of the channel must exactly equal the rate at which sediment 
is supplied to the channel for regime conditions to prevail. This requirement, 
then, places a second restriction on the physical characteristics of the chan- 
1el to be selected. Some of the many sediment-transport equations for bed 
oad that have been proposed are listed in Table I. This table is not com- 
prehensive nor is it intended to rate or judge the many expressions that can 
be found in the literature. In fact, the problem of picking a sediment-trans- 
port equation is an unhappy one that must be faced by each individual designer. 

i 

* 
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All of the equations are slightly different in form; all are based on limited 

data, either entirely or largely from the laboratory; each fits some data bet 
ter than the others; none fits all the data (which is why there are so many); 
none has gained the wide acceptance shared by the relatively few flow equa- 
tions. If unable to decide with confidence on one “best” equation, the design 
may use more than one, or even all, in the same manner as the hydrologist 
wont to do with the many flood formulas at hand. 

In Table I, it should be noted that the equations suggested by Einstein(12) 
and Kalinske(13) in their original papers have not been used, but instead the 
modifications suggested in the sediment chapter of Engineering Hydraulics. 
In these modifications the original parameters are used, but the forms of th 
equations are changed to fit the data for the high rates of sediment transpor 
That these modified equations are not strikingly different from the other en 
pirical sediment-transport equations should not be surprising, because in 
their development the assumption is made that the velocity at the edge of th 
laminar sublayer vs (whether or not actually existent) can be substituted 
for the velocity at the level of the moving grain of sediment Vg. Since vg 
is, as usual, equated to 11.6 /7/p_ and bears little relation to Vg except 
that both are small and both are at levels near the bed, this is tantamount tc 
assuming that qs is a function of rT. 

In order to compare the different transport equations more easily, they 
have been put into a form indicating the dependence of qs on the mean velo 
ity V and the depth of flow y by neglecting the critical shear Tc and usin 
the Manning equation to eliminate the slope. The similarity between not onl 
these but all the commonly used sediment-transport equations can thus be 
made apparent. One might also suspect, along with Gilbert, (14) that the ex- 
ponent as well as the coefficients in the transport equation should be variabl 

If a sizable bed load is to be expected, the simplified forms can be used t 
estimate the transport characteristics of the canal. However, since it is dif 
ficult to measure the bed load or to predict the bed-material load entering 
through the headworks of proposed canal, a further modification of the form. 
may be helpful. For example, if the first equation in Table I is written for 
the sediment discharge ratio, and then the ratios of the various parameters 
are taken for the proposed and an existing “similar” channel, relative rathe: 
than absolute values can be estimated: 


ds/q (proposed) By (p)_n4(p) V3(p) y75/3 (p) 


(4) 
qs/q (existing) B,(e) n#(e) v>(e) y5/3 (e) 


It is interesting to note that Chezy proposed his flow equation in this ratio 
form, (15) and that the accumulation of n values together with descriptions « 
channels is essentially only a somewhat more sophisticated manner of han- 
dling the problem. 

A further rearrangement of the simplified transport equations will put 
them in the form of the Kennedy equation and the variations thereof as show! 
in Table IL The similarity between the modified transport equations and the 
Kennedy type of equations indicates that the physical basis for the latter equ 
tions is the requirement for sediment-transport Capacity. That only Inglis — 
has considered this factor explicitly is probably due to the fact that the equa 
tions were derived from data of canal systems in which the sediment- 4 
discharge ratio did not vary greatly. 


Re Soyo ee 
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On the basis of personal judgment, aided to some extent by the sediment- 
transport equations of one form or another, the requirement for sediment 
capacity must be satisfied, just as the requirement for flow capacity was 
satisfied, and a second restriction on the combination of size, shape, and 
slope of channel is imposed. 


The Bank-Erosion Criterion 


If the banks of the canal were vertical and non-erodible, the final selection 
of channel characteristics satisfying only the first two requirements could be 
made either arbitrarily or on the basis of economics. An alluvial channel 
will have sloping banks of erodible material. Flow conditions along the banks, 
therefore, must be such that sediment will not deposit on the banks and that 
bank material will not be scoured. The first of these requirements will 
govern for canals designed for a minimum slope; the second for canals de- 
signed for minimum area. Satisfactory solution of the non-silting problem 
requires an understanding of sediment-transport mechanics in much greater 


detail than is presently available, especially of the suspended-sediment load. 


Knowledge concerning the non-scouring problem is somewhat more advanced, 


but the bank-erosion criterion is still the least satisfactory of the three re- 


lationships needed for design. 

In the case of banks composed of granular material (e.g., sands and 
gravels) erosion will occur if the forces tending to move the particles are 
sufficiently large to overcome the resisting force of gravity. In the case of 
banks composed of cohesive material (e.g., clays) the force of cohesion will 
also resist movement and in fact may be much greater than the resisting 
force of gravity. 

In associating the flow conditions in the canal to the bank erodibility the 
canal shape enters the problem on both sides of the relationship. The slope 
of the banks will obviously affect the susceptibility of the banks to erosion 
and the cross-sectional shape will affect the distribution of the boundary 
shear. 

Lane and Carlson(20) have presented the relationship for the ratio of the 

critical tractive force on a bank slope Tcs to the critical tractive force Teh 
as a function of the angle of the side slope and the angle of repose of the ma- 


terial. That is, 
2 
eset tan<¢ 
K= rep = 008 4/1 - tonzg 


where ¢ is the angle of the side slope, and 6 is the angle of repose. Other 
‘studies under the direction of Lane(21) have attempted to obtain the distribu- 
tion of shear in a trapezoidal section. This knowledge of the effect of cross- 
‘sectional shape allows data or equations on the critical tractive force for bed 
material to be transferred to the banks. 

The two formulas for critical tractive force most widely quoted are those 


“of White(22) 


Tob = 12 D (6) 
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and of Shields(10) 
Tb = (Ss-!) D f(IL6 D/8) (7) 


Any of the relationships for Te of the equations in Table I can also be used— 
at the discretion of the designer. It should be noted, however, that these re- 
lationships may result in a low rate of movement rather than zero movement, 
A criterion based on zero movement is, perhaps, unduly restrictive, because 
a small amount of scour can usually be tolerated and because suspended sedi- 
ment will tend to add cohesive material to the banks. Application of the rela- 
tionship from the first equation in Table I will result in values of permissible 
canal velocities approximately the same as given by Fortier and Scobey; (23) 
application of the relationships of White and Shields will result in smaller 
permissible velocities. 

Very little is known concerning the resistance to movement of clays—ex- 
cept that some clays can stand up under very high boundary shear. One might 
speculate that the resistance to boundary shear! should be related to the in- 
ternal shearing characteristics of the clay mass in an undisturbed saturated 
condition. For naturally deposited material, which is not likely to be com- 
pletely homogeneous, the weakest part of the mass would govern. Erosion of 
seams or pockets of weak material would, of course, expose the adjacent 
stronger material to greater forces. The only published work concerning the 
resistance of cohesive materials is contained in Fortier and Scobey, and even 
there insufficient information is given about the character of the material to 
serve as more than a rough guide for the designer. 

This last relationship, together with the first and second, fixes the width, 
depth, and slope of the channel. No method is available to determine the best 
geometrical shape of the cross section—except examination of other channels. 
In fact, canals are usually dug to a shape easily handled by construction ma- 
chinery, and then allowed to modify that shape by a moderate amount of scour 
and fill. 


CONCLUSIONS 


Three relationships are needed to specify a stable alluvial channel. The 
three relationships may be those directly stating the three requirements for 
the stable canal as presented herein, or they may be any three independent 
forms derivable therefrom. The three independent regime-theory criteria 
implicitly contain these same factors. 

The difficulties of designing a stable canal lie not in using the general re- 
lationships, but in assigning numerical values to coefficients and even in 
choosing a formal equation. The essence of design procedure is shown in 
Fig. 1. Choice of the equation form and numerical values to describe the 
sediment-transport capacity will fix a relationship between the mean velocity 
of flow and the depth of flow as indicated by the heavy curve in the schema- 
tized diagram. Choice of the equation form and numerical values to describe 
the flow capacity will result in a family of curves relating velocity and depth 
for various canal slopes as indicated by the dashed curves. The dot-dash 
curve of constant tractive force describing the bank competence where it in- 
tersects the transport-capacity curve will specify one of the family of flow- 
capacity curves. Velocity, depth, and slope are thereby fixed, and, of course, 
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the width can be obtained from the continuity principle, 

The interrelationship shown above constitutes the general application of 
sediment-transport mechanics to stable-canal design. Only a complete, com- 
prehensive, rigorous theory, however, could substitute for a full scale model 
canal. It is by observation of similar canals that the designer can assign the 
needed numerical values, and sediment-transport mechanics can aid in as- 
sessing the degree of similarity of different canals. The parameters of the 
various equations can be used as comparative criteria and the equations 
themselves as approximate scaling relationships. If the proposed and exist- 
ing conditions are treated as ratios, the discrepancies between absolute 
values as computed by the equations should be greatly reduced. 

Practically speaking, few if any canals are completely stable. In planning 
remedial and maintenance work the principles of sediment-transport 
mechanics can be used in the same manner as in design, the principal dif- 
ference being that a canal almost completely similar is then readily available 
for observation. 
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Table I. Similarity of Bed-Load and "Kennedy" Formulas 
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MEASURING EVAPOTRANSPIRATION FROM ATMOSPHERIC DATA 


George S. Benton,! A.M. ASCE and Jack Dominitz2 
(Proc. Paper 1035) 


SYNOPSIS 


A method of evaluating evapotranspiration is presented which is based on 
the mass balance of water vapor in the atmosphere. The method is applied 
to various regions of the North American continent for the calendar year 
1949 and the results are compared with hydrologic data and with an empirical 
method of estimating evapotranspiration. 

It is shown that although accurate values of evapotranspiration can be ob- 
tained for large areas such as the entire continent, the accuracy of the meth- 
od is reduced as the size of the area under consideration is decreased. The 
potentialities and limitations of the proposed method are evaluated. 


INTRODUCTION 


Evapotranspiration—the loss of water from natural surfaces by evapora- 
tion and transpiration—has resisted measurement more strongly and is con- 
sequently less well known than most hydrologic parameters. The addition of 
moisture to the land by precipitation and its depletion through surface runoff 


. may be measured, but the rate at which the soil moisture supply is lost to the 


atmosphere under natural conditions cannot be measured directly by existing 
instruments. 

Yet, a quantitative understanding of the process of evapotranspiration is of 
considerable importance to the meteorologist and the hydrologist. The mete- 
orologist would like more information concerning the water balance of the at- 
mosphere, as well as the turbulent exchange processes that take place in the 
microclimatic layer. On the other hand, the hydrologist would like to be able 
to maintain a balance sheet of soil moisture content to assist him in evaluat- 
ing drought hazard or irrigation requirements, or to improve methods of 


river forecasting or reservoir management. Evapotranspiration is also of 
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importance to the military engineer, for off-the-road trafficability is strong- 
ly affected by the water content of the soil. 

Since instruments do not yet exist for measuring actual evapotranspiration 
under field conditions, scientists have turned to expressing evapotranspira- 
tion from natural regions by various crude formulae. Many such formulae 
have been developed. Some are based on simplified physical theory (e.g., 
Richardson and Montgomery(9))3 but most are purely empirical (e.g., 
Meyer;(7) Lloyd(6)), Naturally, one of the greatest difficulties in devising 
and using such formulae is the lack of objective evapotranspiration measure- 
ments to use as control data. Tank estimates of evapotranspiration and soil 
moisture determinations in experimental areas provide the bulk of such data, 
although in many regions the mean annual evapotranspiration can be computed 
with reasonable accuracy from precipitation and runoff measurements. 

A recent empirical formula for evapotranspiration which has received 
widespread attention was devised by Thornthwaite.(11) Thornthwaite devel- 
oped a formula for potential evapotranspiration, which he defined as the maxi- 
mum amount of water that could be lost from a region under a given set of 
climatic conditions. He also developed a simple water budget for estimating 
actual evapotranspiration from potential evapotranspiration. Thornthwaite’s 
method generally yields estimates of evapotranspiration which are reason- 
able in their seasonal distribution and annual total (Carnahan and Benton). (3) 
However, such empirical estimates cannot be properly evaluated until accur- 
ate measurements of evapotranspiration become available for natural re- 
gions. 

One possible method of determining evapotranspiration, which has not 
heretofore been adequately explored, involves an evaluation of the atmos- 
pheric water balance over the continents. The latter subject has been dis- 
cussed by many authors, but little quantitative information was available until 
Holzman(5) presented his analysis of the relationship between the air mass 
and the hydrologic cycles. The concepts advanced in this paper were later 
extended and developed by Benton, Blackburn, and Snead.(1) In recent years 
our understanding of the atmospheric water balance has been greatly in- 
creased by the extensive use of electronic computers in meteorological re- 
search. As an outgrowth of this program, the flow of water vapor over the 
North American continent was analyzed and discussed by Benton and 
Estoque. (2) 

The availability of extensive meteorological data and the ability to process 
this information numerically in an efficient manner, led to a growing realiza- 
tion that our knowledge of atmospheric vapor flow might be sufficient to allow 
the direct computation of evapotranspiration from a statement of the mass 
balance of atmospheric water, developed in a manner analogous to the fa- 
miliar hydrologic equation. Possibly the first attempt to evaluate evapotrans- 
piration in this manner was presented by Moller, (8) who used a sparse net- 
work of meteorological stations over northwest Europe. A more extensive 
study of evapotranspiration over the North American continent was presented 
by Benton and Estoque. (2) 

The investigations of Benton and Estoque showed that for land areas of 
continental extent, estimates of evapotranspiration obtained from atmospheric 
data are at least as accurate as those which can be obtained from hydrologic 
estimates or from empirical relationships. The major objective of the 
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present study is to evaluate the extent to which atmospheric data can be suc- 
cessfully used to determine evapotranspiration from smaller areas. The at- 
mospheric balance equation is derived in the following section, and applica- 
tions of this equation to specific areas of various sizes are then presented. 
Computed evapotranspiration is compared with estimates obtained using 
Thornthwaite’s empirical procedure, as well as with annual evapotranspira- 
tion totals obtained from precipitation and runoff records. 

It is shown that although the procedure for evaluating evapotranspiration 
from atmospheric data is physically sound, practical difficulties arise. 
These difficulties, related to problems of statistical sampling, are compara- 
tively unimportant when|large areas and long time intervals are considered 
but become more critical as the size of the area and time interval is reduced. 
Nevertheless, the atmospheric method offers a fundamentally simple ap- 
proach to the evaluation of evapotranspiration which may become of increas- 
ing practical value as meteorological observations and techniques of elec- 
tronic data processing are improved. 


Development of the Atmospheric Balance Equation 


The conventional equation, expressing the basic relationship between 
evapotranspiration and other hydrologic parameters for any time interval, is 
the hydrologic equation: 


P=E+R+AS (inches) (1) 


Here P, E, R, and AS are precipitation, evapotranspiration, net runoff, and 
change in storage of surface and sub-surface water all expressed in inches 
depth, respectively. Ground water flow into or out of the region can usually 
be neglected if a large area is considered. Unfortunately, equation (1) cannot 
be solved for E for periods less than a few years in length, since AS is not 
negligible for short time periods. 

A similar mass balance equation may be developed for the atmosphere, 
expressing the relationship between evapotranspiration and meteorological 
parameters. The following expression is obtained for an arbitrary time in- 


‘terval and region: 


meee nN (gm/em*) (2) 


Here F is the net outflow of atmospheric water vapor (F = outflow minus in- 
flow); E is evapotranspiration; P is precipitation; and AN is the change in the 
precipitable water content of the atmosphere. Two terms, AN and F, require 
further explanation. 

The change in precipitable water content, 4N, is defined as the difference 
between the atmospheric water content at the end of the time period and that 
at the beginning. It is analogous to the change in storage term in the hydro- 
logic equation, although fortunately the moisture holding capacity of the at- 
mosphere is so low that it has none of the numerical importance of the latter 


‘term, Precipitable water, N, is defined as 


N= f Be dp (gm/cm”) (3) 
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Here oe is the specific humidity, in gm/kg; p is the pressure in millibars 


which varies from zero to its surface value, pg; and g is the acceleration of 
gravity, in cm/sec2, Over an arbitrary time interval T : 


AN = Nj-N, (gm/em) (4) 


where Ng is the precipitable water at the end of the interval, and Nj is the 
precipitable water at the beginning. 

The net outflow of water, F, from the arbitrary area is more difficult to 
evaluate. Let the vector w be the mean transfer of water vapor at an arbi- 
trary point in the atmosphere. Then: 


we ae (2) Cc dt (gm/cm-mb-sec) (5) 


Here c is the horizontal wind vector, in cm/sec, expressed as a function of 
time t. Note that w is the flux of vapor in gm/sec through a cross-sectional 
area one centimeter in width and one millibar in thickness. The mean trans- 
fer of water vapor integrated vertically through the atmosphere is: 


W= A ‘wp (gm/cm-sec) (6) 
noe i 


where pj, is now the mean surface pressure. Finally the net outflow of water 
vapor, F, from an arbitrary area A is given by: 


F=FpA-wdh — (em/em?) (2) 


where n is a vector of unit length normal to an infinitely small segment of 
the perimenter, dA, and pointing out of the area; and the line integral is 
evaluated around the perimeter of the area, Note that the transfer of liquid 
water is neglected. It can be shown that the resulting error is negligible in 
all save extremely unusual local situations. 

Substituting equation (7) into (2): 


E=P+AN+— $A-Wa) (gm/em”) (8) 


All of the variables on the right hand side of equation (8) can be evaluated 
provided a rawinsonde or radiosonde network is available to provide esti- 
mates of wind velocity and humidity through the atmosphere. Such a network 
is operated by U. S, Weather Bureau and Armed Services personnel, and par- 
ticipating weather stations are scattered throughout the United States and 
North America, Equation (8) thus provides a physically sound method of 
computing evapotranspiration from natural areas. Furthermore, if it is 
possible to evaluate the various terms in the atmospheric balance equation, 
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change in storage of moisture in the soil could also be determined by 
uvstituting the appropriate value for evapotranspiration in the hydrologic 
quation (1). 

It must be pointed out that although equation (8) is almost exact, and in- 
lves no hazardous assumptions, other difficulties do arise. The net inflow 
r outflow of water vapor through the atmosphere is for most areas and time 
tervals a small difference between two large numbers, which must be de- 
rmined with considerable accuracy if satisfactory results are to be ob- 
ined. This limits the size of the area and the length of the time interval 
hich can efficiently be considered. The following sections describe the ap- 
lication of equation (8) to actual data. 3 


Description of Data 


The calendar year 1949 was chosen for analysis because the necessary 
ee topical data were readily obtainable from the research project with 
vhich this study was associated. Monthly evapotranspiration was evaluated 
sing equation (8) for each of four sub-continental areas, represented by four 
olygons with meteorological stations at each vertex. These polygons were 
jhosen to correspond approximately to the following regions: Canada, the 
jentral United States, north-central United States, and the Ohio River Basin. 
‘wenty-two radiosonde stations were utilized in the study, as listed in 
fable 1. At each station measurements of the elevation of constant pressure 
jurfaces and of humidity were recorded twice a day, 0300Z and 1500Z (i.e., 
areenwich Civil Time), at each of three levels in the atmosphere: 850, 700, 
nd 500 mb. The resulting data were entered on punchcards. 

Geostrophic winds were used in evaluating the mean transfer ofswater 
apor, W, at each station. The geostrophic wind is a theoretical wind which 
an be computed from the slope of constant pressure surfaces even in the ab- 
ence of actual wind observations. Geostrophic winds are comparatively un- 
ffected by local topography and by small scale turbulence since they are 
omputed from the average slope of the pressure surface over a large region, 
ften several hundred miles in width. For these and other reasons, they are 
enerally used in computing atmospheric transfers of energy, momentum, 
nd other physical quantities. In the present study, geostrophic winds were 
omputed graphically from maps routinely prepared by the Weather Bureau- 
ir Force-Navy Analysis Center in Washington, D. C. These winds were en- 
sred into the punch-cards along with pressure and humidity data. 

Once the basic atmospheric data were entered into the punch-cards, the 
lowing extensions were performed. For each observation, the eastward 
nd northward components of the geostrophic wind velocity were evaluated 
nd each was multiplied by the water vapor content. These two components 
f water vapor transfer were summarized and averaged for each station- 
ionth at each level to give w (equation (5)), and the integrated water vapor 
‘ansfer at each station, W, was then evaluated graphically from the com- 
onent relationship: 


Po P, 
w=[ w dp =1/, ap+}/ wy, dp  (gm/cm-sec) _— (9) 
‘0 ° fo} 


se 
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RADIOSONDE STATIONS USED IN THIS STUDY 


Station 


Miami, Fla. 

New Orleans, La. 
Brownsville, Texas 
El Paso, Texas 
Little Rock, Ark. 
Atlanta, Georgia 
Greensboro, N. C. 
Portland, Me. 
Buffalo, N. Y. 

St. Cloud, Minn, 
Lander, Wyo. 
Medford, Ore. 

Santa Maria, Calif. 
Port. Hardy, .0.. G0. 
Edmonton, Alta. 
Churchill, Man. 
Stephenville, Newfoundland 


Upper Frobisher, Baffin Is. 


Yakutat, Alaska 

McGrath, Alaska 

Barrow, Alaska 

Coppermine, Dist. of Mack. 


TABLE 1 


August, 195 


Lat. 
North 


25° 49' 


30 


00 
55 
48 
44 
39 
05 
39 
56 
35 
48 
23 
56 
4l 
34 
21 
32 
44 


=) 


58 
20 
47 


Lor 
Wes 


80° 
90 
97 
106 
92 
84 
79 
70 
78 
94 
108 
ize 
120 
127 
113 
94 
58 
68 
139 
155 
156 
115 
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here x and y are directed towards the east and north respectively. Vapor 
transfer above 400 mb was neglected. One possible error in this evaluation 
wolves the layer between 850 mb (the lowest level of observation) and the 
surface. Previous work has shown that the maximum flux of vapor may oc- 
ij in this layer. (2) 
As a last step in computing the net vapor flux out of each region, the ap- 

— line integral of the component of W normal to the boundary of each 
olygon was evaluated as indicated in equation (7). A linear variation of the 
tegrated water vapor transfer was assumed between stations on the perime- 
er, although such an assumption can lead locally to non-negligible errors. 
finally, the remaining quantities required for the determination of evapo- 
ranspiration were computed. The precipitable water content of the atmos- 
here at the beginning and end of each station month was estimated graphical- 
y in a manner similar to that of vapor transfer, and the mean change in 
recipitable water for each region was then found by averaging weighted sta- 
ion values using the method developed by Thiessen.(10) Thiessen’s method 
as also used in evaluating precipitation from station measurements. With 
yet vapor flux, change in precipitable water, and mean precipitation deter- 

ined, evapotranspiration was computed for each month for each of the four 
egions under consideration using the atmospheric balance equation. 
__ As indicated previously, evapotranspiration was also estimated by two 
pther methods. The first of these was Thornthwaite’s empirical equation. 
Although Gilbert and Van Bavel (4) have shown that Thornthwaite’s equation 
ields somewhat better results when daily temperatures and precipitation 

re used, the present study was limited for simplicity to the use of monthly 
ata in accordance with the procedure outlined in the original reference. (11) 
hirty-seven stations scattered through North America, were employed in 
yrder to get adequate coverage. Thiessen’s method was again used to find 
lhe mean value of evapotranspiration from station values. 
: A further check on the values of the total evapotranspiration for any year 
may be obtained by evaluating the runoff, neglecting the change in storage, 
\fter equation (1). It is anticipated that such estimates would usually be with- 
n two inches of the real values.4 River discharge data were compiled from 
J. S. Geological Survey records.(12) Since only United States records were 
‘eadily available, three polygons (central United States, north-central United 
‘tates, and the Ohio River Basin) were checked. This involved the examina- 
ion of sixty natural watersheds located within these polygons. 


Presentation of Results 
North American Continent 


As noted previously, the field of water vapor flux over the North American 
ontinent has been described in detail by Benton and Estoque. (2) As part of 
heir study, the water balance of the continent was computed and measure- 
aents of evapotranspiration using atmospheric data were compared with hy- 
rologic estimates and with precipitation. For the convenience of the reader, 
he monthly distribution of evapotranspiration for the continent is reproduced 


. Calculations of the change in storage for the year 1949 for ten sample 
stations, using a method developed by Thornthwaite,(11) indicates that 
two inches is usually a conservative estimate for large areas. 
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in Figure 1. For a discussion, the reader is referred to the original refer- | 
ence. 


Canada 


The distribution of the monthly evapotranspiration for the area as illus- 
trated in Figure 2 is shown in Figure 3. A sample tabulation of the various 
quantities in equation (8) is shown in Table 2. The relationship between the 
evapotranspiration evaluated from atmospheric data and that determined from 
Thornthwaite’s empirical expression is very similar to that found for the 
North American continent. Values computed by use of the atmospheric bal- 
ance method are higher in the winter; reach a maximum later in the year; 
and result in a larger annual total. These differences may in fact be caused 
by certain assumptions made in the empirical formula. Thus, Thornthwaite 
assumed that evapotranspiration ceases whenever the mean monthly tempera- 
ture falls below freezing, and that potential evapotranspiration is fully real- 
ized as long as any water remains in storage in the soil. Neither of these 
assumptions is realistic. They would result in errors similar to the differ- 
ence found in Figure 3. 

Due to the unavailability of river discharge data for most of northern 
Canada, it is difficult to estimate the annual evapotranspiration using runoff 
measurements. It may be noted, however, that the annual evapotranspiration 
determined from atmospheric transfers of water vapor was 17.0 inches, 
whereas precipitation totalled 21.7 inches. This would require runoff for the 
year to total only 4.7 inches, if change in storage were neglected. Also, an 
increase in storage might be expected since precipitation in the last few 
months of 1949 exceeded that in corresponding months of 1948. Assuming the 
value of the precipitation to be correct, it becomes evident that evapotrans- 
piration may be slightly overestimated. This apparent discrepancy and its 
cause will be investigated in a later section. 


Central United States 


The monthly distribution of evapotranspiration for this area (Figure 2) is 
shown in Figure 4. In general, the relationship between the curves obtained 
from the present method and from Thornthwaite’s empirical formula remains 
the same, although the atmospheric balance gives consistently higher monthly 
values, A check on the annual total may be made since the runoff is known. 
For the year, precipitation amounted to 38.0 inches and runoff to 10.9 inches. 
Annual evapotranspiration should therefore be 27.1 inches plus a correction 
due to changes in storage between the beginning and end of 1949, This change 
was probably negative for this area primarily because November, 1949, was 
a very dry month (U. S. Weather Bureau).(13) A reasonable estimate of 
evapotranspiration might therefore be between 27 and 29 inches for the year. 
The computed value of 34.6 is high, and cannot be explained by inadequacies — 
in the precipitation and runoff data. Nevertheless, the seasonal distribution 
of evapotranspiration is good. 

It may be noted in passing that the empirical method yields an estimate of 
the annual evapotranspiration equal to 24.7 inches. This figure is apparently 
too low. The smooth curves in Figure 4 should not, therefore, be regarded 
as accurate. The correct evapotranspiration probably lies somewhere be- 
tween the two sets of estimates. 
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North-central United States 


The area chosen is presented in Figure 5 and the computed monthly dis- 
tribution of evapotranspiration is compared with that of the empirical method 
in Figure 6. Although both curves display the same general characteristics, 
the difference found in the examples presented above are again apparent. 
Thus, the empirical method yields lower values and an earlier maximum. 

The annual evapotranspiration may be roughly estimated from the hydro- 
logic equation. Precipitation for this area was found to be 44.8 inches while 
runoff was 16.2 inches. This results in an annual total of approximately 29 
inches plus a correction due to change in storage. It is estimated that this 
change in storage was negative since a large portion of this polygon (the Ohio 
River Basin) experienced a fairly large negative change.(3) The annual 
evapotranspiration may, therefore, be between 29 and 31 inches. The com- 
puted value of 32.8 inches may be slightly too high; however, it is probably 
more accurate than the low estimate obtained by use of Thornthwaite’s 
method (25.4 inches). It is evident that although the individual monthly values 
of evapotranspiration obtained from the atmospheric moisture balance may 
not be too precise, the seasonal distribution and the annual total are reason- 
ably adequate. 


Ohio River Basin 


A volygon representing the Ohio River Basin (Figure 7) is next considered. 
This area was chosen primarily because much pertinent work had been done 
by Carnahan and Benton(3) in evaluating hydrologic parameters, including 
evapotranspiration, for this area and year. The comparison between the two 
methods (atmospheric transfer and Thornthwaite) is presented in. Figure 8. 
The resemblance between the two sets of data is poor and the trend of the 
values obtained from the proposed method during the first six months seems 
to be opposite to what is expected. The maximum, although of the proper 
magnitude, occurs late in the summer season. The annual evapotranspiration, 
however, seems to be of the proper magnitude, 28.5 inches, compared to the 
“best estimate” of 28.2 inches, found by Carnahan and Benton employing many 
different methods. Therefore, although the monthly data and seasonal distri- 
bution are unsatisfactory the annual value seems to be about correct. 

Whether the discrepancies in the monthly data are caused by the fact that this 
is the smallest of all areas considered (826,000 km2), will be discussed in a 
later section. It is interesting to note, however, that this region is not very 
different in location or size from the north-central United States, for which 
much better results were obtained. This might suggest that some systematic 
error, rather than a random sampling error, is involved. 


Statistical Analysis of the Atmospheric Balance Equation 


The atmospheric balance method as developed earlier, is based on sound 
physical principles and theoretically should give correct values for the rate 
of evapotranspiration from any selected area. From the limited number of 
examples furnished in the previous section, however, it is evident that cer- 
tain difficulties are encountered in practice which must be removed before 
this method can be applied to sub-continental regions. To clarify these dif- 
ficulties, statistical analysis of an idealized system has been used to deter- 
mine the factors affecting the accuracy of the atmospheric balance equation. 


BENTON - DOMINITZ 1035-15 


ASCE 


IOJ saze1S psituQ [e1}Uua0-Yy210U ay JUaSaada1 03 pasn eaIy 


‘cg ‘bry 


S3LVLS Q3LINN 
TWYLN30- HLYON 


August, 1956 


HY 4 


~yq10U ay} OJ uoTJeITdsuerjoOdeas Jo uot Ingraistp ATy.uOM ayT *9 “Bry 


J30 AON 190 1d3S SNV INC NAC AVW YdV YVW 834 we 


eee 


ay 
\ 
\ 


Pa es a 

atl soe 
ame 
= 


SHIOMYJUJOY] 19457 
DJDQ juaseld 


ase 
eras 


© 
SSHONI NI NOILVYIGSNVYLOdVAS 


1035-17 


BENTON - DOMINITZ 


ASCE 


“uot e1Tdsuerj0dead 


Buyienyeas OJ urseg IaATY OTYO 241 JUasaadex 02 pasn eoay *) “Bry 


NISVG Y3AIY 


HY 4 August, 1956 


1035-18 


OTYO 942 OF udT{JeAardsuezzOdeaa Fo wot nqraIstp ATYyIUOW aYT “9 “Brg 


930 AON 190 1d3S 9NV INC NAP AVW YdV YW 934 Neo 


Y 


OILVYIdSNV 


© 
An 


SSHONI NIN 


ASCE BENTON - DOMINITZ 1035-19 


Such analysis was limited to the divergence term since the possibility of er- 

ror is greatest in evaluating the net flux of water vapor and all errors in 

_ computing the other terms are negligible in comparison. 

Consider an idealized area, A, on whose perimeter, P', n stations are 

_ located with equal spacing, Each station reports N observations at 7 levels, 
each representing equal pressure increments through the atmosphere. Such 
a model is illustrated in Figure 9. 

Let the standard error of the mean value of any component of w (the mean 
atmospheric transfer of water vapor at any level and station) be given by s. 

Integrating the mean transfer through the height of the atmosphere: 


W= . wh, (gm/cm-sec) (10) 
isl 


where hj is the height, in millibars, of the atmospheric layer represented by 
the value of the mean transfer, wj. Since all levels are spaced equally, hj is 
a constant and will hereafter be denoted as h. 

_ The standard error of estimate of any component of W (the integrated va- 
por transfer) is made up of two parts: error due to sampling, since only a 
finite number of levels are used through the atmosphere; and error due to 
inaccuracies in the computed values of w. For averages over a moderate 
period of time (e.g., a week or more), w varies smoothly through the atmos- 
phere; and hence if several levels are used the representative sampling error 
will be very small. Neglecting the sampling error, and assuming for sim- 
plicity that s is the same at each level, the standard error of any component 
of W becomes: 


| r=sh/7 (gm/cm-sec) (11) 


| After converting the divergence of water vapor over an area into a line in- 
tegral, and after assuming linear variation of W between stations on the 
| perimeter: 


n 
| Vew = Sr Sy fs (W, + Wt, (gm/cm”-sec) (12) 


where W..4= W, and A, is a unit vector normal to the line L, and directed 


' 1 
| outward from the polygon. Since the stations were equally spaced around the 
perimeter, Lj is a constant L. It will also be assumed that the standard er- 
ror of estimate of each component of W does not vary from station to station. 
The standard error of estimate of V°W is then given by: 


s6.= =e (gm/em*-sec) (13) 


Sampling error has again been neglected: i.e., the assumption of.linear vari- 
ation of W between stations has been assumed valid. 
Introducing the perimeter P', and the number of levels 7 from equation 
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: 
: 


Fig. 9, An idealized equilateral polygon used in evaluating errors 
in water vapor divergence, Vertical observations are for 


equal increments of pressure. 


CE BENTON - DOMINITZ 1035-21 
11), the standard error in V-W becomes: 
| id 
| 
| by eae | (gm/cm?-sec) (14) 
a7 


1 ay H is the depth of the moisture bearing atmosphere. The 0.95 confi- 
ence interval for V-W is given by: 


Vw se shar (gm/cm?-sec) (15) 


Although this formula was computed for an ideal case, and certainly not 
or any of the polygons treated in the text, valuable information may, never- 
heless, be obtained from it. Thus it can be seen that the accuracy of the 
vergence term may be increased by employing a greater number of stations, 
evels, and observations. It is also evident that the size and shape of the area 
o play an important role in obtaining a true estimate of the divergence. An 
timum region would be a large area whose shape approximates a circle. 
__ Let us now consider the problem of obtaining a numerical value of the 
onfidence interval from expression (15). The only unknown term is s, the 
tandard error of estimate of the components of w (the mean transfer of 
ater vapor at any station-level). It is obvious that s includes three types of 
rror: representative sampling error, random instrumentation error, and 
-onsistent error. The last type of error is comparatively unimportant, since 
ts percentage effect will not increase upon differentiation. Thus we may 
mite: 


s &/s*+ me (gm/cm-mb-sec) (16) 


i ie S 1 is the contribution due to representative sampling error, and so is 
contribution due to random instrumentation error. 

By estimating the average standard deviation of wind speed and mixing 
ratio within a twelve hour period (the interval between successive observa- 
ions) it was found that a reasonable figure for sy was 0.7/ /N gm/cm-sec. 
Also, if the observational error in measuring wind speed is about 1.5 m/sec 
ind in measuring mixing ratio is about 0.3 gm/kg, sg will be approximately 
1.8//N gm/cm-mb-sec. Thus the value of s is about 1.1//N. 

Subject to the stringent assumptions inherent in the derivation, the confi- 
lence interval for V-W thus becomes:® 

: 


VeWt or Ww a (gm cm*-sec) (17) 
| 


vyhere N is the number of observations at each station-level. 


. For example, it must be noted that expression (17) would give a confidence 
interval equal to zero if an infinite number of stations, levels, or observa- 
tions were employed. It is clear, therefore, that expression (17) is only 
_approximate as a result of the idealization of the problem and must be in- 

_terpreted carefully. : 


‘@ 
‘ 
¢ 
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The Accuracy of Evapotranspiration Estimates 


In the present section, the confidence interval for divergence of water 
vapor has been computed for four sub-continental regions and for North 
America using equation (17), The data are presented in Table 3. It is evi- 
dent that the P/A ratio is the dominant factor. As an approximate rule, the 
accuracy of the method varies inversely as the square root of the area. 

Table 4 compares the values of the annual evapotranspiration obtained by 
the proposed method with those obtained from other methods. In four out of 
eight comparisons the confidence interval does not include the estimate de- 
rived by use of other methods. Although this is not conclusive, since the 
other methods are themselves subject to considerable question, the possibil 
ty must be considered that the confidence intervals are too low. This is 
likely, since the assumptions made in deriving equation (1 7) were not alway 
realistic. 

In particular, two assumptions are open to question. First, it was as- 
sumed that there were no sampling errors introduced due to interpolations 
between levels or in extrapolation to the surface. However, water vapor 
transfers were computed only at 500, 700, and 850 mb, and as has been note 
earlier the extrapolation to lower levels is very critical, since maximum 
transfers sometimes occur below the 850 mb level. Although every effort 
has been made to achieve accurate extrapolations, there is no assurance tha 
the results are completely satisfactory. 

A second assumption which may often be invalid is that the variation of 
water vapor transfer from one station to another can always be adequately 
described by a linear function. Assuming the average spacing of stations to 
be more or less constant, errors due to non-linearity would tend to increas 
as the size of area (and hence the number of stations) is decreased. It is en 
tirely possible that errors of this type also appear in Table 4, 

It is of interest to note that inaccuracies in the evaluation of the atmos- 
pheric balance equation which arise due to lack of adequate representative 
sampling in the vertical and around the periphery of the area can largely be 
eliminated by use of more levels and a larger number of radiosonde stations 
The improvement in evapotranspiration estimates to be obtained by the use 
of more extensive meteorological data has yet to be determined by actual 
computation. However, it is entirely possible that the success achieved for 
the North American continent can be extended to areas as small as the Ohio— 
River Basin and to time intervals as short as ten days. 


CONCLUSIONS 


This study has presented a method of evaluating evapotranspiration based 
on sound physical principles. The method has been applied to various sub- — 
regions of the North American continent; and the estimates obtained have 
been compared with hydrologic data where available, and with Thornthwaite’s 


are considered. However, the technique may become unsatisfactory when 
Substantially smaller areas (e.g., of the size of the Ohio River Basin) are 


One obvious possibility is that this decreased accuracy for smaller areas 
is related to statistical sampling, and the fact that the divergence of water 
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vapor represents a small difference between two large numbers (the total in- 
flow and outflow). Statistical analysis of an idealized system leads to a 
formula for the confidence interval for evapotranspiration which verifies that 
the error shouid increase as the area is decreased. However, numerical 
evaluation indicates that for smaller areas this estimate of the confidence 
interval is probably too small. This may, of course, be due to inadequacies 
in the expression for the confidence interval. For example, it was assumed 
that errors in extropolating moisture transfer from 850 mb to the surface 
and in assuming linear variation of transfer between stations were negligible. 

Despite the difficulties encountered, the potentialities of this method can- 
not be discounted. Improved accuracy can surely be obtained by use of more 
meteorological stations and more levels of observations. The use of ob- 
served winds for smaller areas is also a possibility, although the adequacy 
of moisture divergence fields obtained from such transfers is open to ques- 
tion. In any case, the atmospheric balance method has yielded monthly esti- 
mates of evapotranspiration for large areas which are at least as accurate 
as those obtained from any other known method. With further refinements 


_ of meteorological observations and techniques of data processing, this 


method of measuring evapotranspiration may prove to be of considerable 
practical interest to the engineer. 
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SYNOPSIS 


Chemical coagulation is a basic water treatment process, which prepares 
the water for subsequent treatment operations. Coagulation includes addition 
of chemicals, their distribution, floc formation and floc build-up. The con- 

ditioning process that leads to floc formation and build-up is “flocculation.” 
Flocculation is usually accomplished by prolonged agitation at low velocity. 

Commonly used coagulants are alum and ferric sulfate. Alum’s ionization, 
hydrolysis, and reaction with alkalinity are of interest but the process of 
coagulation cannot be explained on the basis of simple chemical relationships. 
Colloidal particles play an important role in coagulation. Some of the signifi- 
cant characteristics of colloidal dispersions are large specific surface area, 
high adsorptive capacity, and particle charge. Colloidal particles commonly 
encountered in natural waters carry a negative charge, whereas those formed 
during the initial phase of the coagulation process often have a plus charge. 
Ions possessing a charge opposite to those borne by the colloidal particles 
tend to precipitate the latter. An ion’s precipitating power increases tre- 
mendously with its valence. 

Important factors in the coagulation process are pH, dissolved solids, 
‘temperature, time-concentration, coagulant aids, exchange capacity, and 
mixing. The last named is essential for efficient coagulation and consists of 
two stages: rapid mixing for a short period followed by slow mixing for a 
much longer time. Flocculation is accomplished in the latter mixing stage. 

Laboratory studies of the coagulation process are essential. Flexibility 
is an important aspect of plant design. 


INTRODUCTION 


Clarity is one of several highly-regarded characteristics of water. This 


: Note: Discussion open until January 1, 1957. Paper 1036 is part of the copyrighted 
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is especially true if the water is to be used for municipal purposes. Some 
industries also require process water virtually free of turbidity. The treat- 
ment of waste water—sewage and industrial wastes—involves concentration 
and separation of suspended matter, that is clarification of the waste water. 

One of the oldest methods of water treatment is the removal of suspended 
solids by subsidence.(1) However, it is a well known fact that the nature of 
the suspended matter in most natural waters is such that subsidence, un- 
assisted, will not, even in combination with rapid sand filtration, render the 
water suitable for general use. Chemical coagulation, a water treatment key- 
stone, is an almost universal requirement. Coagulation prepares the water 
for sedimentation and greatly increases the efficiency of the latter. Coagula- 
tion plus sedimentation prepare the water for filtration. Coagulation ranks 
high on the list of things that must be done to render water safe, clear and 
palatable. 

Much has been learned about the whole process of coagulation especially 
in the last 25 to 30 years; yet, today the practice of coagulation has certain 
of the empirical characteristics of an art. The process is complex; it is 
sensitive; and natural waters are endowed with rugged individuality. Weare 
probably just beginning to obtain a fairly clear picture of a process that has 
been practiced for centuries. There is still room for exploration. 


Coagulation - Flocculation 


The coagulation process includes “addition of chemicals to the water, 
their uniform distribution through the water, and the building up of a readily 
settleable floc, usually by prolonged agitation”.(2) The process of condition- 
ing that leads to floc formation and build-up is termed “flocculation,” the 
subject of this discussion. Obviously it is difficult to separate “flocculation” 
from the overall process of “coagulation.” Any discussion of the former can- 
not help but include the latter. Therefore much of this discussion relates to 
the process of coagulation which includes flocculation. 

The aim of flocculation is the promotion of coalescence of the precipitates 
derived from the coagulating chemicals and the natural suspended matter into 
fast settling precipitates termed “flocs” in water treatment parlance. As will 
be mentioned subsequently in this paper “conditioning” or “promotion-of- 
coalescence” is commonly accomplished by slow mixing by hydraulic or, 
more frequently, mechanical means. An important step in the coagulation 


process, flocculation appears to merit the attention it has received especially 
in recent years, 


HISTORICAL NOTES 


Mankind because of his need for water and his curiosity regarding it has, 
through the ages, exhibited much interest in its qualities and the means by 
which they might be altered the better to serve his current needs. Conse- 


quently it is not surprising to note that coagulation of water was practiced in 


ancient times. Apparently there was appreciation of the need for mixing as 
early as the sixteenth century for an Italian physician who visited Egypt then 
described the use of almonds for coagulation as follows: 

“, .. they smeared the edges of the vessel with five sweet almonds, prop- 
erly crushed, and grasping the almonds in the hand suddenly plunged hand 
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and arm into the water up to the elbow and moved elbow and hand vigorously 
and violently this way and that through the water stirring it up until they had 
made it far more turbid than before. Then they withdrew the arm from the 
vessel leaving the almonds in the water. They let it clarify and it was prop- 
erly clear in three hours time.”(1) 

In the United States coagulation was used in connection with a New Jersey 
public water supply in 1885 and studies of aluminum sulfate as a coagulant 
were undertaken at Rutgers University during the same year.(1) Since 1885 
coagulation “has been, with few exceptions, part and parcel of rapid filtration 
in America.”(1) 


Some Chemical Aspects of Coagulation 


Chemicals(3)(17) 


Numerous chemicals have been considered for use as coagulants. How- 
ever, relatively few are in common use today. These are listed in the follow- 
ing table. 

TABLE I 


CHEMICAL COAGULANTS 


Name Approximate Formula 
Aluminum Sulfate Alg(SO4)3 15H,O 
(Alum, Filter Alum) 

Ferric Sulfate Fe9(SO4)39H20 
(Ferrifloc, Ferrisul) 
Ferric Chloride FeClg or FeCl36H20 
Chlorinated Copperas FeCl3 and Fe2(SO4)3 
Ferrous Sulfate FeSO47H20 
(Copperas) 
Lime 

Hydrated lime Ca(OH)9 

Quicklime CaO 
Activated Silica SiO, 


An important characteristic of aluminum and iron salts is their ability, 
under appropriate conditions, to form gelatinous-appearing precipitates or 
“flocs.” A more precise term would be “chemical flocs” to distinguish them 
from “biological flocs” produced in sewage treatment processes. Of the 
chemicals listed in Table I aluminum sulfate (alum) is probably the most 
widely used, It is cheap and generally available. Lime is frequently applied 
merely as a coagulant aid but in some waters it alone can produce effective 
coagulation. Lime is a widely used softening and coagulating agent for the 
hard waters encountered throughout much of the mid-western section of the 
United States. Consequently its inclusion in the list of commonly used coagu- 
lants is believed justified. Activated silica is a newcomer and is generally 


1036-4 HY 4 August, 1956 


considered a coagulant aid. However, at least one authority (3) has classed it 
as a coagulant. 


Chemical Reactions 


Although alum has certain inherent chemical peculiarities it is widely 
used and therefore alum has been selected to illustrate certain coagulation 
reactions. The commercial product is approximately Alg(SO4)g3 15H20. 
(Detailed specifications have been prepared by the American Water Works 
Association).(4) When alum is dissolved in water it ionizes yielding alumi- 
num ions (Al+++) bearing a high positive charge and negatively charged 
(SO4-~) ions. This may be illustrated by the equation 


Aly (SO4)3 coed 2Alt+t++ + 380477 
Alum also hydrolyzes readily as follows: 
Alg(SO4)3 + 6HgO — 2A1(OH)3 + 3H2SO4 


In most natural waters the sulfuric acid formed by this reaction is prompt- 
ly neutralized by the natural alkalinity of the water. (This is assumed to be 
calcium bicarbonate.) 


3Ca(HCO3)9 + 3HpSO04 — 3CaSO4 + 6COg + 6H2O 


(The acid may also be neutralized by the addition of an alkaline agent, com- 
monly lime.) 

The overall equation, assuming 15 molecules of water of crystallization, is 
often written as follows: 


Alg(SO4)3 15 H2O + 3Ca(HCO3)2 —~ 2A1(OH)3 + 3CaSO4 + 6CO9 + 15H2O 


Aluminum hydroxide (more accurately the hydrous oxide of aluminum) is 
an amphoteric substance. This means that under acid conditions (low pH) it 
behaves as a base and under alkaline conditions it behaves as an acid. Con- 
sequently it will be dissolved by too great a surplus of either H+ or OH™ ions. 


Al(OH)3 + 3H* — Alttt + 3H\O 
(Solution at low pH) 

AL(OH)3 +OH™ —+ AlO)” + 2H,0 
(Solution at high pH) 


Precipitates formed by iron coagulants do not exhibit amphoteric behavior, 
at least not the extent shown by those derived from alum. 

These reactions are intended to serve as simplified illustrations and do 
not purport to indicate the exact chemistry of the coagulation processes nor 
the composition of the alum floc.(5) The chemistry involved is much more 
complex than these simple equations suggest. However, certain observa- 
tions of a qualitative nature may be made based on a consideration of these 
equations. There are: 


a) Alum yields both plus and minus ions of fairly high charge. 

b) Alum produces a precipitate. (The composition varies). (5) 

c) Alum is an acid salt. Its hydrolysis can result in a reduction in the 
natural alkalinity of the water and an increase in the carbon dioxide content. 
The result is an increase in the hydrogen ion (H+) concentration or a reduc- 
tion in the pH value of the water. 
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d) Because of its amphoteric nature the solubility of alum floc is sensitive 
to the pH value of the water. pH adjustments may be required to achieve op- 
timum conditions for precipitation. 


Similar equations leading to generally similar conclusions could be written 
for most of the iron coagulants (ferric sulfate, ferric chloride, chlorinated 
copperas) Ferrous sulfate (copperas) behaves similarly; however, its reac- 
tion is further complicated by the necessity of oxidizing the iron from the 
“ferrous” to the less soluble “ferric” form. As previously noted an impor- 
tant difference between alum and the iron coagulants is the fact that the pre- 
cipitates or flocs formed by the latter are not amphoteric. 


Colloidal Particles 


The coagulation process cannot be described in terms of simple, stoichio- 
metric relationships. Coagulation is intimately associated with the science 
of colloidal materials. Toa great extent understanding of the mechanism of 
coagulation is dependent upon knowledge of the characteristics and behavior 
of matter in the colloidal state. Colloid chemistry has been characterized as 
“the twilight between chemistry and physics.” (6) 

Colloidal dispersions may be thought of as consisting of a homogeneous 
medium such as water and particles dispersed therein. The term “colloid” 
does not refer to a particular kind of matter but to a somewhat indefinite 
state of subdivision. In general colloidal particles are not microscopically 
resolvable nor can they be readily separated from water by filtration methods. 
The author is not aware of any satisfactory, short definition of colloidal sys- 
tems. No true dividing line exists between what may be termed suspensions, 
colloids, and solutions. Colloidal dispersions have been referred to as 
“colloidal suspensions” and as “colloidal solutions.” Actually there is a 
“continuous transition” from suspensions through colloids and on to so-called 
true solutions. (6) 


Particle Size 

Colloidal particles are often arbitrarily considered to be those falling in 
the approximate diameter range 1.0 to 100 millimicrons (10-7 to 10-5cm) al- 
though some authorities suggest a range of “over 1.0 millimicron to some- 
thing less than 500 millimicrons.”(6) Such distinctions are purely arbitrary 
and are based on expediency. 

For comparative purposes the following table of approximate particle 
sizes is presented. 


TABLE II 


PARTICLE DIAMETERS (7) 


Material Diameter 
(Order of Magnitude ) 
1 cm upward. 
Gravel qd 
Coarse Sand 0.1 em (10° em. ) 
Fine Sand -Oy1,cm (10° “cm. ) 
Silt 10_j,cm- 
Bacteria 10° ‘cm. (1 micron 
Clay 107 cm. 
Virus 10 gem. (1 millimicron) 
10 “cm. 


Oxygen molecule 
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Surface Area 

Some of the distinctive properties of colloidal particles can be explained 
by their large specific surface area. Assume that water contains 200 parts 
per million of a material having a specific gravity of 2.5 and that this mate- 
rial is colloidally dispersed in the water as spherical particles having a di- 
ameter of 1.0 millimicrons (10-6cm), Each 1000 gallons of water would con- 
tain about 6 x 1020 particles and their surface area would be approximately 
45 acres. Obviously surface activity plays a prominent role in colloidal be- 
havior. 


Charge 

Colloidal particles normally bear a charge. This charge probably can be 
obtained in various ways, including ionization of a part of the colloidal par- 
ticle. The charge may also be obtained as a result of the colloid’s surface 
activity; i.e., by preferential adsorption of various ions from the water in- 
cluding hydrogen (H+) or hydroxyl (OH~) ions.(8) The stability of many col- 
loids can be explained by the fact that the particles carry like charges, ~ 
Similarly charged particles repel each other; thus, continued dispersion 
rather than precipitation is favored. The degree of stability of a particle is 
determined by a value termed the “zeta potential” which is related to par- 
ticle charge and the distance through which the charge is effective. (8) 

Particle charge can be demonstrated experimentally using colloidal ma- 
terial found in natural waters. For example Saville,(9) working with a 
colored water, showed that the color was caused by negatively charged col- 
loidal particles. Saville actually demonstrated electrophoresis, a character- 
istic property of colloidal particles having to do with their migration under 
the influence of an electrical potential. Colloidal particles encountered in 
natural waters such as clay and organic color are commonly negatively 
charged. On the other hand colloids formed in the initial stages of the coagu- 
lation by aluminum and iron salts generally carry a positive charge. 


Precipitation of Colloidal Particles 


Colloidal particles can be precipitated by ions of opposite charge. The 
relative precipitating power of ions increases tremendously with increase in 
the ion’s charge or valence, For example sulfate (SO4--) is more than 100 
times as effective as chloride (Cl~) in precipitating positively charged col- 
loidal particles. The aluminum ion (Al***) is roughly 500 times as effective 
as the sodium ion (Na*) for precipitating a negatively charged colloidal par- 
ticle. The greater the concentration of dissolved mineral solids (ionic 
strength) the more readily precipitation occurs. (An increase in ion concen-. 
tration lowers the zeta potential), Therefore it is generally easier to accom- 
plish coagulation in waters of high mineral content than in those containing 
low concentrations of mineral matter. 

Inasmuch as ions responsible for charge on colloidal particles are often 
hydroxyl (OH-) or hydrogen ions (H+) colloidal-particle precipitation should 
be influenced by the pH value of the water, Water treatment experience sup- 
ports this statement. The process of coagulation is often quite sensitive to 
pH changes. Recognition of the importance of the pH factor in coagulation 
may mean the difference between good and poor results. 


ASCE HANEY 1036-7 


Coagulation Process 


Mechanism of Coagulation 


Black(3) whose studies of coagulation phenomenon have contributed greatly 
to our knowledge of the subject has divided the process into three phases. 

a) The coagulant furnishes trivalent aluminum or ferric ions (high plus 
charge) capable of neutralizing the negative charge on colloidal particles of 
color and turbidity. 

b) If the coagulation reaction occurs below about pH 7.0 the result of the 
initial action of the coagulant is the formation of positively charged particles 
of colloidal dimensions. Black has termed these “micro-flocs.” Their plus 
charge at low pH is probably due to adsorption of hydrogen ions and surplus 
aluminum or ferric ions. These positively charged particles have adsorptive 
capacity and neutralizing power and can continue to remove negatively 
charged particles of turbidity and color. 

c) The third or flocculation phase involves further adsorption, charge 
neutralization and finally agglomeration of the micro-flocs into larger and 
larger particles which will settle readily. 

It is apparent that two types of colloids are involved in coagulation, name- 
ly, those present in the water undergoing treatment and those formed by the 
coagulating chemicals added to the water. In general the former are nega- 
tively charged whereas the latter commonly carry a plus charge. However 
reversal of the charge on particles formed by the initial coagulation reactions 
may occur at approximately pH 7.0 and above. This is probably due to ad- 
sorption of hydroxyl ions. 


Factors Affecting the Process 


pH Value 
The pH value of the water affects the solubility of precipitates formed by 


iron and aluminum, time required for floc formation, and charge on colloidal 
particles. Effect of pH value on the time required for floc formation has been 
studied by Theriault and Clark;(10) Bartow, Black and Sansbury;(11) Black, 
Rice and Bartow;(12) Baylis(13) and numerous other investigators. Based on 
- theoretical considerations and certain experimental evidence Theriault and 
Clark concluded that when other influences are not considered optimum con- 
ditions for alum-floc formation will be found at about pH 5.5. From a prac- 
tical standpoint other influences are important but it is often the case that 
alum coagulation is best accomplished at fairly low pH values (5.0 - 7.0). 
This is particularly true of soft, colored waters. ; 
The pH range of relative insolubility for aluminum is approximately 5.0 to 
7.0, above pH 4.0 for ferric iron and above pH 9.5 for ferrous iron. 


Dissolved Mineral Matter 
Dissolved mineral matter in the water is a source of ions which can play 


an important part in the coagulation process. Coagulation of the positive 

colloidal particles (micro-flocs) formed by iron or aluminum coagulants is 

influenced considerably by other ions. The importance of any particular ion 
is toa great extent related to its valence or charge. As previously men- 
tioned, the sulfate ion which bears a fairly-high negative charge strongly 
favors precipitation of positively charged micro-flocs— much more so than 
_the chloride ion which carries a lower charge. Various investigators have 
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demonstrated that sulfate definitely extends the zone of optimum floc forma- 
tion in the lower pH range. The magnitude of this effect is somewhat depend- 
ent on sulfate concentration. Sulfate is a constituent of two commonly used 
coagulants, alum and ferric sulfate, and is found in a wide range of concen- 
trations in natural waters. 

At pH values of approximately 7.0 and above the charge on the micro-flocs 
may be negative rather than positive. If such is the case then positively 
charged ions particularly calcium and magnesium (high plus charge) would 
play a part in floc formation. 

In the complex process of coagulation probably all ions have some influ- 
ence. It is therefore apparent that natural waters, which vary widely in 
chemical make-up, will respond in various ways to the addition of coagulating 
chemicals. 

Effects of synthetic detergents on coagulation have been investigated by 
Vaughn(25) and associates at the South District Filtration Plant, Chicago, and 
by Smith(26) and others at the Robert A. Taft Sanitary Engineering Center, 
Cincinnati. 


Temperature 

The principal effect of temperature is its influence on the rate of floc 
formation. Generally the lower the temperature the longer the time required 
to produce good flocculation. Higher coagulant dosages may be required for 
cold water than for warm. Powell(14) has presented curves illustrating in- 
creased chemical requirements at temperatures below 50° F. Camp’s 
studies(15) showed that decreasing the temperature shifted the optimum pH 
to slightly higher values. Thuma(16) has reported recently on the adverse 
effects of low water temperature on water softening. 


Time-Concentration 

Various authorities(8)(17) have called attention to the fact that the coagula- 
tion process is a time-concentration phenomenon. High concentrations of 
coagulant produce results in shorter time than low concentrations. 


Coagulant Aids 

Discovery that partially neutralized sodium silicate solutions would sub- 
stantially improve coagulation was a notable advance in water treatment.(13) 
The silica activation process yields a negatively charged colloidal particle 
which not only hastens the rate of coagulation but also results in the forma- 


tion of larger floc particles. Activated silica is also beneficial in connection | 


with water softening by the lime-soda ash method and has been successfully 
used as a substitute for the conventional coagulants, (18) 

Lime, soda ash and acids (sulfuric and phosphoric) are coagulant aids in 
the sense that they can be used when necessary for pH adjustment. 

Suspended particles can assist flocculation by acting as nuclei for floc 
formation, Coagulation in some instances may be improved by the addition 
of small amounts of finely divided, insoluble material. Clay, bentonite, and 
activated carbon have been used for this purpose.(17) Nuclei may.also be 


provided by returning sludge from previously treated water. y 


Application of chlorine ahead of the coagulation process has been reported 
to be advantageous.(14)(19) Chlorine may alter or destroy materials which 
inhibit floc formation and render other materials more susceptible to coagu- 
lation. 

Polyelectrolytes are a recent development of considerable interest. These 


—_—_ 


ASCE HANEY 1036-9 


are synthetic organic compounds and are said to be useful in flocculating a 
wide variety of suspensions. One application has been in the treatment of 
suspensions encountered in ore processing. Apparently they can be used 
alone or in combination with conventional coagulants. 


Exchange Capacity of Turbidity Particles 

Investigations by Langelier and associates(20)(21) have indicated that an 
important factor in the coagulation process is the cation exchange capacity 
of the turbidity particles. When exchange capacity is fairly low, flocculation 


| requires a binder material such as alum or iron floc. An amount of coagu- 


lant proportional to the exchange capacity must be added to provide colloid- 


| destabilizing cations (Al+++ or Fe+++) and an additional amount is required 


to furnish binder material. Formation of the latter by coagulant hydrolysis 
occurs only after the turbidity particles have been destabilized by the coagu- 
lant cations, 

Mixing 

The importance of appropriate mixing cannot be over-emphasized. With- 
out it efficient coagulation is virtually impossible. It is important to recog- 
nize and distinguish between the two types of mixing commonly employed in 
coagulation. These are “rapid” or “flash” mixing and “slow” mixing or floc- 
culation. Rapid mixing requires only a few seconds using a high-speed mix- 
ing device. Flocculation on the other hand requires a much longer period of 
mixing at fairly low velocity. 

The first and probably one of the most important phases of coagulation in- 
volves the action of aluminum or ferric ions provided by the coagulant. It is 
therefore important that rapid dispersion of the coagulant through the water 
be achieved. The function of a rapid mixing device is that of quick and 
thorough dispersion of the coagulant. 

The initial phase of coagulation produces finely divided particles and op- 
portunity for coalescence must be provided if a readily settleable floc is to 
be obtained. Coalescence requires contact between particles and there is 
greatly increased opportunity for such contact if the water is slowly mixed. 
The result of slow-mixing is floc “build-up” or “flocculation.” (The process 
is often referred to as “conditioning.”) 

Flocculating velocities employed vary but often range from about 0.9 to 
1.5 feet per second. (22) Too high a velocity will prevent the formation of 
large floc particles or break up those that do form, Flexibility in mixing 
speed is desirable in most instances. Time allowed for slow mixing varies 
but an allowance of 30 to 45 minutes is common practice. Camp(23) has re- 
lated the rate of flocculation to velocity gradient and detention time. 

Following the flocculation phase of the coagulation process turbulence 
must be avoided. The floc particles are fragile and can be broken up in pas- 
sage from the flocculation chamber to the settling basin. The transition from 
flocculation to sedimentation should be as direct as possible and be accom- 
plished at low velocity. If conduits are required velocity in them should in 


general be less than 1.0 foot per second. 


Coagulation Control 


In view of the complexity of the coagulation process and the well-known 
variability of natural waters, it is not surprising that coagulation poses a 


" 
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difficult control problem for the water-plant operator. For practical operat- 
ing control there is no substitute for trial determinations of coagulant dosage 
conducted with laboratory stirring equipment.(24) Optimum dosage, pH ef- 
fects, and use of coagulation aids can be readily studied in the laboratory and 
it is essential that this be done. The results of laboratory studies, to be use- 
ful, have to be interpreted in terms of actual plant performance. For this, 
there is no substitute for experience. In this connection one further point 
needs emphasis. Laboratory data indicating the need for adjustments are of 
little value to an operator if the means of making the required changes are 
not available to him. Obviously there is often need for the provision of flexi- 
bility in the way of chemicals and facilities for their application, points of 
chemical feed, and mixing. The coagulation process cannot be forced into a 
rigid pattern, Although coagulating chemicals have been used in large-scale 
water treatment for more than a half century, we still have much to learn 
about their behavior and the nature of the substances dispersed in natural 
waters. 
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A FLOW CONTROLLER FOR OPEN OR CLOSED CONDUITS 
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(Proc. Paper 1037) 


SYNOPSIS 


: 


A general principle of flow control, utilizing a non-linear resistance, is 
veloped that is applicable to both open and closed conduits. The control is 
initely adjustable for discharge over its design range, and holds the flow 
Ose to the predetermined discharge over its head range. The control is 
adily converted into a flow meter with a linear head-discharge curve and 
justable range and sensitivity. Experimental work is presented to confirm 
ne theory. 


INTRODUCTION 


- The combination of a disc moving within a profiled throat, and a nonlinear 
esistance to support the disc against the pressure drop is shown to yield a 
rinciple of flow control that is infinitely variable and theoretically exact (i.e., 
ith no change in discharge or drag coefficients). By reversing the throat 
ection and measuring the head across the disc, the device is converted into a 
ow meter. The underlying principles of single orifice flow control are out- 
ned briefly before developing the specific relationships and design data. A 
10del has been constructed and tested in the Civil Engineering hydraulic 
uboratory of the University of Michigan. The results are presented for both 
1¢ flow control and the flow meter. 


Single Orifice Flow Control Concepts 


An orifice flow control must act to hold the discharge, Q, constant in 
Q = Cy av2g h (1) 


y changing the orifice area, A, as the head, h, varies over its design range. 


ote: Discussion open until January 1, 1957. Paper 1037 is part of the copyrighted 
‘Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
82, No. HY 4, August, 1956. 


Prof. of Hydraulics, Civ. Eng. Dept., Univ. of Michigan, Ann Arbor, Mich. 
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Cq is the discharge coefficient, and g the acceleration due to gravity. Assu 
ing the discharge coefficient constant, the area must vary inversely as the 
square root of the head. For a control meant to operate at a single dischar 
only, this may be accomplished by means of a linear spring assembly, Fig. 
Head across the device displaces the disc, thereby altering the net area. 1 
throat may be profiled so that the discharge is held constant. 

If one attempts to make a variable flow control from this device, by vary 
the position, z, of the spring base BB, Fig. 1, the dimensionless head-dis- 
charge curves appear as in Fig. 2, with Qg the design discharge, ho the mir 
mum design head, and 7) the ratio of maximum design head to minimum des 
head. 

By selecting another approach, although still retaining a linear spring, t 
throat profile may be constructed so that the discharge is the same at mini 
mum and maximum heads, for any setting, z, Fig. 1, within its size range. 
The discharge at intermediate heads varies from that at the end points, hov 
ever, as indicated in Fig. 3. To obtain the area required, as a function of x 
(Fig. 1), let Y be the displacement of the disc for head change from ho to hy 
i.e., from minimum to maximum design head. The area must change from 
A, to Ag as given by 


AV ho = AgV hoy ( 


which must hold for any setting, z, of the spring base BB. As the area of 0: 
fice opening is a function of x only, A = f(x"), 


f(x')Vhg = f(x’ - Y) Von ( 


Since this relation must hold for any initial x', a law of variation of A with: 
of the following form is indicated 


A = Age ~* ( 


where Ao, @, are undetermined constants. Substituting into Eq. (3) 


‘ a(x'- Y) — 
Age aes Vv ho = Age hoy 
or, simplifying 


a= ing ( 
2Y 

where ln is the natural logarithm. 
The variation of A with x' then becomes 


x'1n7 
A = Age ( 


with Ap the opening between disc and throat when x = 0, ie., Ap = Amin- 
To improve the accuracy of this control for the same head and discharge 
range, two linear springs may be introduced into the assembly, as in Fig. 4 
altering the first spring so that the disc is displaced 5 by the first spring o 
as the head varies from ho to hoVN and is displaced the remaining distance 
vale 


as both springs act with the head changing from hoVn to hoy. This results 
a dimensionless head-discharge curve, Fig. 5, in which the control is accur 
for three heads and the maximum error has been reduced. : 
By taking n linear springs, the first acting over the distance Y/ n, the fir; 
three over 2Y/n to 3Y/n, etc. with the corresponding head ranges 4 
4 
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hy to hon!/2, non!/2 to hon?/, hon?/™ to hon 3/2, ete. 


| the discharge is maintained the same at n + 1 heads with the maximum varia- 
| tion in discharge decreasing as n increases. By increasing n indefinitely, a 

| theoretically exact flow controller is obtained. The spring system then be- 
comes a non-linear spring. The relationships for this case are developed in 
ba following section. 


) Basic Equations for Nonlinear Resistance Flow Control 
| The basic equations for the flow controller have been derived by the author 
}in a different manner in a previous paper.! One form of the nonlinear resist- 
ance flow control is shown in Fig. 6. The spring element is of the restrained- 
tip cantilever type, in which the two leaf springs are forced together by a pro- 
filed backing that shortens their effective length. The spring element is a 
distance y = Y apart at minimum design head and contacts at y = 0 for maxi- 
/'mum design head. By rotating the milled knob, the discharge setting, z, can 
'be varied. From the geometry of Fig. 6 it is evident that 


x' =y+z (7) 


‘where x' gives the disc position (the area of opening between disc and throat 
is a function of x' only); and y is a function of the pressure drop across the 
|dise only, i.e., it is a function of h only, for a given fluid. 

For any setting z 


| 

Avh = (A+ AA) vh+Ah (8) 
}or, since A = f(x") 

: Vh f(x') = Vvh +Ah f(x" + Ay) (9) 


because a change in Ah causes a change in Ay and Ax' = Ay for constant z. 
| Expressing the area of opening, A, as an exponential function of x', Eq. (4) 


Cixi) =! Mee (10) 
where « is as yet undetermined. Substituting Eq. (10) into Eq. (9), and solving 
for a, 

[, 4b 
es a 
Ay 
Expanding the logarithmic term into a series and taking the limit as Ay ap- 
proaches zero, 
— f[Lah _ {| ah@ | 
OL, Se SS Se scene aviaiel is) 4 

lee a ; y z2h?ay 

or 
= Cea dh (11) 
2h dy 


Since « is a function of x' only, and y is a function of h only, Eq. (11) shows a 
function of x' equal to a function of y'. This is impossible because of Eq. (7), 


1. “A Quick Response Variable Flow Control Device,” V. L. Streeter, Instru- 
ment Society of America, Vol. 2, No. 2, Feb. 1955, pp. 48-51. 
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hence both sides of Eq. (11) must be constant. Integrating Eq. (11) 
-24y =Inh+ete 


Using the end conditions: h = hg, y = Y; h= hon, y = 9, both w and c are de- 


termined, 
bo 
haere (12) 


Yo buh 
LAREN Se (13) 


and 


Equation (10) becomes 
x : In ” 
Av=A e (14) 


Now, substituting into Eq. (1), using h from Eq. (13), A from Eq. (14) and x' 
from Eq. (7) 


z in 
Q = Cg Amin¥ 2€ bon e—2Y— (15; 


which yields the discharge for the given arbitrary setting z, and shows that 
is independent of head h. 

The three equations, Eqs. (13), (14), and (15) determine the characteristic 
of the flow control. The nonlinear resistance must follow the head-deflectior 
relation of Eq. (13), and the throat profile must be made so that the net area 
of opening is given by Eq. (14). By providing for the base of the resistance t 
be adjusted axially, the quantity z may be changed and variable flow control 
results, as given by Eq. (15). 


Equation for Throat Profile 


If an axially symmetric throat section is desired, the throat profile, as 
given by the cylindrical coordinates (r, x) in Fig. 7, must yield the net area 
for any x' required by Eq. (14). Using Pappus theorem this area is 


An 20 Foe V(r = te) ogee sen)? 
x'1n¥ 
= Amin ©  2¥ (16 


or 


x'1n7 


12 2 / 
(r= Fo.) %e°* 2Y eur or PLN (roe ro)” Pe wae t= (17 


For a given x' this equation yields a curve of r as a function of x. Drawing 
segments of these curves for various values of x', Fig. 8, shows that their 
envelope yields the desired throat profile. The partial differential of Eq. (1’ 
with respect to x' is taken and solved simultaneously with Eq. (17) to yield tl 
cylindrical coordinates of the throat profile. Dividing the partial of Eq. (17) 
with respect to x' into Eq. (17) yields, 


gy ore | — [ort : 
ote ral | Y ) (18 
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Substituting Eq. (18) into Eq. (17) gives the desired equation, 


B - rin? _ |)_[CretedIn 
ae nf V2 J oteglnt pany 4 egy (19) 


| This equation permits ne direct computation of x for any r in terms of the 


design quantities Y, n, r', and ro. 
There is an inherent eawees to the value of r in Eq. (19), obtained by 
Setting the radical term equal to zero, 


Imax =To + .== (20) 


| 
For larger values of r the envelope solution breaks down. Forr=r from 
“Eg. (19) ses 


| 
Xnax = = “lh (2B. va 5%) -+| (21) 


The corresponding value of x', from Eqs. (18) and (21), is 


Kar = BY (2st BEV fG) . 


ly See, re 


Discharge Limitations 


For the flow control to operate the disc must stay within the range 
0<x'<X},x- For control over the whole head range this restricts z to the 
range 0 < Z < Xmax - Y. To find the maximum discharge over the whole head 
range, Z = X‘max - Y is substituted into Eq. (15) using Eq. (22), 


Y 
Gaexi Fuca 27 Veiho nr (2r9 + p>) (23) 


making use of the fact that Amin = m(r'2 - T°). The minimum discharge is 
obtained from Eq. (15) by setting z = 0, 
Qmin = &q Amin ¥28 hoy (24) 


Equations (23) and (24) permit the selection of the three design quantities ro, 
r', and Y for given head and discharge range. It is of interest to note that 
the ere and minimum discharges are not related. The minimum throat 
radius r' does not come into Eq. (23). 

For flow control over partial head ranges, the limitation is that Z varies 
from -Y to 0 in the lower extended zone and from X max - Y to x'max in the 
upper extended zone, with the disc (x') remaining between o and x\max at all 
times. The spring travel is limited to portions of the range Y to 0, which re- 


stricts the head range to part of ho to ho?. 
For z less than zero, the “threshold” value of y, to maintain x' = 0, is 


y=-z 


: 
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Substituting this value of z into Eq. (15), using Eq. (13) for y, 

Q = Cg Amin V 26h (25) 


This is the lower limit line of the lower extended zone. 
For z greater than x'max - Y, and with the disc at its limiting position 


x" By X'max, 
Z=Ximax- Y 
Again, using Eqs. (15) and (13), the upper extended zone is limited by 
X max 197 
Q = Cq AminV2gh e” 2Y (26) 


The total flow control range is shown in Fig. 9 in dimensionless form, with 
the z range indicated on the right side. 


Design of Nonlinear Resistance 


The nonlinear resistance to displacement of the disc, as required by Eq. 
(13), may be expressed in terms of force exerted on the disc. This force, F, 


1S 
F = Cg@ 1,2 wh (27) 


where Cq is the dimensionless drag coefficient, and w is the unit weight of 
liquid. The minimum design force, Fo, is then 


Fo = Ca7ro* who (28) 


Cq is assumed to be constant over the range of travel of the disc. Substituting 
for h and hg in Eq. (13) 


Y Fon 
oe Sot 
y int In = (29) 


The resistance may be obtained in many ways, three of which are dis- 
cussed here, namely 1) fluid resistance, 2) dead weight resistance, and 
3) spring resistance. 


Fluid Resistance. Referring to Fig. 10, a bellows with area, Ap forces 
manometer liquid into a manometer leg of varying cross-section. The area, 
a, of manometer is determined as a function of s. The necessary equations 
are Eq. (29), 
-Ap dy = aydsy =a ds (30) 
and 
F = Ap (Wo - w) (s + sy) (31) 


where aj is the constant cross-sectional area and Wo is the unit weight of 
manometer fluid. F and y may be eliminated from the last three equations to 
yield a and s in terms of sq 


a (32) 


ye Sai ae Sat 
tln” S1At'p 71 = 
ALY Soe bY | 
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and sa Inn 
S=s5,6 MY = s, (33) 


where So is the value of s for Sj =Oand F=F >. The discharge setting is 
altered by changing z. For vertical installations the weight of moving parts 
would change the setting z for a given discharge, but would not otherwise dis- 
turb proper functioning of the controller. 


Dead Weight Resistance. A weight suspended on a cable or metal tape such 
that it is lifted out by a cam can easily be arranged to yield the nonlinear 
force-displacement relation required by Eq. (29). Figure 11 illustrates one 
method, with the cam rotating about a transverse axis. To determine the 
shape of the cam, reference is made to Fig. 12, where the cam rotates through 
the angle, a, from o when F = Fo, to Q@, when F= F,7. The moment of the 
weight, W, about the shaft for any angle & is 
-y int 
We CoS(@-X) = Feo = Fo poy if (34) 


With po constant, 


2 
alte 


since pyp&= Y- y. Letting W = Fo, Eq. (34) reduces to 


alnn 
pcos (8- &) Be BA, Cron ee (35) 


For any a, Eq. (35) is a straight line in polar coordinates (p,@). The envelope 
of the family of straight lines is the cam profile equation. Taking the partial 
differential of Eq. (35) with respect to & and solving simultaneously with 

Eq. (35) to eliminate a@ yields 


p= vir(eye ™ pan (36) 


which is a logarithmic spiral. By dividing the partial differential of Eq. (35) 
with respect to a by Eq. (35) 


tan (9- a) = 13? (37) 
showing that @ -a@ is aconstant. @ varies from tan~! (Inn/o9) to &p + tan=! 


(Inn/o,) in Eq. (36). 
If eran is to be designed to support the submerged weight, G, of the mov- 


ing parts of the controller, 


a inn 7 
Fo po © A ope Gee Wp.cos (6 X) ian) 
Taking the partial differential with respect to a and dividing by Eq. (38) 


agtan(@- %) = a ay (39) 
cree 
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For a series of values of & between 0 anda, corresponding values of @ are 
found. Using these values the corresponding values of p are found in Eq. (38), 
yielding polar coordinates of the cam profile. 


Spring Resistance. A restrained-tip cantilever spring may be designed to 
provide the required resistance law by following a method due to Mr. S. P 
Clurman.? Two identical, flat springs of constant width and thickness are 
placed as shown in Fig. 13. As the force, P, is increased the springs are 
forced against the profiled backings; thereby reducing their effective lengths 
and increasing their stiffness. The springs are fastened together at each end 
by a spacer block that prevents rotation. As there are actually four springs, 
each one takes half the load and is deflected through half the total displace- 
ment. The equation of the backing profile is determined as follows: 

Letting G = KF, be the submerged weight of moving parts of the controller 
for vertical stem applications, (G = 0 for horizontal applications) 


-y ln 
P= F-G =F, (yey - ®) 


The sign of G has been taken for the case where the weight of moving parts 
acts in the opposite direction to the pressure force on the disc. With 6 the 
deflection of one spring from the position for minimum design head 


1 Y ; 
9 ie ag ace a Ya a In (ee + k) (40) 


where P' is the load on a single spring (P/2). For the minimum design load 
the spring is cantilevered from its full length, ), with the backing touching at 
Xo = 0, Fig. 14. The variables xo, yg are coordinates of the backing profile. 
With no load on the spring it will be displaced from the 5 = 0 position a 
distance 

Fo(i-k) £9 

oa Le 


where E is the modulus of elasticity and I the section moment of inertia about 
the neutral axis. 

With (Xo, Yo) the last point of contact between spring and backing for a 
given active load P', the deflection from no-load position is 


3 2 42 
Fo(1-k)2 (2-X0)*% d“yo 2 dYo 
Serle EI = <a “dx52 + 3 (L-X9) ax + Yo. (41) 


The right-hand side of the equation is obtained by consideration of deflection 
of cantilever beams, as follows: In Fig. 15, by superposition of the two simple 
loadings 


6 + 


dyo _ P'(L-x9)” — M(2-x0) 
dks 2 Galt (42) 


Defl. = P(L-x0)3 _ M(L-x0)? 
3 EI 2 EI 


Replacing M by EI dyo/ dx52 and eliminating P' in the two equations 


2. “The Design of Nonlinear Leaf Springs”, by S. P. Clurman, Transactions, 
American Society of Mechanical Engineering, Feb. 1951, pp. 155-161. 


) 
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2 
1 2day 2 dy 
Defl. = (i-x, ) Oo + KE=XA 7S 
. ax? 4 ;. dxo 


which, with the addition of yo, is the right-hand side of Eq. (41). 
6 and F, may be expressed in terms of the stiffness of the spring. From 


Eq. (40) 
261n 
p'=Fo (e YY ~- x) 
2 
Then 
. 2d 1n7 
OP tty 
T a ae: 4 (43) 
which is the spring stiffness. For 6 =0 
dP’ Foln7 =oK = 12 EI (44) 
oF gH 0 nd L 
The stiffness of a spring varies inversely as the cube of its length, hence 
Ng aad Sak See 
dé ~ (£-x9) 
Solving Eq. (43) for 6, using the last two equations, 
eat Dh RS eee we 
ne 1n7 a8 estes | 2 1n7 L-Xo (45) 
From Eq. (44) 
Fo (1-k) 
24 EI tay 7H Bk 7 


Substituting the last two equations in Eq. (41) 


3Y Y 1 d2 2 
re oe ean G'S) = exe + = (h-x,) 72 + Yo (46) 
2 1ny L-Xo 2 1n7 6 dx5 3 dx, 
Inserting the dimensionless quantities 
sc a yr yo 
ey 
2- = 
OCS Kt ies) = ee) oe 4 (lex) SZ + 67 
n(1-x) (1-x) a + 4( ax 


In} 1n7 
This equation is converted into a linear equation with constant coefficients by 
the substitution of 


t = In(l-x) 
Then 
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and Eq. (47) becomes 


2— a 
ary. 25 SY are? t1-k-3k) 48 
az 5 1 oe 6y ins ( (48) 


The solution of this equation for the conditions 


y= 0 20S 89556 
y ok ‘ae 


is 


3 
oy Inq = -k- 3 - 32% 3G + ®) e2> _ 2xe3$ 


Substituting back for x, 
2y ing = -k - 3 - 3 1n(1-x) + 3(d+k) (1-x)?- 2k(1-x)? (49) 
This is the dimensionless equation for spring backing profile. 


From Eq. (45) the maximum value of X = x,/8 is, for § = Y/2, 


26 thas -1/3 
Rena we gin (50) 


The cross-section of the spring must satisfy Eq. (44). 
Consideration of stresses in the spring indicates the maximum moment at 
the spring extremity. From Fig. 16 the moment Mj is given by 


M) = P L(1-X) - My (51) 
Mg may be found from Eq. (42), by obtaining dy,/dx, from Eq. (49) 


dyo i '¢ dy . x 3 = = =\2 
aXo = Z an = SP inn Inq Ee 3(1+2k) (1-x) + xa 
and 
M cag 8 Eee Me By Peed bag @ = ane 
pe L-8)_ = Zpyl-8)*- Zin E - $142k) 1-8) +6K0-9) (52) 
From Eqs. (40) and (45) 
" = Fo 
Poe ae ) 
and from Eq. (44) 
Ege. Fob 
ing 12 EI 


Substituting the last two equations into Eq. (52) and simplifying 


6 Pee Faia . 
Mo= = ae oh ae 2K0-%) 


Replacing Mg in Eq. (51) by this value, 


Fax 3 
My = “9° Reese dt Ws 2K0-9)| 


M, is a maximum for x = 1 -7 -1/3 


: 
: 
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Fol 
Se (a7 Seloaky +43) (53) 


The maximum fibre stress, S,,, is 


ee Niger, 0 Fo & 2/3 -1/3 
See ee 87 - 1 - 2k7 ) (54) 


\where t is the spring thickness and b the spring width. Eqs. (44) and (54) per- 


mit the thickness and width to be expressed in terms of Sm» E and the design 
quantities 7, Y, k, and Fo, 


t, = 4 Sn 5 1n7 
Bos cSikY ee d8 1 = 2k7 1/3) (55) 


y] _- 
ta 
2 
205 275 F- CEY) 2/3 2 -1/3,3 
| be = 2 Gd ies Cinq)? (37 1 2k 7 ) (56) 


‘The restrained-tip cantilever spring is relatively simple to construct and 
|has very small hysteresis loss. Its main objection, from the standpoint of a 
iclosed conduit controller, is its size. 


Flow Meter with Adjustable Sensitivity and Range 


__ If the throat casting in the flow control unit is reversed, as in Fig. 17, and 
b manometer is connected to the two ends of the throat, the pressure drop 
across the disc, as given by R the gage difference, is a linear function of the 
discharge. The relation among the variables x', y, z' is 


: 


x 4-y-= Z* (57) 


where z' + z is a constant. Substituting into Eq. (1), using Eqs. (13), (14) and 
(57) to eliminate x', y, and A, 


z'1n7 
| Q = Cq Amin / 2e e h (58) 
ij ho? 


‘showing Q to be a linear function of h. 
The sensitivity and range is changed by adjusting z, or z'. Equation (58) 
may be arranged into the dimensionless form 


z‘lng7 
ee: ee Ts (59) 
Cq Amin¥ 2Zghon 7 No 


which is plotted in Fig. 18. 
The upper threshold of the metering range is given by 


X max Iny 
Q epee ong (60) 
t i Cq Amin¥ 22h} it No 


The lower threshold of the metering range is 
i 


Rie, 
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Q = | h (61) 
CqAmin V26 hoy hoy | 


Large settings, z, result in large gage differences for a relatively small dis- 
charge range, while small settings permit metering over a much wider flow | 
range. The sensitivity of the meter, as given by Q/h, is an exponential func- 
tion of the setting Zz. 


Experimental Program 


An experimental flow control unit was designed to operate in a concrete 
flume in the floor of the Civil Engineering Hydraulic Laboratory. The unit 
was placed into a circular opening in the hérizontal shelf of a barrier, Figs. 
19 and 20. A 45° v-notch weir was placed a short distance downstream from 
the controller. Hook gages were mounted upstream and downstream from the 
barrier, so that the head across the controller, as well as the discharge, 
could be determined. The disc has a displacement of 0.75 inches as the head 
varies from a minimum of 2 inches to a maximum of 8 inches. The disc di- 
ameter is 7.50 inches and the minimum throat diameter was designed to be 
7.60 inches. The minimum throat section was later measured to be 7.605 
inches. 

With the design quantities Y = 0.75 inches, h, = 2 inches, n = 4, rg = 3.75 
inches, and r' = 3.80 inches, cylindrical coordinates of the throat profile wer 
computed from Eq. (19). The restrained tip cantilever spring, Fig. 19, was 
used for the nonlinear resistance, and was computed from Eqs. (49), (55), anc 
(56) for = 3 inches and k = 0, i.e. not allowing for weight of the disc assem- 
bly. The submerged weight of disc assembly was balanced by use of styro- 
foam floats placed on a downstream extension of the shaft. The amount of 
buoyancy was adjusted so that y = 0.75 inches when h = 2 inches. Originally 
the drag coefficient for the disc was assumed to be unity. Based on this co- 
efficient the spring bottomed at h = 6.5 inches. By adding extra width to the 
spring the resistance was increased so that the disc displacement was 0.75 
inches for head change from 2 inches to 8 inches. 

The styrofoam floats also served another purpose by stabilizing the disc 
for the larger discharges. Without stabilization the disc tended to oscillate 
through its whole stroke for discharges greater than about 50 g.p.m. A steel 
can was anchored to the floor of the flume concentric with the disc and floats 
with about 0.25 inches radial clearance between can and floats. The resist- 
ance to forcing water into and out of the can provided adequate stabilization 
for all the experiments. Limitations in flume height restricted the range of 
heads for the larger discharges. 

Fins were placed in the throat section to guide the disc. These were foun 
to cause excessive friction, and were cut back radially to avoid contact with 
the disc. A teflon bearing in a spider at the small end of the throat provided 
excellent guidance with a minimum of resistance. 

Test results for the flow controller are shown in dimensionless form in ~ 
Fig. 21, with a summary of results in Table I. 

The throat casting was inverted to convert the unit into a flow meter. He: 
across the disc was measured with hook gages and discharge by means of th 
weir. Test results are shown in dimensionless form in Fig. 22, with a sum- 
mary of results in Table II. 
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TABLE I 
SUMMARY OF RESULTS---FLOW CONTROLLER 


Number of 


turns on Maximum 
adjusting Head Range Discharge Divergence 
screw inches g.p.m. percent 
0 2.32-8.27 20.0 1.4 
3 2.17-8.03 24.5 1.6 
6 2.22-8.12 30.0 2.0 
9 2.12-8.72 35.8 2.2 
12 2.34-8.06 44.8 1.3 
15 2.20-8.88 54.5 2.8 
18 2.29-8.46 66.6 1.9 
21 2.14-7.91 80.8 2.4 
24 1.92-7.59 99.4 2.8 
27 1.88-7.17 123.8 3.9 
TABLE II 


SUMMARY OF RESULTS---FLOW METER 


Number of 
anaes rae 
screw Rinches percent 
(¢) 4.229 2.8 
3 5.225 4.2 
6 6.559 3.3 
9 8.096 3.2 
12 9.876 3.2 
15 11.972 2.4 
18 14.773 3.1 
21 18.334 3.0 
24 22.728 1.5 
27 28.372 3.7 
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Discussion of Experimental Data 


As only one model was constructed, there has been no opportunity to take 
advantage of experimental findings to increase accuracy. If a change in dis- 
charge coefficient occurs due to change of disc position within the throat, it 
may be compensated by redesigning the throat so that CgA varies with x'as 
prescribed by the theory. In the present model no instrumentation was pro- 
vided to determine the exact position of the disc for a particular setting. The 
setting z could be altered either by rotating the disc on its shaft (4 threads to 
the inch) or by rotating the adjusting knob over the spring (13 threads to the 
inch). The turns referred to in Figs. 21 and 22 are for the adjusting knob 
above the spring. 


Flow Controller. For the flow controller the “0 turn” disc position at h = 8 
inches was about at the minimum throat section. The average value of dis- 
charge for this setting, 20.05 g.p.m. was taken aS Qmin- As the exact positior 
of the disc in the throat is not known for this setting the area is unknown and 
likewise the discharge coefficient. 

Separate measurements of discharge coefficient were made with the disc 
in known, fixed positions, and without the upper supporting framework for the 
spring assembly. These results indicated a coefficient of about 0.62 within 
the throat and about 0.63 with the disc slightly above the throat casting. 

In each series of observations for flow controller or flow meter measure- 
ments were taken about equally for increasing heads and for decreasing heads 
Hysteresis effects appear to be very small, much less than the s-shaped tren 
of the data, which apparently is caused by the spring. The width of spring wa: 
varied in the model, as well as the amount of styrofoam, until the displace- 
ment Y was 0.75 inches at h = 2 inches, and was zero at h = 8 inches. The s- 
shaped trend in the data may be explained by the spring backing being incor- 
rect to the extent that the spring was too stiff for the head range 2 to 5 inches 
thereby holding the disc higher in the throat and increasing the flow area 
slightly. Also the spring was not quite as stiff as required by the theory for 
the range 5 to 7.5 inches and permitted the disc to be lower than required by 
theory, thereby decreasing the discharge. 

Figure 21 is plotted with log Q/ Qmin 2S ordinate so that the same relative 
percentage error scale occurs for all settings, and so that the average dis- 
charge lines would be equally spaced. A change of three turns represents a 
displacement Az of 3/13 inch. 

Log Q/ Qmin iS plotted against number of turns in Fig. 23, where Q is the 
average of the discharges for a given setting. The solid line is determined by 
least squares and yields the equation 


Q = 0.994 e0-06704t (62) 
Qmin 


Comparing the exponent with z Inn/2Y, z = 13.78t. As z = 13t by direct meas. 
urement, the discrepancy can best be explained as due to a change in dischar; 
coefficient. Tests with a fixed disc yielded a change in discharge coefficient 
when the spring support stands were in place that conforms generally to the - 
relationship z = 13.78t. ; 

The disc was stabilized by placing the steel can around the styrofoam 
floats for the settings 15 to 27 turns. 


STREETER 1037-15 


The “0 turn” position for the flow meter was near the minimum 
roat section for h = 2 inches. The discharge Q,,jn was taken as 20.05 g.p.m. 
5 in the case of the flow controller. The data as plotted in Fig. 22 generally 
mfirms the relationship that discharge varies linearly with head for any 
btting. A slight s-shaped trend in the data, similar to that of the flow con- 
oller, but reversed, can be explained by the spring being too stiff for small 
s and not stiff enough for the larger heads. 
| The slopes m = Q/h (Q in g.p.m., h in inches) for the best radial lines 
rough each set of observations have been computed by the method of least 
uares and are shown in Table Il, with the maximum percentage divergence 
‘any Single observation. The taking of data for the flow meter was rather 
ficult due to the large volume of storage upstream from the meter. Water 
S admitted to the flume from a constant head tank and the level in the flume 
me to equilibrium very slowly. Measurements were usually taken before 
eady flow was established, which could introduce minor errors in the ob- 
ervations. The flow meter is inherently stable and did not require any form 
; Stabilization. 
_ It was necessary to shift the disc on its shaft between the runs 0 - 12 and 
}- 27. Unfortunately, it was not adjusted so that the change from 12 turns 
15 turns was actually 3/13 inch. This is evident in the plot of log m/ Mo 
ainst number of turns in Fig. 24, where mg is the value of Q/h for 0 turns 
o = 4.229). By constructing the best line through the points for 0 - 12 turns, 
' for 15 - 27 turns, the relations between t and z were found to be z = 12.97t 
nd z = 12.87t respectively, which confirm the known relation z = 13t. This 
\dicates a substantially constant coefficient of discharge over the whole 
ange of tests on the flow meter. 


CONCLUSIONS 


Theory and design data for both the flow controller and the flow meter have 
2en developed, with several forms of nonlinear resistance. An open channel 
odel was constructed using a restrained-tip cantilever spring for nonlinear 
2sistance, and data was taken for ten settings each as a flow controller and 
3a flow meter. The results generally confirm the theory and point the way 
the construction of accurate and relatively simple control and metering 
juipment. 
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NOTATION 


A1,A, = area of opening between disc and throat 
a,aj = manometer cross-sectional areas 
ap = area of bellows 
b = width of leaf spring 
Cq = drag coefficient 
Cy = discharge coefficient 
E = modulus of elasticity 
e = base of natural logarithms 
F = pressure force on disc 
F, = minimum design force on disc 
f = functional relation 
G = submerged weight of disc assembly 
g = acceleration due to gravity 
h = head drop across disc 
hy = minimum design head loss across disc 
I = moment of inertia of spring section about neutral axis 
K, = Stiffness of spring of length Q 
k = ratio G/Fy 
Q = length of spring for load F,/2 
Q' = length of disc support shaft 
In = natural logarithm 


M = moment 


ae} 
ul 


force on spring assembly 
P' = force on single spring element 
Q = discharge 
R = manometer gage difference 
r = radial coordinate of throat 
r' = minimum throat radius 
r, = disc radius 
Sn = maximum fibre stress 
S,S,,8; = dimensions of manometer 
t = thickness of spring 
W = weight 
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w = unit weight of fluid 
X = axial coordinate of throat profile 
x' = position of disc in throat 
Xp = coordinate of spring backing 
x = dimensionless length x/Q 
Y = value of y forh=h, 


y = position of spring 


y = dimensionless quantity yo/Y 
Yo = coordinate of spring backing 

z = discharge setting 

a@ = constant, angle 

B = angular deflection 

6 = displacement of single spring element from minimum design head 
| position 
| € = dimensionless variable 

7 = ratio of maximum design head to minimum design head 
@ = polar coordinate 


P, Po = polar coordinate 
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Fig.!| Linear Spring 
Flow Controller 
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Controller 
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Fig-I| Profiled Cam Fig.l2 Notation for Cam 
Equation 


Fig.13 Restrained Tip Cantilever Spring 
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Me 
Position for Minimum Design Head 


Fig.14 Notation for Cantilever Spring 


~ dy P'(L-xX0)* 
ps R= Faxed 


Fig. 15 Deflections of Cantilever Spring 
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Fig.16 Free body Diagram of Spring 
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Fig. 19 Elevation Section of Controller 


ASCE 


STREETER 


1037-25 


KRiaw Cantrollear in Diara 


20 


Fis. 


1037-26 HY 4 August, 1956 


8 é 8 oa, 27turn 
3 BES °° 
7 00 © 0©0| 9° Oo 24 turns 


Fig.2! Flow Controller Test Data 
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SEVEN EXPLORATORY STUDIES IN HYDRAULICS 


Hunter Rouse, M. ASCE 
(Proc. Paper 1038) 


FOREWORD 


Immediately following the Sixth Hydraulics Conference, held at Iowa City 
in June of 1955, the Iowa Institute of Hydraulic Research conducted an inten- 
sive three-week course in experimental techniques, in which twenty-one 
graduate engineers from various organizations participated. The purpose of 
‘the course was not only to acquaint the class with new types of instrumenta- 
tion, but also to instill in them an appreciation for the methods of research. 
During the first part of the course, eight of the elementary experiments in 
mechanics of fluids taught at Iowa to engineering students were carried out 
and reported upon by the class in teams of three. These experiments includ- 
ed curvilinear flow with cavitation in a closed conduit; curvilinear flow with a 
free surface; viscous flow in a pipe; turbulent flow in smooth and rough pipes; 
form effects at conduit transitions; diffusion of a submerged jet; fall velocity 
of spheres; and flow past a circular cylinder. Messrs. D. W. Appel, P. G. 
Hubbard, E. M. Laursen, and D, E. Metzler of the Institute staff acted as 
instructors during this introductory period. 

The second part of the course involved eight experiments utilizing ad- 
vanced methods of measurement: the Pitot cylinder for constructing a two- 
dimensional flow pattern; resistance wires, pressure cell, and Sanborn re- 
corder for wave studies; the electrolytic tank for pressure distribution at a 
three-dimensional inlet; an electrical analog computer for the cross-thrust 
on neighboring cylinders; the hot-wire anemometer for submerged-jet turbu- 
lence; stagnation-tube surveys for boundary-layer shear on a plate; multiple- 
resistance gage for scour at a bridge pier; and a strain-gage dynamometer 
for lift and drag on a foil. The foregoing instructors were assisted in these 
experiments by Messrs. M. Hug, S. C. Ling, T. T. Siao, and A. Toch. 

In the third part of the course, which is of primary concern herein, the 
members of the class were encouraged to undertake short investigations of 
an original nature, to the end of determining new and publishable information 
through application of the foregoing techniques, The following seven topics 


Note: Discussion open until January 1, 1957. Paper 1038 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
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1. Director, Iowa Inst. of Hydr. Research, State Univ. of Iowa, Iowa City, 
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were suggested: 


Development of the non-circulatory waterspout 

Penetration of a jet into a counterflow 

Surface profiles at a submerged overfall 

Coefficient of contraction for a submerged jet 

Maximum amplitude of oscillation of a cylinder due to its vortex wake 
Deflection of a jet by a normal wall 

Influence of slope and roughness on the free overfall 


Experimental facilities for the investigations were improvised from equip- 
ment at hand, and a member of the Institute staff was assigned for consulta- 
tion to each of the seven teams. In the short time available, each team per- 
formed the required experimental measurements and computations, and six 
of them submitted preliminary reports to serve as a record for reference 
purposes. The seventh team, under the leadership of Mr. J. C. I. Dooge, was 
able to give somewhat more time to its project, and the report that he finally 
submitted was well beyond the preliminary stage. Such amplification of the 
other reports as proved necessary to prepare them for publication in the fol- 
lowing form was carried out subsequently by the advisers under the writer’s 
editorship, and the illustrations were brought to their present state of uni- 
formity by Mr. Metzler. Although none of the articles even approaches the 
final word upon the subject, it is believed that each will provide both useful 
information and an incentive to further research. 
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DEVELOPMENT OF THE NON-CIRCULATORY WATERSPOUT 


J. Davidian,? J. M. ASCE, and J. E. Glover3 
(D. W. Appel,4 J. M. ASCE, Adviser) 


Statement of the Problem 


The most familiar occurrences of waterspouts in nature are those accom- 
ying cyclonic storms over the oceans. A vortex is formed by the inter- 
tion of air masses, into the low-pressure core of which water is drawn, 
metimes to great heights. The circulatory motion of the vortex is evident- 
r an essential part of the phenomenon, Less well known is the non- 
irculatory waterspout, which may form in the low-pressure zéne produced 
7 radial flow toward an intake in either of two fluid layers of different den- 
ty. In plumbing systems the phenomenon is known as back-siphonage, 
ter as well as air being drawn into the unsubmerged end of an open supply 
pe that is momentarily under vacuum. The purpose of this study was to 
termine experimentally the rate of flow necessary to develop a waterspout 
st to the point at which the second as well as the first fluid is carried into 
xe pipe. This represents the critical limit for safeguarding water-supply 
mes from contamination, and for the general separation of liquids of nearly 
f same density by decantation. 
_ The geometric variables in the problem are the diameter D of the suction 
pe and the vertical distance y from its open end to the interface between 
ie fluids. If one assumes for the purpose of a preliminary analysis that the 
fects of viscosity and surface tension are negligible, the only significant 
uid properties are the density of the upper (or moving) fluid P and the dif- 
rence between the specific weights of the two fluids Ay. As the flow vari- 
ale, the mean inlet velocity V at which a waterspout forms has been cho- 
2n. When these variables are combined by using the techniques of 
mensional analysis, the two dimensionless parameters obtained are y/D 
id pV2/Ay D. The first of these is a geometric variable, whereas the latter 
presents the ratio of inertial and gravitational forces and hence is a form 
the Froude number. This analysis indicates that a single series of mea- 
irements of the critical velocity for the formation of a waterspout over the 
actical range of relative heights of the suction pipe should suffice to deter- 
ine the function 


z= ¢ (25) 


wever, to check the validity of the assumptions made in the analysis and to 


Hydraulic Engineer, U. S. Geological Survey, Atlanta, Georgia. 
Hydraulic Engineer, Waterways Experiment Station, Vicksburg, Missis- 


Sippi. 
Research Engineer, Iowa Institute of Hydraulic Research, State University 


of Iowa, Iowa City. 
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demonstrate the generality of the results, two series of tests were conducte 
In the first series a waterspout was formed in water under air, and in the 
second the waterspout was formed in saline water under fresh. 


The Experiments 


To simplify observation and extend the range of the first series of tests, 
three suction pipes 0.5, 1.5, and 3.0 inches in nominal inside diameter were 
used interchangeably. Each of these in turn was mounted vertically above a 
open tank of water 3 feet in diameter and 3 feet deep in which the elevation ¢ 
the water surface could be regulated by means of supply and drain lines and 
measured by a hook gage (see Fig. 1). The inflow of air necessary to form 
the waterspout was produced by connecting the pipe to a large tank from 
which the air had been exhausted to a pressure of 10 inches of mercury abst 
lute. 

As shown in the figure, a differential manometer was connected between 
stagnation tube and a wall piezometer in the horizontal portion of the suctio: 
pipe for measurement of the centerline velocity at the gaging section. This 
was related to the mean velocity V = Q/ 1m7D2 at the inlet section by means 
of the Bernoulli equation, through the assumption of adiabatic expansion be- 
tween the inlet and the gaging section. The pressure at the inlet was read o 
a barometer and the pressure at the gaging section was obtained from a 
Bourdon gage. The coefficient of the flow meter was determined by calibra 
tion with water. 

With the water surface set at the desired elevation for a test, the control 
valve was slowly and continuously opened until the waterspout that develope 
just began to carry water into the pipe. At this instant the two gages were 
read at the flow-measuring section. The value in the suction line was then 
closed, the vacuum tank pumped down, and a different water level establish 
in preparation for another test. 

In the second series of experiments a glass-walled tank 1 square foot in 
cross section and 2 feet high was used, This permitted visual observation 
of the interface between clear fresh water above and dyed salt water (speci 
gravity = 1.1) below. The arrangement of the suction pipe and vacuum tank 
was the same as before. The rate of flow, however, was determined by tim 
the fall of the free surface in the tank, which was quite feasible with the lov 
rates required, Because of the limited size of the tank, only the 0.5-inch a 
1.5-inch suction pipes were used. 


The Results 


All experimental data are presented graphically in Fig. 2 in the dimensi 
less form suggested by the preliminary analysis. It can be seen that the te 
using different sizes of pipe as well as entirely different fluids produced re 
markably consistent results. Some scatter was to be expected, particularly 
for the 1/2-inch pipe in the two-liquid test, because of the difficulty in regu 
lating the rate of flow to obtain critical conditions. The tendency of the poi 
to follow a single curve (approximated by the power relationship 


ab 
y/D = 0.42[ V/y DAY/p ]?) indicates that the assumption of negligible visco 
and capillary effects was justified for such an exploratory study. A 
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otograph of the waterspout in air is shown in Fig. 3, and one of the saline 
terspout in fresh water is presented for comparison in Fig, 4. 

In conclusion, the limitations on the applicability of these results should 
emphasized. The tests were run with fluids of low viscosity in relatively 
rge tanks. Ripples, pendulations in the level of the interface, and the prox- 
ity of walls or other solid obstacles can be expected to increase the gap 
at is necessary to prevent formation of a waterspout. Thus, the results 
esented represent the minimum spacing required. It is evident that the 
ue of y/D of 2 sometimes used in design is not sufficiently great to in- 

e against the formation of a waterspout at high rates of siphonage. 


= 
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Fig.l. Sketch of 
experimental installation 


/0 Saline warerspout 
in tresh water 
e D=0045 fr 
°0 D=0/25 ft 


Worerspout in air 
° D=0045 ft 
e D=Q/25 ft 
2202 94-77 


100 


fig.2. Composite lot of measured data 
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Fig. 3. Surface profile of waterspout in air. 


Fig. 4. Saline waterspout (dyed) in fresh water. 
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PENETRATION OF A JET INTO A COUNTERFLOW 


J. Arendt,5 H. A. Babcock,® and J. C. Schuster? 
(H. Rouse,! M, ASCE, Adviser) 


If a fluid jet is directed upstream against a uniform flow of the same fluid, 
it is to be presumed that the jet will penetrate some distance into the oncom- 
ing flow and then be turned back upon itself and carried downstream, as 
shown schematically in Fig. 5. The degree of penetration—in fact, the geome- 
try and kinematics of the entire flow pattern—could be expected to depend as 
a first approximation upon the diameter and velocity of the jet and the veloci- 
ty of the counterflow; that is, 


f(x, Pe doivs Vy Yo) =-0 
which may be reduced to the non-dimensional form 


oi(9>5>% > 3) =0 


According to this relationship, two jets will yield similar patterns of flow 
(i.e., similar distributions of v/vg) only if the ratio vj/vo is the same for 
both, under which circumstances the distance upstream to the point of stag- 
nation (or the diameter of the enclosed jet flow, or any other characteristic 
length) will be a specific number of jet diameters. However, available infor- 
mation on submerged jets8 would lead one to conclude that approximate simi- 
larity may be attained for jets of sufficiently small diameter and high velocity 
through the control of a single jet variable rather than two together. This 
variable is the momentum flux: 


= =lag?py2 
a ae 4nd PY, 
Under these circumstances, 


fo(x,1T,V, Vg .p,M) = 0 


5. Assistant Professor, Civil Engineering Department, Fenn College, 
Cleveland, Ohio. 

6. Assistant Professor, Civil Engineering Department, Colorado School of 
Mines, Golden, Colorado. 

7. Hydraulic Engineer, Bureau of Reclamation, Denver, Colorado. 

8. Albertson, M. L., Dai, Y. B., Jensen, R. A., and Rouse, H., “Diffusion of 
Submerged Jets,” Transactions A.S.C.E., vol. 115, 1950. 
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whence, in non-dimensional terms, 
x ad Vv 
eas 5») eae 


in which m=/M/o Io = 


Experiments to permit an evaluation of this functional relationship were 
performed in an air tunnel having a 3x3-foot cross section and a velocity 
range from 15 to 45 feet per second. Measuring equipment consisted of a 
differential water manometer read by vernier to 0.001 foot, a hypodermic- 
needle stagnation tube, and calibrated Pitot cylinders with pairs of holes at 
30° and 180° from each other for determining the velocity magnitude and 
direction. The jet was produced by discharging compressed air from a cop- 
per tube of 0.015-foot inside diameter mounted at the axis of the tunnel test 
section approximately midway along its 12-foot length. 

Preliminary observations of the flow pattern were made by means of 
threads giving a visual indication of the distance that the jet penetrated into 
the counterflow and the radial extent of the zone of mixing. Once a sufficient 
length of penetration relative to the initial diameter of the jet had been estab- 
lished, measurements of the velocity distribution were made at a series of 
cross sections. Whereas a proper investigation would require the precise 
determination of both the magnitude and the direction of the velocity vector, 
in the limited time available it was possible to obtain only the magnitude. 
Even this, however, permitted an approximate evaluation of the flow pattern 
by the following technique. 

Assuming as a first approximation that v;, > v, diagrams are prepared of 


| 2mr vx asa function of r, as in Fig. 6. Since Q=27/rvx dr, use of a 


_planimeter then permits the preparation of diagrams for Q asa function of 


: r. Now a stream tube (which is formed of streamlines) represents by defi- 


nition a surface through which no fluid passes. It follows that Q must be the 


same at successive cross sections of the same stream tube, and as well that 
the incremental rate of flow AQ between coaxial tubes must also be constant, 
It is hence possible to approximate the location of streamlines in an axial 
plane by noting the values of r for the same flow rates 0, 4Q, 2AQ, etc., 
at successive longitudinal distances x. From the direction of flow indicated 
by these streamlines, a second approximation may be obtained, which is usu- 
ally sufficient. 

The experimental results were evaluated in this general fashion, the sole 
difference being that non-dimensional rather than dimensional notation was 
used throughout. In other words, the quantity 27rv/mvg was integrated 
with respect to r/m to yield Q/ m2vo , and for the construction of stream- 
lines the increment AQ/m2vg = 0.2 was chosen. The resulting streamline 
configuration is shown in Fig. 7. This may be presumed to indicate, at least 
as a first approximation, the pattern of flow that would be produced by any jet 
of high enough relative velocity and small enough diameter for the latter to 
be inappreciable in comparison with the length of penetration (i.e., Xs/d must 
be greater than perhaps //25). Under such circumstances, it follows that, 


since x,/m = -2.7, m =V M/e/vg , and M = tnd2pv,? ; 


x iM MA 
ee 


enn 
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Although the measurements leading to the pattern shown in Fig. 7, lacking 
as they were in both precision and completeness, must be regarded as ex- 
ploratory only, they point at once to fruitful directions of more extended 
study. Three are readily apparent: first, verification of the assumption that 
m is a widely useful characteristic of jet performance by tests under a con- 
siderable range of conditions; second, determination of the pattern of second- 
ary flow—i.e., turbulence—in the same generalized manner; third, evaluation 
of the effect of viscosity in terms of a Reynolds number such as vgm/y , and 
perhaps eventually a completely analytic study of the pattern for the limiting 
case of laminar flow. 
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SURFACE PROFILES AT A SUBMERGED OVERFALL 


L. F. Ingram,9 R. E. Oltman,10 A, M. ASCE, and H. J. Tracy! 
(E. M. Laursen,12 A, M. ASCE, Adviser) 


Depending upon the interrelationship of a number of variables, the free 
surface at a submerged overfall will assume one of several basically differ- 
ent forms. In addition to the coordinate position, the pertinent variables in- 
clude the velocity and depth of the oncoming flow; the depth of the tailwater in 
the forebay; the channel slope, roughness, and cross-sectional dimensions; 
and the dimensions of the forebay. For the purpose of this exploratory study, 
the channel was made smooth and horizontal and the drop as great as possi- 
ble, changes in boundary alignment were limited to the vertical, and the sur- 
face configuration was observed in the crest vicinity only. The function under 
consideration was thus reduced to one between the following major 
parameters: 


Herein x and y represent horizontal and vertical distances from the over- 
fall crest, yg and yt are the initial and tailwater depths, and F is the 
Froude number V// gyo. 

The experimental investigation of this function was conducted in a glass- 
walled flume 1 foot wide, 2 feet deep, and 10 feet long. Falsework in the 
flume created a 1-foot drop, the horizontal portion of the channel being 2.5 
feet long and having a quarter-circle transition at its upper end. Discharge 
through the flume was measured by a calibrated elbow meter and an air-water 
manometer. The depth of the approaching flow could be controlled by a sluice 
gate at the beginning of the horizontal channel and measured by a point gage 
at a section 0.6 foot (i.e., 4y9) upstream from the drop. The depth of the tail- 
water was controlled by vertical slats at the end of the flume, 7.5 feet beyond 
the drop, and was measured by floor piezometers connected to open manome- ~ 
ter tubes. The water-surface profile was measured with a point gage gradu- 
ated to 0.001 foot and mounted on a carriage sliding on level stainless-steel 
rails. 

The first step in the standard operating procedure was to adjust the control 
valve and sluice gate until the desired rate of discharge and depth of approach 
were obtained. In all runs the approach depth was maintained at 0.15 foot, 


9. Hydraulic Engineer, Waterways Experiment Station, Vicksburg, Missis- 
sippi. 
10. Chairman, Branch Training Committee, U. S. Geological Survey, Wash- 
ington, D.C. 
11. Hydraulic Engineer, U. S. Geological Survey, Atlanta, Georgia. 
12. Research Engineer, Iowa Institute of Hydraulic Research, State Universi- 
ty of Iowa, Iowa City. 
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and Froude numbers of 1, 2, 3, and 4 characterized the different series of 
runs. In each series successively greater tailwaters were established and 
the water-surface profiles measured. Except for the maximum tailwater, 
values of yt/yo which were either positive or negative integers were chosen. 
The maximum tailwater used was that which formed a hydraulic jump just 
upstream from the drop; any slight increase in the tailwater above this value 
resulted in submergence of the sluice gate. Observations were also made of 
the tailwater elevation corresponding to the transitions between a plunging 
and a riding nappe and between a riding nappe and a hydraulic jump; because 
of something akin to a hysterisis effect, these observations were made with 
both a rising and falling tailwater. 

The water-surface profiles as measured are presented in dimensionless 
coordinates in Fig. 8. A more graphic representation of the typical flow pat- 
terns can be obtained from Fig. 9, wherein photographs of the four runs with 
the highest tailwaters at F = 2 in Fig. 8 are shown. Three distinct regimes 
of flow can be seen. At the low tailwaters the nappe plunges below the sur- 
face. The tailwater elevation in this regime of flow affects the nappe prin- 
cipally insofar as it determines the pressure under the nappe: the lower the 
tailwater, the more the nappe is depressed. Above a certain tailwater eleva- 
tion, depending on the Froude number, a distinctly different flow pattern ob- 
tains in which the nappe does not plunge, but rides the surface downstream 
from the drop. As the tailwater elevation is further increased, the surface 
waves which characterize this regime of flow become progressively larger 
until the leading wave begins to break. Above this second critical tailwater 
elevation, a hydraulic jump forms with its toe slightly upstream from the 
drop. A small increase in the tailwater elevation is then sufficient to cause 
the jump to move entirely into the channel. For Froude numbers appreciably 
less than 2 the jump is of the undular rather than the breaking type, and the 
upper regime becomes indistinct from the intermediate. At a Froude number 
approaching unity, at which limit the sluice gate could no longer be used as a 
control, it is obvious that even the undular form of jump could not occur; tail- 
water elevations appreciably above the line of critical depth then resulted in 
complete submergence of the overfall and consequently in flow with a Froude 
number less than unity. 

The relative tailwater elevations at which the transitions between flow 


_ regimes occurred are shown in Fig. 10 as functions of the Froude number. 


Although no instability of the nappe was noted, there is evidently a significant 


_ difference between the values for rising and for falling tailwater. Because 


the flow tended to be stable, the critical elevations for falling tailwater are 


invariably below those for rising tailwater. 
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Fig 6. Measured variation in surface profile with 
Froude number and depth of tailwater 
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Fig. 9. Photographs of surface profiles for F = 2. 
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COEFFICIENT OF CONTRACTION FOR A SUBMERGED JET 


D. D. Curtis,13 M. ASCE, J. E. Martinez,14 J. M. ASCE, and V. Vazquez15 
(H. Rouse,l M. ASCE, Adviser) 


In computing the rate of discharge for a submerged orifice, it is custom- 
ary to assume that the coefficient of jet contraction is the same as it would 
be if the jet were free. This assumption apparently yields satisfactory re- 
sults, else it would not have been made so consistently for the past century 
or more, It is not conceivable, however, that the flow patterns of the free 
and submerged jets would be identical or even closely similar. The present 


_ experiments were hence undertaken to compare both the geometric and the 


_ Kinematic characteristics of a free and a submerged jet issuing from the 


same orifice under the same differential head. 


The tests were conducted in a glass-walled flume having a length of 25 


feet, a width of 2.5 feet, and a depth of 3 feet. A 2-inch sharp-edged orifice 
was placed 1.25 feet above the floor of the flume in the middle of a vertical 


steel plate separating the flume from a 6-foot stilling tank at its upstream 
end, Closure of the flume at its downstream end by a plate 2 feet high pro- 


vided not only a skimming weir for maintaining a constant depth of submer- 


gence but also a volumetric tank for preliminary determination of the head- 


discharge relationships. A gage carriage on rails, point and hook gages 
reading by vernier to 0.001 foot, a differential manometer with the same pre- 
cision of indication, a hypodermic-needle stagnation tube, and a stop watch 
were the measuring instruments used. 

Figure 11 shows the results of the calibration runs for both the free and 
the submerged jet, from which it will be seen that essentially the same coef- 
ficient of discharge—and hence essentially the same coefficient of contrac- 
tion—applies to both cases under the same differential head. The actual 
contraction of the free jet under a head of 2 feet was determined by measur- 
‘ing the surface profile in the vertical plane, computing the elevation of the jet 
axis therefrom, and eliminating the effect of gravitational deflection by sub- 
tracting from the measured profile elevations the corresponding deviation of 
the axis from the horizontal. The superposed results are shown in non- 
dimensional form in Fig. 12, together with computed values for the case of 
irrotational flow.16 Although the profiles differ somewhat in form, the limit- 
ing ordinates are very nearly the same, from which it is found that 


Ce = (0.78)2 = 0.61. 


13. Professor and Head, Department of Mechanics, Clemson Agricultural 
College, Clemson, South Carolina. 

14. Assistant Professor, University of New Mexico, Albuquerque, New Mexico. 

15. Assistant Professor, Monterrey Institute of Technology, Monterrey, N. L., 


Mexico. 
16. Rouse, H., and Abul-Fetouh, A. H., “Characteristics of Irrotational Flow 


through Axially Symmetric Orifices,” Journal of Applied Mechanics, Vol. 
17, No. 4, 1950. ; 
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From the stagnation-tube traverses of the submerged jet a family of 
velocity-distribution curves, selected ones of which are shown in Fig. 13, was 
obtained. From this sequence it is seen that the shear around the periphery 
of the jet both decelerates the flow within the jet and brings into motion the 
fluid immediately surrounding it; the region in which this effect occurs 
necessarily expands steadily with distance from the orifice. The extent of 
the jet contraction must, as a result, be defined rather arbitrarily. Now the 
most logical designation of the jet profile is that within which the rate of flow 
is the same at successive sections. The quantity 27rv at each section was 
hence plotted against the radius r , the integral Q = 27 rf rvdr of the area 
enclosed by each curve thereafter being plotted against r. From each of 
these curves the magnitude of r for which Q equals Q, the rate of orifice 
discharge, could readily be determined. A series of the resulting values will 
be seen in Fig. 14, together with the profile of the free jet from Fig. 12. 

Comparison of the trend of the points with that of the curve will show that 
the profile of the submerged jet attains minimum proportions approximately 
a half diameter from the orifice but then gradually expands, whereas the free 
jet—although it theoretically approaches its limiting diameter asymptotically— 
actually undergoes little further change beyond a half diameter from the ori- 
fice. The minimum dimension of the submerged-jet profile, moreover, is 
somewhat larger than that of the free jet, the coefficient of contraction being 
correspondingly higher—i.e., Cc = (0.79)2 = 0.63. Since the discharge coef- 
ficients of the two jets have already been seen to be essentially the same, one 
is forced to conclude that the rate of flow through the orifice is essentially 
independent of the tangential forces exerted upon the jet by the fluid body into 
which it emerges. In other words, the action of the turbulent shear in causing 
the submerged jet to expand has no appreciable effect upon the pressure field 
that controls the rate of efflux. This is also evident from the experimental 
indication in Fig. 13 that the velocity (and hence the pressure) in the irrota- 
tional core of the jet remains unaffected by the shear in the mixing zone. 


\ 
vi Soe elScet 


Q 
(cfs) 
| 06 


ROUSE 1038-19 


One Q4 06 08 10 
Ah (1) 
fig. //. Results of calibration runs 


Measurement 


— Computation ‘ee aa seciee 
10 


Fig./2. Comparison of measured and 
computed profiles of free jer 


1038-20 HY 4 August, 19 


“0 002 004 006 008 O10 a2 Or 
r (ty) 
/719./3 Variation of velocity distribution with distance 
along axis of submerged jer 


= ° Submerged set 
GC) Ct 
O 
0 OD x LO 
d 


fig./4. Comparison of profiles of free 
and submerged Jers 


‘a ROUSE 1038-21 
| 
| 


MAXIMUM AMPLITUDE OF OSCILLATION OF A 
CYLINDER DUE TO ITS VORTEX WAKE 


) 
|B, Guyton,17 A. M. ASCE, B. W. LeTourneau,18 and R. L. Player, Jr.19 
(L. Landweber,20 Adviser) 


INTRODUCTION 


When a circular cylinder is situated with its axis normal to a stream, ed- 

iés are shed periodically from alternate sides of the cylinder, generating a 
-called Karman vortex trail.21 This gives rise to a fluctuating velocity 

eld and hence also a fluctuating pressure field which exerts a periodic 
teral force on the cylinder. 

_ The frequency of generation of the eddies by a non-vibrating cylinder has 
en determined experimentally.22 It was found, however, that, when the 
linder is allowed to vibrate, the eddy frequency and other properties of the 
ke are strongly affected by the oscillation;23 i.e., the oscillation modifies 
e force which causes it. 

_ An approximate analysis of the vortex-induced vibration of a cylinder, tak- 
Mg into account the interaction between the vibration and the wake, was given 
Steinman,24 who devised an ingenious graphical method for displaying the 

odification of the vibration frequency due to the vibration. The section of 

iS paper entitled “Effect of Vibrations on Wake and on Eddy Frequency” con- 

udes with the statement 


“In these vortex-induced oscillations, there is no ‘catastrophic 
range’ of increasing amplification with unlimited increase of stream 
velocity. The greatest amplitude attainable is little more than the 
diameter, and may occur only in a limited critical range just above 
the reasonance velocity.” 


| 
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), Research Engineer, Iowa Institute of Hydraulic Research, State University 
of Iowa, Iowa City, Iowa. . 
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by the motion of circular cylinders through a fluid,” British Advisory 
Committee for Aeronautics, R. and M. 917, 1924. 

. Landweber, L., “Flow about a pair of adjacent, parallel cylinders normal 
to a stream,” David Taylor Model Basin Report 485, July, 1942. 
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Steinman defines amplitude as the maximum displacement of the center of 
the cylinder from its undisturbed position. However, he furnishes no evi- 
dence for his statement concerning the greatest attainable amplitude. 

The purpose of the present study is to enquire further, by both analytical 
and experimental means, into the question of the maximum amplitude of a 
vortex-induced vibration. 


Properties of the Karman Vortex Trail 
In the range of Reynolds numbers in which laminar separation occurs on 


the cylinder, it is found experimentally23 that the frequency n of vortex 
generation by a vibrating cylinder is given by 


nD _ 0,185 
ee (1) 


where D is the diameter of the cylinder, V is the free-stream velocity, 
and A isthe ratio b/bg of the width of the vortex street with and without 
vibration. If a is the amplitude of the oscillation, we have the experiment: 
ly obtained relation 


i 1+ 1.545 (2) 


The force per unit length F acting on the cylinder is an oscillating quar 
tity whose amplitude is proportional to pV2h , where /p is the density of 
the fluid; or, since b = Abo. and since, for the case of laminar separation, 
is found experimentally that bo = 1.3 D, it follows that 


where § is a constant of the order of unity. 


Analysis of Cylinder Vibration 


The vibration of the cylinder will now be analyzed on the assumption tha 
it is a linear harmonic oscillator satisfying the equation 


e. : 
nS +c H+ ky = F cos 2nnt (4) 


where m is the apparent mass per unit length of the cylinder (actual mass 
plus virtual mass), c is a damping coefficient, and k a stiffness coefficie 
When the vibration is steady, at frequencies near resonance, there results. 
for the maximum amplitude 


2nmnc 


or, from Eqs. (1) and (3), 
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A780 VD" 
va M (6) 
ay = 0.86 hasan 


But, from Eq. (2), 
a, = 0.65 D(Ay - 1) (7) 


) Hence, upon equating the two expressions for ayy in Eqs. (6) and (7), one ob- 
tains a quadratic equation in } y whose solution is 


| 2 =7(2+y2- 5.5422) (8) 


_ The largest permissible value of 5.38 pVD/c is seen to be unity, correspond- 
_ing to which A y = 2 and, from Eq. (7), ay, = 0.65 D. This indicates that the 
_Maximum attainable amplitude is about 2/3 of the diameter. 


Experiments 


| An attempt was made to support experimentally either Steinman’s state- 
‘ment or this simplified theory. Since the theory as developed is based on the 
two-dimensional case of a rod vibrating without deformation, experimental 
confirmation should have been sought with a stiff cylinder supported by 
springs or elastic cantilever beams which would vibrate at right angles to 
the stream without deformation. However, time did not permit the design and 
construction of such an apparatus, and experiments were conducted instead 
on the vibration of a string. Since the string is flexible and deforms ina 
combination of characteristic harmonic shapes when vibrating transversely, 
and in addition is bowed in the direction of flow, the analysis presented above 
is no longer truly applicable, although it may serve as a first approximation. 
No analysis for the case of the string was attempted. 

For the assumed condition that the foregoing theory was applicable to the 
string, resonance would occur when the natural frequency of the string was 
equal to the frequency of generation of the vortices—i.e., 


Wisin et 
0,185 KD 2s Tm 9 n =1, es 3, ecoce 


where T is the tension in the string and s_ its length. 

Attempts were made to secure resonance using both 1/16-inch piano wire 
and 0.054-inch nylon string in a small air tunnel, but without success. This 
may have been due to the fact that the lowest fan speed available resulted in 
an air velocity in excess of 20 feet per second. Resonance was achieved, 
however, by stretching the same nylon string across the 14-inch diameter 
test section of the Institute water tunnel. The velocity was varied from 0 to 
15 feet per second, and the vibration of the string was clearly observed 
through a glass window. The fundamental mode and the first and second har- 
monic shapes, respectively, were clearly distinguishable as the velocity was 


increased. 
f 


7 
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The velocity was measured with a mercury manometer reading differen- 
tial pressure across the bell entrance to the test section. The amplitude of 
vibration, relative to the string diameter was determined by measuring with 
a surveyor’s transit the distance between the upper and lower extreme loci 
of the string as it vibrated. The transit was located on the same level in or- 
der to eliminate extraneous drag vibrations, and the tripod was rubber- 
insulated from the vibrations of the floor. The frequency was measured by 
observing the vibrations with a Strobotac and recording the highest synchro- 
nized frequency. The tension of the string was measured with the water at 
rest by pulling normally to the string center with a spring balance, applying 
a known force, measuring the resultant sag, and calculating the stress by the 
usual methods of elementary statics. It was found difficult to apply a large 
tension to the nylon, because the tied end connections had a tendency to slip. 


RESULTS AND CONCLUSIONS 


The experimental results are presented graphically in Fig. 15, in which 
the dimensionless amplitude a/D and the frequency parameter or the 
Strouhal number nD/V are plotted against the Reynolds number VD/V. The 
Strouhal numbers observed were of reasonable but somewhat high magnitude 
and decreased slightly with the Reynolds number, as would be predicted by 
the relation nD/V = 0.185 and the amplitude data. The amplitude data showed 
considerable scatter, some of which was due to the rudimentary instrumenta- 
tion and some to the inherent nature of the experiment. The amplitude rose 
rapidly from zero at very low speed to a range of a/D from 0.25 to 0.55 at 
the fundamental resonance Reynolds number of 3000 predicted by using a cal- 
culated fundamental natural frequency of the string of 212 cycles per second 
corresponding to the 8.4-pound tension. The maximum observed dimension- 
less amplitude was 0.58 at velocities near fundamental resonance. A mean 
curve is plotted through the data. The results tend to support the foregoing 
simplified theory, since the predicted and observed maximum relative ampli- 
tudes are quite close: 0.65 and 0.58, respectively. 

The visually observed modes of vibration are also indicated on the graph. 
It is to be noted that, once the maximum amplitude was attained, it remained 
relatively constant and independent of the mode of vibration as the velocity 
was increased, This indicates a sort of sustaining resonance, as was graphi- 
cally explained by Steinman.24 The transition between vibratory modes was 
not smooth, but there were ranges of velocity when the vibration was observec 
to “jump” from one mode to another, 
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DEFLECTION OF A JET BY A NORMAL WALL 


Re LE Murray,2° L. E. Nelson, 26 M. ASCE, and J. R. Steven, 27 J, M. ASCE 
(P. G. Hubbard,28 J, M. ASCE, Adviser) 


Although the deflection of a submerged jet by a normal wall is closely re- 
lated to practical problems such as the scouring tendency of a plunging nappe, 
it is of interest primarily as an example of free turbulent diffusion and mo- 
mentum transfer. The investigation of this phenomenon had been preceded at 
the Iowa Institute by a rather complete study of the flow pattern which is 
created by discharging a jet of fluid from an opening in a plane boundary into 
a semi-infinite reservoir of the same fluid.29 In the present case, the object 
was to determine the streamline pattern which results when the flow is de- 
flected by a wall normal to the initial direction of the jet. 

The experimental equipment for the study consisted of two carefully 
rounded nozzles 0.025 and 0.05 foot in diameter which could be mounted, in- 
terchangeably, flush with the surface of an octagonal table 30 inches across. 
Air was supplied to the 6-inch pipe preceding the nozzle by a centrifugal fan, 
and the discharge was controlled by a variable-area intake to the blower. A 
masonite board mounted 2.5 feet above the table provided the normal wall for 
deflection of the jet. Velocities throughout the region of interest were mea- 
sured by means of a hot-wire anemometer mounted on a mechanism for mak- 
ing longitudinal and lateral traverses. The pressure distribution on the nor- 
mal wall was measured by an alcohol-kerosene differential manometer which 
could be connected to a series of piezometers arranged radially outward from 
the stagnation point. 

Because the flow system was axially symmetrical, it could be described 
completely in terms of the velocity v(x,r) or streamlines y(x,r) ina 
meridian plane, x being measured from the plane of the issuing jet and r 
radially outward from its axis. From the available results of the earlier jet 
studies already mentioned, it was known that the determinants of the flow pat- 
tern would be the nozzle diameter D , the momentum flux of the jet | 
M = Qpvj , the fluid density p , and the distance a between the bounding 
planes. In order to reduce the scope of the study to correspond to the time 
available, a was made so large compared to the jet diameter that the jet 
could be considered as a point source of momentum. 

A dimensional analysis of the problem, beginning with the variables 


25. Assistant Professor, College of Engineering, University of Santa Clara, 
Santa Clara, California. 

26. Hydraulic Engineer, U. S. Engineer Office, Kansas City, Missouri. 

27. Assistant Professor, Department of Civil Engineering, College of the a 
of New York, New York, 

28. Research Engineer, Iowa Institute of Hydraulic Research, State University 
of Iowa, Iowa City. 

29. Albertson, M. L., Dai, Y, B., Jensen, R. A., and Rouse, Hunter, “Diffusion 
of Submerged Jets,” Trans. ASCE, Vol. 115, 1950. 
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already defined, can be outlined as follows: 
v= f(x, eae M) 


or, in dimensionless terms, 


Wipe = %(a+ 5) 


An equally useful and more easily visualized method of presenting the re- 
sults is to plot streamlines in an axial plane, using the definition 


ya [vy cdr = Q/en 
0 


where vx is the axial component of velocity and Qy is the discharge 
through a circle of radius r. A pattern of streamlines can be obtained from 
this relationship by evaluating graphically the radii corresponding to equal 
increments of y and then plotting curves of constant y or of the dimen- 
sionless value 


oe 
apa reed) = 022) 


Because the hot-wire anemometer used in measuring mean velocities was 
oriented so that its axis was normal to an axial plane, it indicated the result- 
ant velocity v rather than the axial component vx required in the computa- 
tions for the streamlines. This meant that two or more successive steps had 
to be taken in determining the location of the streamlines in the region where 
they were not parallel to the axis of the jet. The first approximation was ob- 
tained by assuming the measured velocity to be axial, integrating graphically, 
and plotting a streamline pattern on this basis. With 6' as the angle between 
the axis of the jet and the tangent to one of these streamlines, then a second 
graphical integration could be made using v' = v cos 6'. This process con- 
verged rapidly to the correct pattern, and even the second approximation was 
unnecessary except where the streamline curvature was pronounced. 

Near the wall, of course, the streamlines curved rapidly and became prac- 
tically radial where r/(a-x) >3. For this region a modified method of locat- 
ing streamlines was used. If vr represents the radial component of velocity 
and Qg-x the rate of flow through a cylindrical surface of radius r and 
length a - x , then the stream function can also be defined as 


y= Be -f vs(- ax) 


This takes the dimensionless form 
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Where correction of the values for inclination of the streamlines is neces- 
sary, the process is similar to the one already outlined, except that the 
quantity v' = v sin 9' is substituted in the numerator. In the region 

x/a > 0.8, r/a<0.3, both methods were used as independent checks upon the 
accuracy of integration. 

The final result for this study is shown in Fig. 16. Identical curves were 
obtained using each of the two available nozzles, so that the assumption of a 
point source of momentum was justified. In the latter circumstance the 
stream tube enclosing the discharge from the nozzle becomes so small that 
it cannot be shown on such a plot, the flow pattern hence seeming to consist 
almost entirely of streamlines for the fluid entrained by turbulent diffusion. 
The dashed line indicating the nominal boundary of the jet has the slope 1:5, 
which applies to a jet without the normal wall, and which also defines the 
limits in this case quite well. 

As is evident from the diagram, the effect of the wall is felt primarily in 
the region x/a > 0.8. In spite of the concentration of the flow about the axis, 
the measured form of the streamlines in this region is almost identical with 
that computed on the assumption of irrotational flow. For the latter, the 
streamline equation is30 


ay = mr*#(a-x) 


wherein the coefficient m for the present experiments is approximately 


52 (M/p fon : 

Dimensionless values of the pressure measured on the normal wall are 
also plotted in Fig. 16, together with that obtained from the foregoing rela- 
tionship for irrotational flow with a stagnation value of 27 M/a2. Integra- 
tion of this pressure distribution yielded the unexpected result that the total 
force on the wall is approximately 4/3 as great as the initial momentum flux 
of the jet rather than exactly equal to it. This apparent dynamic unbalance 
was probably due to the pressure reduction on the lower boundary accompany- 
ing the radial acceleration of the air drawn in by the jet. Because this nega- 
tive pressure occurred over a relatively large area, its local intensity was 
so small that it could not be measured reliably. 

On the basis of the experimental results of this study, two important con- 
clusions can be drawn concerning the deflection of a jet by a normal wall. 
First, the wall has a negligible effect on the flow in an initial zone having a 
radius of at least one-fifth and a length of roughly four-fifths the boundary 
Spacing. Second, in the region of pronounced streamline curvature near the 
stagnation point on the wall, the flow pattern may be assumed to agree with 
that of irrotational motion over a radial distance of about a/8. 


30. Forchheimer, P., Hydraulik, Leipzig and Berlin, 3rd Ed., 1930, pp. 28 
and 29, 
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INFLUENCE OF SLOPE AND ROUGHNESS ON THE FREE OVERFALL 


J. W. Delleur,31 J. M. ASCE, J. C. I. Dooge,32 A. M. ASCE, and K. W. Gent33 
(E. M. Laursen,!2 A, M. ASCE, Adviser) 


This study is concerned with the depth of flow which occurs at the end of a 
long channel discharging freely over a fall under various conditions of chan- 
nel slope and bed roughness. The case of a horizontal channel ending ina 
free overfall was investigated in detail by Rouse (Civil Engineering, April 
1936). The effect of channel slope was studied by Fathy and Amin (A.S.C.E. 
Proceedings, Separate No. 564, December 1954), A listing of a number of 
publications on various aspects of the free-overfall problem appeared in the 
February 1955 issue of Civil Engineering. 


Momentum Analysis 


By examining flow conditions, assumed to be two-dimensional, between the 
point of overfall and an upstream point x at which the pressure distribution 
may be taken as hydrostatic (see Fig. 17), one can express the momentum 
equation for a direction parallel to the bed in the following form: 


Ait. Meas Po+M, + W.sin @ - F, (1) 


If 6 is taken as positive downward in the direction of flow, the above equation 
will apply to both downward and adverse slopes. Upon evaluating each term 
in turn, one obtains 


U4 
is R =/ pdy = K, 74 cos @ 
ih Kis J 'p dy/(7y,?cos 9/2); 


Ur 
Me = pf V2 dy = Kz pV; ug 


where K2 = [ “ye dy/ (Ve? Ue) 5 E 
R = 7 3 cos ) 


31. Assistant Professor of Civil Engineering, Purdue University, Lafayette, 
Indiana, 


32. Senior Design Engineer, Hydrometric Section, Civil Works Department, 
Dublin, Ireland. 


33. min ee Engineer, H. G. Acres & Co., Ltd., Niagara Falls, Ontario, 
anada, 
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_ since pressure is hydrostatic; 


Yx 
M,. = ef Vidy = Ksp\y’u, 
| where K,; = / V? dy/(Ny? Tap 
W_ sin® = ¥ sin ef ydx = Ky sinO yy L 


L 
| where K, -[ y dx/(y,L) ; and, finally, 


Fs [ tax ‘ if’y Sak = Kyu.S.L 


where S; is the resistance slope at depth y , Sy is that at point x, and 
| Ks = fy S; Ax/(UxS.b)- Equation (1) now becomes 


2 
KF cos 8 + K, OV," Up = 
Saioe 
= 1% 0088 + ypu, + Ky zy,Lsin@-Ksyy.S,L (2) 


Division of each term by y and elimination of velocities by the continuity 


equation Vsy¢ = Vxyx =4 = Bye? cos@ then yields 


y a Pd yx? je” 
K, 7 cos 8 + Kz cos 8 = = cos 8 + K,= cos @ + y.L(K, sin 8-K,S,) (3) 


This equation is quite general, being identical with Eq. (1). 

: If the slope does not exceed 1:10, cos 9 can be taken as unity and sin 0 
written as Sg. Furthermore, if the flow is turbulent, Kg will be slightly 
greater than 1.0 and Kg will be between 1.0 and Kg , because the flow is 
convergent. Since these terms are on opposite sides of the equation, one can 
take Ko = K3 = 1.0 without appreciable error. These two changes give 


Ys? Y,? ies he A 
Sie ary = z + vr. ts Wo, So~Ks = “ 


If So and S, are not too large, it is reasonable to assume (subject to ex- 
perimental verification) that the last term in Eq. (4) can be neglected to give 


2 3 2 3 5 
2 UF 2 Yx 


This is applicable to turbulent flow on both positive and adverse slopes under 
various roughness conditions. For steep slopes of sufficient length the depth 


en. 
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of the approaching flow is equal to the normal depth, so that | 


Ut pales Uo* 2 ie 
dar nays BMA 


which can be written as 


Up /Ue = Yo/Ye ax (So/5.)? (6) 
2K Ae 2+ (oA)? — 2(S0/S,)+ 1 


or 


L1- KiCur/ue)*] Ut /Yo x 2%) = 2 90 z “Ape (7) 
1+ us /Yo Yo 


For mild or adverse slopes yx can be taken as equal to yc , so that Eq. (5) 
becomes 


aut = 2x, (Ye) (8) 


Experimental Equipment and Procedure 


The experimental investigation was carried out in a glass-walled tilting 
flume 30 feet long and 2 feet wide. Discharges were measured by means of 
an orifice meter and were selected to give simple numerical values of the 
critical depth for parallel flow. In the preliminary tests seven values of the 
critical depth ranging from 0.1 to 0.4 foot were used, while the final tests 
were confined to discharges corresponding to 0.2 and 0.4 foot. Ten slopes in 
the range +0.03 to -0.02 were tried. Two roughness conditions were investi- 
gated: (a) the relatively smooth stainless-steel bottom of the flume, and 
(b) an artificial roughness consisting of 1/4x1/4-inch brass bars fixed trans- 
versely along the bottom of the flume. Previous tests by Walter Rand of the 
Institute staff had indicated that with this artificial roughness in position the 
virtual bottom of the channel was at the mid-height of the roughness bars. 
Accordingly, the crest of the overfall was arranged as in Fig. 18, so that the 
crest level would be in line with the virtual bottom of the approach channel. 

The standard experimental procedure was as follows. The discharge was 
set to give the required value of the critical depth for parallel flow and kept 
constant while the slope was varied throughout its range. For the horizontal 
case and for the maximum positive and negative slopes, enough point-gage 
readings were taken to determine the complete surface profile, For the in-_ 
termediate slopes, enough readings were taken to determine yg , L , and the 
slope of the water surface at the overfall crest and at the point x. 

The scatter in the results obtained seemed to be within the limits to be ex- 
pected for the experimental set-up involved. The fluctuations in the air- 
water manometer used in measuring the discharge could easily result in 
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an error of 2% or more in the discharge and hence a corresponding error in 

the value of the critical depth and critical slope. For flows near the criti- 

| cal and for the steeper slopes the water surface was quite disturbed and the 
depth as determined by point-gage readings could be in error by 1 or 2%. 


DISCUSSION OF RESULTS 


The experimental results can be first used to justify the dropping of the 
last term in Eq. (4) - namely, L y,(K4Sg-KsSx). From the detailed profiles 
it was found that, as yg increased towards yx , Kq4 varied from 0.90 to 0.98 
and Ks varied from 1.25 to 1.05. It was noted that the following approxima- 
| tions based on a parabolic profile gave good results: 


K5= JE ax/y? zs 2 + (y/us)* 


L/y,* a 


| _ kudx us/y 
| Secor ae Rae 


For steep slopes the term under examination is L yx Sg (Kq - Ks), which 
varies from -0.015 x2 to -0.005 x" For mild slopes its value varies from 
/—-0.02 x" to -0.04 Yx2 and in sixteen of the twenty cases measured lies be- 


tween -0.025 Yx2 and -0.035 Yx2. The term is clearly of the same order (and 
opposite sign) as the approximation involved in taking K3 = 1 and may be 

‘neglected in a first approximation to the solution. 

| The critical slope Sc for two-dimensional flow is given by 


In evaluating the critical slope in the present tests, the value for the resist- 
ance coefficient f for the rough case was taken from an equation derived by 
Walter Rand for the same artificial roughness placed in the same flume, 
while for the smooth case C was taken from a formula proposed by R. W. 
Powell. The two formulas are 


. 0 Yo + A 
Ss th: = — 


_ The values of yg measured in the final series of tests are plotted in 

Fig. 19 in terms of yf/yc and Sg/S¢ . It will be seen that the points repre- 
‘senting widely varying conditions of roughness and slope define a single line. 
This is to be expected in the case of steep slopes, as the plotting is in terms 
of the factors occurring in Eq. (6). The values for mild and adverse slopes 
also plot fairly well about a line, but the small differences involved make it 


‘ 


j 
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difficult to determine whether the scatter which occurs is systematic or not. 
In the case of adverse slopes, less scatter is obtained if yg/yc is plotted 
against Sg +Sc (which will be negative when Sg < - Sc), but further experi- 
mentation would be required to settle this question. Accordingly, the results 
for steep, mild, and adverse slopes are all plotted in Fig. 19 in terms of 
yt/yce and S/S. . 

The corresponding relationship based on Eqs. (6) and (8) cannot be deter- 
mined unless values for the pressure variable Kj are known or assumed, 
The experiments indicated that Kj decreased in value as the slope passed 
from adverse to horizontal and from horizontal to mild and then to steep. A 
good correspondence with the experimental points (note the plotted curve) is 
obtained if the following values of Kj are assumed: 


S/S, < -5 K, = 0.60 
5 < S/S,< 1 K, =0.50 + -/ 1-S,/S, /8 
a S/S, K, = 0.30 


However, the validity of assuming these values must still be checked by 
direct pressure measurements. 
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A NEW DEVELOPMENT IN FLOW MEASUREMENT: 
| : THE DALL FLOW TUBE 


Andre L. Jorissen,* M. ASCE 
(Proc. Paper 1039) 


: 
: 


Several new types of flow-metering devices of the pressure-difference or 
Larential type for pipe flow measurements have been developed in recent 
rs. 

| Of particular interest to civil engineers is the Dall Flow Tube which, as 
by Fig. 1, is a modified Venturi Tube. Originally built and studied in 

gland! its characteristics were first discussed before an American audi- 

e by I. O. Miner of Builders-Providence, Inc.2 

Essential differences between the Dall Flow Tube and the classical Venturi 

be are: 

1) The steepness of both the converging cone and of the diverging or 

pressure-recovery cone, the former having an included angle of 

| approximately 40° and the latter an included angle of approximately 15°, 

| 2) The presence of a dam at (a) and the characteristic location of the up- 
stream pressure connection, similar to the type commonly used in 

__ Europe and known as “corner taps.” 

3) The sharp edges at (b), (c), (d) and (e). 

4) The slot between (c) and (d), through which the downstream or throat 

pressure connection is made. 

5) The sudden enlargement at (f). By truncating the recovery cone and 

because of the steepness of the converging and recovery cones, a very 

compact instrument is obtained, the length of the whole device being of 

the order of two pipe diameters only. 


Experimental data indicate that the differential pressure produced by a 

ll Flow Tube is higher than that of a conventional nozzle or Venturi Tube. 
ie ratio of Dall Flow Tube differential to Venturi Tube differential varies 
pm 1,91 for a diameter ratio B of 0.53 to 2.43 for B = 0.8. 

——eEE——EEEES 

te: Discussion open until January 1, 1957. Paper 1039 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
82, No. HY 4, August, 1956. 

‘Head, Dept. of Hydr. and Hydr. Eng., Cornell Univ., Ithaca, N. Y. 

By H. E. Dall, Research Engineer for George Kent, Ltd. 

I, O..Miner: “The Dall Flow Tube,” Trans. A.S.M.E., 78, no. 3, April 1956. 
B is the ratio of throat diameter to pipe diameter: B = d/D. 


1039-1 


Peper ts 


1039-2 HY 4 August, 195 


An unexpected characteristic of the Dall Flow Tube is that the nonrecove: 
able pressure head, or head loss of the instrument is extremely low. Fig. 2 
gives a comparison, established by Miner, between the head loss, expressed 
in per cent of the differential, for the Dall Flow Tube and two types of Ven- 
turi Tubes. In the standard Herschel-type long form Venturi Tube, the in- 
cluded angle of the recovery cone is 5°5', In the short form Venturi Tube, 
this angle is appreciably larger. Furthermore, the recovery cone is trun- 
cated to obtain a more compact instrument. The solid lines on Fig. 2 indi- 
cate the head loss in per cent of differential as a function of the diameter 
ratio B. The dotted line gives the loss through a Dall Flow Tube producing 
the same differential as that of a Venturi Tube with the value of 8 given by 
the abcissa. 

The advisability of expressing the head loss in per cent of the differentia 
may appear questionable because of the high value of the differential itself. 
Another comparison, based on the ratio of head loss in Dall Flow Tubes to 
head loss in Herschel type Venturi Tubes, long form for equal flow rates (at 
pipe Reynolds Numbers larger than 350,000) indicates that, for a given value 
of B, the head loss for the Dall Flow Tube is smaller when the diameter 
ratio is greater than 0.552. No satisfactory explanation has been given yet 
of the low pressure loss of the Dall Flow Tube. Several studies are in pro- 
gress, trying to obtain a better understanding of the fluid mechanics princi- 
ples involved in the flow through a Dall Flow Tube. 


Dall Flow Tube Discharge Coefficients 


The relationship between rate of discharge and differential head is ex- 
pressed in the most general case of compressible fluids by the classical 
equation for inferential flow meters: 


CY: 2g.h 
q= ——~!* 


Vy 1-64 


in which 
q = volumetric rate of discharge 
2 
a = throat diameter = a 


C = coefficient of discharge 

Y = expansion factor 

B = diameter ratio = d/D 

g = acceleration of gravity 

h = differential pressure head 


It is well known that the coefficient of discharge C is a function of the 
type and proportions of the instrument and of flow characteristics which co! 
bine in the expression of the pipe Reynolds number: 


vD 
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V = average velocity of flow in the pipe 
D = pipe diameter 


v = coefficient of kinematic viscosity of the fluid 


The expansion factor Y takes into account the compressibility of the 
fluid. It is a function of the diameter ratio 8, of the ratio of throat to inlet 
pressure and of the adiabatic exponent k. For incompressible fluids, Y has 

the value of unity. 
| Values of C, as obtained by Miner from his experiments are given by the 
| solid line on Fig. 3. Those values are for Reynolds numbers larger than 
} 350,000. 

The author has had the opportunity to calibrate several Dall Flow Tubes in 
the Hydraulic Laboratory at Cornell University. One of them was of large 
size: 24.02 in. x 17.544 in. (8 = 0.7304). This tube was calibrated by mea- 
suring flows between 5.62 and 20.50 cfs volumetrically in a standpipe 6 ft. in 

\diameter and 55 ft. high providing a tank of 1500 ft3 useful capacity. A suit- 
jable diverter equipped with an electrical timing device was used. Differen- 
tial pressure heads ranged from 0.282 to 3.80 ft. of water. The results of 
that calibration are given on Fig. 4. 

The curve shows the characteristic trend indicated by Miner, the coeffi- 
jcient of discharge decreasing for increasing Reynolds number and reaching 

a constant value: Ce at high Reynolds numbers. 

The constant value of the coefficient of discharge is Ce = 0.661, which 
agrees very well with Miner’s data from Fig. 3. 

Other calibrations were part of a research project on inferential meters 
conducted in the Hydraulic Laboratory at Cornell University with the spon- 
sorship of Builders-Providence, Inc. The tests were all conducted with 8 in. 

|pipe size and the flows measured volumetrically. The results of those tests 
are summarized in Table I giving the dimensions of the Dall Flow Tubes 
tested, the value of Reynolds number at which the coefficient of discharge be- 
comes constant: (Rp)n and the constant value of the coefficient of discharge: 
Cc. The last three lines apply to the instrument known as the Dall Insert 
|Nozzle which is a Dall Flow Tube inserted between flanges in the pipe line. 

| Jt should be noted that the values of (Rp)n differ somewhat from the value 
‘of 350,000 indicated by Miner. The difference, however, is small except for 
‘two of the largest diameter ratios. 

__ A comparison of the experimental values of Cc with the data from Fig. 3 
lindicates a very satisfactory agreement except for the Dall Insert Nozzle with 
ithe largest diameter ratio. Part of this difference may be accounted for by 
ithe fact that inferential flow meters of large diameter ratio are more sensi- 
itive to the conditions of installation than those of small diameter ratio. 


| 


Effect of Two Short-Radius Elbows in Orthogonal Planes 


Miner has investigated the effect of various upstream disturbances on the 
ehavior of the Dall Flow Tube, namely the effects of a single elbow, of a 
gate valve partly opened, and of an increaser or a decreaser. 

_ The author studied the effect of two short-radius elbows in orthogonal 
lanes. 
The test run was an experimental pipe, 8 in. in diameter, under a static 

d of 80 ft. 


_ Tests were conducted both without and with straightening vanes at the 
3 ; 
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upstream end of the approach pipe, immediately after the two elbows. These © 
straightening vanes consisted of four radial baffles approximately 16 in. (2D) 
long. 

The short-radius elbows themselves were preceded by sufficient lengths of 
straight pipe to insure that they alone would affect the behavior of the instru- 
ment tested. 

The results of the tests are summarized in Table II. La is the length of 
approach pipe between the two elbows and the tube tested, D_ the pipe diame- 
ter, C the coefficient of discharge and a the difference (in per cent) between 
the coefficient of discharge for a particular condition and that obtained for the 
same tube with a length of approach pipe 27 ft. 7/8 in. (41D) long. 

It should be noted that any value of @ smaller than 0.2 per cent is insig- 
nificant, being smaller than the possible experimental error. 

The study indicates that the effect of two short radius elbows in orthogonal 
planes on the performance of Dall Flow Tubes is only slightly larger than 
their effect on standard Venturi tubes. The disturbance appears more signifi- 
cant for larger diameter ratios. The effect of straightening vanes is to raise 
the coefficient. 

Comparison of the values of the coefficient of discharge obtained without 
and with straightening vanes gives the following differences: 


0.15 per cent for the 8 in. x 4 in, Venturi Tube 
0.47 per cent for the 8 in. x 4.4 in. Dall Flow Tube 
0.31 per cent for the 8 in. x 5 in. Venturi Tube 
0.90 per cent for the 8 in. x 5.2 in. Dall Flow Tube 
0.73 per cent for the 8 in. x 6 in. Venturi Tube 


Tests are now under way to study the effect of two elbows ona Dall Flow 
Tube of very large diameter ratio. 


Expansion Factor Tests 


A determination of the value of the expansion factor Y is essential to the 
use of Dall Flow Tubes with compressible fluids. Since no theoretical analy- 
sis is available to this effect, the values of the expansion factor must be ob- 
tained experimentally. This was done by the author, using air flows. The 
procedure was as follows: 

Various Dall Flow Tubes were first calibrated in the Hydraulic Labora- 
tory, using water. The Dall Flow Tubes were installed in series with Venturi 
Tubes of known characteristics. Suitable lengths of approach pipes were pro- 
vided both for the Venturi Tubes and for the Dall Flow Tubes to eliminate the 
effect of approach conditions. After the water calibration, each installation 
was installed in an air compressor line. Since the expansion factors for 
Venturi Tubes are known4 the Dall Flow Tube expansion factors could be de-_ 
termined by comparison. The results are given in Fig. 5. Values of the ex- 
pansion factor Y are plotted as a function of x/k. 


x = ratio of differential pressure to absolute inlet static pressure 
k = adiabatic exponent 


In each case, the Dall Flow Tube expansion factor may be compared with that 


4, A.S.M.E. Power Test Code, Part 5, Chapter 4. 
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of a concentric orifice of equal diameter ratio.5 The conclusion is that under 
the conditions investigated the values of the expansion factors for Dall Flow 
Tubes and concentric orifices are quite comparable. 


CONCLUSIONS 


The Dall Flow Tube, as now constructed, thus offers the two-fold advantage 
of a great compactness and a low head loss. It is only slightly more sensitive 
to the conditions of installation than standard Venturi Tubes. 

For Dall Flow Tubes with throat to inlet diameter ratios 8 not more than 
0.75 (usual limit for Venturi Tubes), the coefficient of discharge can be pre- 
dicted within 1%. For larger values of 8, direct calibration under upstream 
installation conditions is recommended, unless larger tolerances on coeffi- 
cient value are acceptable. 


Ae a 
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M.E. Power Test Code, Part 5, Chapter 4. 


HY 4 August, 1956 


1039-6 


ee ee i a ee eee eee eee 


*queqysuod seul0deq 


0 e31eyoSTp Jo 4yUeTOTJJe0o ey. YOTYm 4e ZequmU sptousey edtd so enTeA (2) 


000 £0S*7 
000 £007 
000 ‘0017 
000 £004 
000 ‘0047 
000 £00€ 
000‘ 00d 


(€) (Z) 
u dy 


ezreyostp Jo syuetotyyeog 


I oTqeL 


*8TZZON YtOsUT TTed (T) 


(T)00°s 
(T)o0°s 


(T)oo°s 


00°38 
00°38 
00°8 
Z0° 72 


1039-7 


JORISSEN 


ASCE 


7786°0 4569°0 g066°0 €€TL°O 


9786°0 6969°0 4066°0 9ETL°O 


Te86°0 0S69°0 0166°0 LeTl*o 
9LL6°0 T889°0O ¢886°0 0604°0 


€L16°0 OOTL*°O 9€86°0 T? | 8/L lz 
086°0 90TL*°O T786°0 82% |8/T-7 -6T 
8946°0 €0TL*O z€86°0 S*StIV/E-9 -ot| 
876°0 9602°0 €€86°0 $°9 17 


736°0 (e) (e) 
soueA FUTUSPYSZTVIZS YNouTM s4S9], 


i ae Me > eS Uy-3 
eqny Frn}Ue, eqny, MOTI Tred 
‘UR gx "UTE | ‘Ut 79 x *UT gs 


soueTd TeucFoyuzAO UT SMOQTS SsNn{pei—ploys omy JO JOTI 
II °Tq', 


August, 1956 


HY 4 


ral 


(a 


1039-8 


ASCE 


HEAD LOSS—PERCENT OF DIFFERENTIAL 


JORISSEN 1039-9 
HERSCHEL VENTURI 
TUBE SHORT FORM 


HERSCHEL VENTURI 
TUBE LONG FORM 


DALL FLOW TUBE WITH DIAMETER 
RATIO THE SAME AS VENTURI 
TUBE DIAMETER RATIO 


DALL FLOW TUBE WITH DIFFERENTIAL 
THE SAME AS VENTURI TUBE 
DIFFERENTIAL Rp > 350,000 


0.4 0.5 0.6 0.7 0.8 0.9 
VENTURI TUBE DIAMETER RATIO 8 


HEAD LOSS COMPARISON 
FIG. 2 


1039-10 HY 4 August, 1956 


DIAMETER RATIO 
FIG. 3 
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0 002 004 006 008 O10 O12 O14 0 O02 004 O06 O08 O10 O12 al4 OE 
x/k x/k 
8X 4.397" DALL FLOW TUBE 870.5499 8"X 4.078" DALL INSERT NOZZLE §=0.5096 


0.02 004 006 008 O10 O12 0.14 © 002 004 006 008 O10 O12 0.14 
x/k x/k 
8"X 5.217" DALL FLOW TUBE 8*0.6527 8"X 5.217" DALL INSERT NOZZLE B+0.6521 


O 002 004 006 008 O10 O12 ale © 002 0.04 006 0.08 O10 O12 0.14 
x/k x/k 
8"X6.336" DALL FLOW TUBE 80.7924 8"X 6338" DALL INSERT NOZZLE 80.7922 


——— MOST PROBABLE LINE 
ea EXPANSION FACTOR FOR CONCENTRIC ORIFICE 


DALL FLOW TUBE EXPANSION FACTOR 
FIG. 5 
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A proper evaluation of the role of the flood plain in river planning involves 
a consideration of all aspects of water resources. Such river problems as 
transportation, power, recreation, water supply, sewage disposal, pollution, 
|flood control and sedimentation affect the welfare of the nation. A study of 
all the problems connected with and involving streams is beyond the scope of 
this paper. Instead, only one aspect of the flood control problem will be dis- 
cussed, namely, the use and control of the flood plains through the medium of 
izoning and the problems encountered. 
_ Regulating or prohibiting the use of the flood plain for dwellings or other 
property which would be seriously damaged by flooding has been proposed at 
various times since 1900. With the acceleration of house construction in the 
late 1940’s, the use of zoning to control the flood plain has increased in the 
|United States. The concept of flood plain zoning has its greatest use as a 
leeventive measure in areas which are undeveloped and sparsely developed. 
The presence of buildings and other structures of various types may make it 
more economically feasible to provide protection to the area rather than to 
evacuate the area. The nature of the zoning procedure takes many forms and 
jis administered and enforced by numerous governmental units and agencies. 
The range of flood plain zoning varies from control of encroachments on the 
istream itself to the prohibition of the use of any portion of the flood plain for 
any activity or structure which will adversely affect flood flows. In some 
jcases, where zoning powers were believed to be insufficient to restrict the 
juse of the flood plain, the areas were purchased as public lands and developed 
for parks and other recreational facilities. 
- Encroachment on the flood plains of rivers and streams has occurred since 
ithe beginning of civilization. The earliest settlements, requiring water for a 
multiplicity of uses, such as transportation, potable water, protection and 
waste disposal, grew up along the banks of rivers, often on the flood plains. (1) 
These settlements established economic values and social roots which have 
thstood, in most instances, the destructive forces of floods. But floods are 
eriodic and persistent, and necessitate action by man to minimize or prevent 


1 


f lote: Discussion open until January 1, 1957. Paper 1040 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
82, No. HY 4, August, 1956. 

| * Hydraulic Engineer, Bergen County Engineering Department, Hackensack, 
New Jersey. 


? 1040-1 


1040-2 HY 4 August, 1956 | 


the damages caused by excessive storm runoff. Sometimes this action has 
been directed towards protecting these values and roots without determining 
if the protection, when all aspects are considered, is warranted. The net re- 
sult is a continuous contest of man opposing the forces of nature and in many 
instances the results are far from satisfactory. 

The expense involved, initially and subsequently, of providing the protec- 
tion against floods has reached the point where the wholesale evacuation of 
the flood plain is being considered as the most economical solution to the 
man-made problem of flood damages. Some areas, benefiting from the ex- 
periences of their neighbors, are attempting to forestall future flood prob- 
lems by maintaining as much of the flood plain as is economically feasible 
for the passage of flood waters. 


The Role of the Flood Plain 


The flood plain plays an important part in the development of a society. 
The environmental features of the flood plain lend themselves to manifold 
uses which have a direct bearing on the growth of cities and the development 
of the area. The proper uses of the flood plain varies from area to area and 
depends on the climate, geography and natural resources, as well as on the 
existing and potential development of the non-flooded area. The solution of 
stream problems is usually directed towards the maximum utilization of the 
flood plain. This may involve the inundation of agricultural lands by storage 
reservoirs on one part of the stream to protect highly developed industrial 
and residential areas situated on the flood plains of another portion of the 
stream. The social, psychological and economical problems arising from 
such a solution are exemplified in the controversy which resulted from 
the proposed Tuttle Creek dam on the Big Blue River in Kansas to protect 
the cities of Topeka and Kansas City. 

Planning, in its best sense, entails full knowledge of the problem, its 
causes, its symptoms and it cures. So with river planning, it is necessary to 
understand the behavior of streams and the characteristic features of their 
drainage basins in order to plan effectively for the maximum utilization and 
protection of the water resources available. In addition, procedures must be 
established to minimize the dangers from floods, an excess of water within a 
relatively short period of time. ; 

A stream is a constantly moving body of water providing drainage for all 
parts of the watershed. Storm runoff, from roofs, lawns, streets, farms and 
other surfaces, eventually finds its way to a natural drainage ditch or man- 
made storm sewer and then into a stream. Past experience has repeatedly 
shown that isolated sections of a stream cannot be altered without regard for 
adjacent reaches of the stream. A drainage basin is a closely correlated 
natural unit and should be studied and treated as a unit. The land along any 
stream may be divided into three zones:(2) 


1) The normal stream channel. 


2) The flood plain, which is usually dry but occasionally flooded. 
3) The high land, which is never flooded. 


In general, the stream channel is encroached upon less frequently than the 
flood plain. Many of the stream channel encroachments are bridges neces- 
sary to provide continuous and uninterrupted transportation. These bridges, 
if properly designed from a hydraulic viewpoint, nullify, to a considerable 
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xtent, any adverse effect on flood flows. Most of the remaining stream 
hannel encroachments occur along banks of streams. Structures, including 

ls, which extend into the existing stream channel frequently restrict the 
jtream to the point where higher flood stages result. The economic and 
hysical obstacles involved in constructing structures adjacent to streams 
ave deterred a greater number of encroachments on streams. The flood 
lain, on the other hand, does not have these obstacles unless it is swampy. 
ince the flood plain is usually dry, sometimes for years, there is consider- 
ble encroachment on this area. The damage to these encroachments during 

es of flood account for a major portion of the total flood damages. There- 

re, if future encroachment on the flood plain is not controlled, flood dam- 
ges will increase rapidly and the cost of providing protection will also in- 
jrease. 

In many places, encroachment in the flood plain has been encouraged by a 
gal doctrine which permits man to invade not only the flood plain but even 
e stream channel itself, despite recurring and increasing flood losses. 

Vhile some of this encroachment has been made in ignorance of the flood 
ard, most of it has resulted from attempting to utilize the resources of 

e terrain, soil and water.(3) Farmers found easy and rich soil to till, 

nufacturers found cheap land, convenient transportation and power, rail- 

oad companies found easy grades on which to lay their tracks. Following 

ne usual tendency, businesses and people moved into the flood plain area to 

é near the center of activity. Largely through inertia, since some activities 

en are now in the flood plain are there uneconomically, many of these ac- 

| have remained in the area. 

_ Regulation of flood plain use has not been extensive because of: 

: 

. 


1) Lack of sufficient flood data. 

2) Lack of land use studies in flood plains. 

3) Objections of flood plain property owners and users to any measures 
limiting land use. 

4) Lack of legal precedent to prevent or restrict land use in the flood 


| plain. 


' Hydrologic and Hydraulic Aspects 


The engineer, administrator and planner, in considering flood plain zoning, 
3 interested principally in the following hydrologic and hydraulic information 
nd data: 


1) Flood seasons. 

2) Duration of flooding. 

3) Flood stages. 

4) Flood velocities. 

5) Extent of flooding. 

6) Minimum and/or economical channel widths. 
7) Frequency and magnitude of the greatest flood. 
8) Potential greatest flood. _ 

‘Knowledge of the flood seasons and probable duration of flooding will 
overn to a large extent what uses may be permitted on a flood plain. Flood 
tages, flood velocities, and the extent of flooding will assist in delineating 
1e flood plain and in determining the ground floor elevations and the type of 


‘: 


é 


: 
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buildings to be permitted in the flood plain. The extent of flooding and flood 
stages are generally required for the zoning plan. The last three items are 
required to arrive at an economical appraisal of the value of zoning a flood 
lain. 
; Nearly all floods are caused by excessive rainfall, which may be augment 
ed by water from melting snow in some areas. Heavy rains result from a 
combination of meteorological factors which may occur at any time in an 
area. Meteorologic and climatic factors vary from area to area and deter- 
mine the type of storms to be expected, the time of occurrence of such 
storms, and the probable severity of each type of storm in each area. The 
two principal types of storms which produce enough rainfall to cause floods 
are the thunderstorm and the stagnant-cold-front storm. Combinations or 
variations of these two types of storms, in addition to simple lifting of warm 
moist air masses over mountains, may result in flood-producing rains. 

The physical characteristics of drainage basins vary considerably from 
one basin to another. The upper portions of many drainage basins have rela 
tively steep gradients. In these areas floods will have a rapidly rising and 
falling flood hydrograph. The flood plains are narrow with little storage and 
velocities are high. Large quantities of debris and sediment may be trans- 
ported for considerable distances. As the streams approach the sea, the 
gradients decrease and the flood hydrograph has a long period of rise, a 
broad peak, and a longer period of fall. Considerable flood plain storage is 
usually present; inundation of the flood plains may be protracted and sedi- 
ment and debris transported from the upper reaches may be deposited on the 
flood plains. The rivers tend to meander and frequently change their course 

It is impossible to separate the flood plain from the river itself in order 
to consider its hydrologic and hydraulic aspects. All streams have flood 
plains along their entire length, although their width may vary from zero, in 
reaches which have perpendicular side slopes, to many miles in nearly flat 
plains, The behavior of the river in flood stages is determined not only by 
the physical properties of the normal channel, but also by the physical 
properties of the flood plain. 

Protection against the maximum flood which has occurred or is likely to 
occur is not considered feasible or desirable in some areas, particularly in 
small drainage basins.(4) Instead, protection is provided against a flood of 
given frequency. The question of whether or not partial protection against 
floods is justifiable is beyond the scope of this paper. The fact remains, 
however, that some flood control projects are being designed and constructe 
for partial protection only. Therefore, a very important consideration in an 
analysis-of a flood situation is the frequency with which floods of certain 
magnitudes are likely to occur. The frequency of floods is usually estimate 
by statistical analysis of the available data. The reliability of the analysis 
depends on the length of flood record, with respect to time, the accuracy of 
the data used, and the extent to which the flood frequency curve is extrapo-— 
lated. The frequency and magnitude of floods, together with an analysis of 
the damages experienced during each flood and the estimated cost of protec: 
tion for each flood, determine to a large degree the flood control design. 


Economic Considerations 


Until the advent of the motor age, it was economically advantageous for 
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y industries to be located in the flood plain.(5) Convenient water trans- 
psing usable processing water, a ready and inexpensive method of dis- 
osing of wastes, and the low cost of land more than compensated for the 
mages resulting from occasional floods. Some of these reasons no longer 
ce an economic advantage. Today, many industries are seeking locations, 


m or near major highways, which have a dependable water supply and attrac- 
ive surroundings. The problem of getting to work is no longer a major con- 
sideration since, with the widespread use of the automobile and good roads, 
employees are not only able and willing to commute to work but prefer the 
ider choice of home sites this mobility gives them. However, some indus- 
ries and utilities do have a valid reason for remaining in the flood plain. 

Flood plains have generally been considered as the ideal location for sew- 

ge treatment plants and wells for water supply. These facilities are neces- 
ary for the continued existence of modern cities. Those industries which 
pend on water transportation for their supply of raw materials, andasa 

means to ship the finished product, regard the flood plain as the only feasible 
ocation for their operations. In these instances, structural adjustments and 
mergency measures would alleviate the extent of flood damages. On the 
ther hand, telephone exchanges, hospitals, penitentiaries and essential in- 
ustries definitely should not be placed in the path of any conceivable flood 

ters. 

_ The rich alluvial soil in some flood plains is considered by some persons 
0 mean that the flood plain was ordained to have an agricultural role. In 
ther instances, the swampy conditions of the flood plain has led to the filling 

in or draining of these swamps to reclaim the land for industrial and residen- 
jal use, a role which is encouraged by the growth of cities and increasing 
and values. 

__ In some areas, the extensive use of the flood plain makes it more economi- 
al to provide protection to that area than to evacuate it. The protection, 

however, is only against a flood of a certain magnitude. There is no assur- 

ce that a greater flood will not occur causing considerable damage because 
he feeling of security generated by the flood control works reduces the pre- 
cautionary measures which were taken before the flood control works were 
installed. This feeling of security is likely to occur regardless of whether or 
not partial or full protection against floods is provided. The existence of the 
protective works, in addition, usually encourages further use of the flood 
plain. 

For most towns and cities, which have local flood problems, the economic 
value of the flood plain does not warrant extensive and expensive control 
structures. Frequently, it is more economical to develop a stream valley in- 
-0 a park rather than to enclose the stream in a limited-capacity culvert. (6) 
The flood plains adjacent to the stream, utilized for parks, playground areas, 
zolf courses and other recreational activities either public or private enhance 
the beauty of the city and contribute materially to the welfare of its citizens. 
The damage experienced during flood periods is minor and usually can be 
aken care of with the regular maintenance force. However, park buildings 
und facilities subject to extensive damage should be placed in areas which 
will not be flooded. On most small streams, the flood plain is inundated less 
han one per cent of the time per year. It is doubtful that this minor inunda- 
ion would affect the lawn areas to any extent if no oil or harmful chemicals 
vere present in the water. 


In determining the feasibility of flood plain zoning from an economic point 
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of view, the following items should be considered: 


1) Existing and potential flood damages. 

2) Cost of protection. 

3) Benefits from protection. 

4) Cost of evacuation. 

5) Benefits from evacuation. 

6) Allowable land use. 

7) Community benefits. 

8) Structural adjustments of existing encroachments. . 

9) Emergency measures such as flood warning services, rescue squad 
and elevating goods and machinery above possible flood heights. 


Social Implications 


In ancient times, the most favored locations for cities were those where 
conditions were suitable for agriculture and for transportation. The rich 
alluvial soil in the flood plains with the convenience of nearby water trans- 
portation satisfied these requirements. These cities have expanded, in most 
instances, to higher surrounding areas above possible flood heights. How- 
ever, the flood plain location of the original portions of many cities has beer 
and still is reflected in many of the social problems facing cities today. 

Periodic flooding of residential areas in flood plains fosters the develop- 
ment of blight. The unrealistic economic values, the high cost of providing 
services, the deterioration of the areas due to unwholesome topographic and 
physical surroundings—all contribute to the conditions existing in many de- 
veloped flood plains. The existence of some of these factors is evident in th 
difficulty in evacuating the area or in providing adequate protection to the 
area. This protection usually requires extensive and costly flood control 
works, the proportioned share of which is usually beyond the financial abilit 
of the people directly concerned. Therefore, the community, at large, is 
faced with three alternatives: 


1) Do nothing to the area and continue to pay for the high cost of 
services required and to absorb the cost of flood damages. 

2) Evacuate the area and retain it as public lands to be used for 
suitable functions. 

3) Provide adequate flood control works to protect the area. 


Protection of the Flood Plain 


The extent to which flood plains are zoned determines the degree of flooc 
damage prevention which will result. The minimum amount of regulation in 
flood plain prohibits land development which causes a construction of the 
natural stream channel. This is sometimes called the “channel-capacity” (3 
regulation since the degree of restriction depends on the capacity of the cha 
nel to discharge a flood of a certain frequency. The design frequency depen 
either on an economical analysis of the flood damages and benefits from pr« 
tection or on the degree of protection desired. Naturally, a flood greater th 
the design flood will inundate the adjacent flood plain. The maximum regul: 
tion of the flood plain limits the use of the flood plain to purposes which do 
not suffer from the action of a flood or to purposes which, if damaged, caus 


ASCE LUNETTA 1040-7 


ho serious inconvenience to persons or interests located outside the flood 
| plain. This regulation is often referred to as the “land-use”(3) regulation, 
which makes possible complete prevention of flood damage when this is con- 
| sidered necessary or desirable. 

Most of the existing flood plain ordinances in existence today are based on 
the “channel-capacity” regulation. The reasons for this are: (1) it is easier 
to determine, (2) its enforcement is simple, (3) legal precedents have been 
_ established, (4) the maximum amount of land is available to owners of flood 
plain property for use and (5) it directly affects only riparian owners. “Land- 
use” regulations are becoming more prevalent despite the restriction placed 
on the land. Often, however, complete restriction is considered unnecessary 
| and some use of the flood plain is permitted provided certain safeguards are 

met. 

Inherent in the establishment of any area as unsafe because of possible 
flooding is the implication that the areas which are not thus delineated are 
| secure from possible damage from floods. The degree to which the danger of 
this implication exists depends on the criteria used in establishing the limits 
of flooding. Flooded areas delineated by past floods may or may not include 
the area which could be potentially flooded by a future flood. It is very pos- 
sible that, as drainage basin characteristics and channel conditions are al- 
tered, a greater flood height than had been previously experienced may be 
expected. In fact, it is always possible to have a greater flood than any pre- 
vious flood regardless of whether or not there is any topographic change 
within the basin. 
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Discussion of 
“SYNTHESIS OF RAINFALL-INTENSITY-FREQUENCY REGIMES” 


by D. M. Hershfield, L. L. Weiss, and W. T. Wilson 
(Proc. Paper 744) 


D. M. HERSHFIELD,! LoL: WEISS, 2 AND W. T. WILSON,? A.M. ASCE.— 
We appreciate Mr. Sammons’ complimentary reference to multi-graphical 
correlation and his suggested improvements. In the evolution of figure 1 the 


_ primary requirement was to estimate 2-year, 1-hour rainfall for localities 
| lacking even daily rainfall data. The daily rainfall parameter was added to 
_ the rest of the diagram for application to localities having daily data. Before 


applying the diagram the writers tested it with independent data from many 


_ parts of the world, and noted that it had rather broad applicability with toler- 


able error. 

In arctic and subarctic regions this diagram was not as successful as an- 
other, shown here in figure 3. The large proportion of snow in the mean an- 
nual precipitation of arctic regions does not contribute much to a rainfall in- 
dex, and days-of-precipitation is a difficult observation to make in snow 


| regions. The two parameters, mean annual precipitation and number of days 


of precipitation per year, combine to give a measure of rain per rainy day. 


| Ina large measure the daily precipitation parameter which the writers added 
to the diagram in figure 1 is redundant to these two parameters. 


it has been observed that over much of the world the average ratio of 2- 
year, 1-hour rainfall to 2-year, 24-hour rainfall is about 0.40. This ratio 
may be 0.20 or lower in maritime regions having a low incidence of thunder- 
storms, while in continental regions where there is a high incidence of con- 
vective activity the ratio may exceed 0.60. Figure 3 shows the relationship 


| between annual extremes of hourly and daily rainfall as a function of thunder- 
| storm incidence. The diagram is based on data from stations in Alaska, 


Canada and northern United States. 
To apply this diagram to localities having more than 25 thunderstorms per 


_ year, the diagram would have to be extended. The shape of the diagram may 


vary slightly from region to region. We believe that this type of diagram, 
where the required data are available, and derived for the pertinent climatic 
region will be at least as good, simpler, and possibly superior to the one the 
writers published in the paper. 

While the foregoing discussion is not specifically a reply to Sammons’ 
discussion, the writers’ paper was presented as a progress report, and it is 


felt that recent progress is pertinent. 


Coming, now, to some of Sammons’ points: The writers have found that 
1.13 is a better average factor for converting annual series to partial- 


‘1. Weather Bureau, U. S. Dept. of Commerce, Washington, D. C. 


2. Weather Bureau, U. S. Dept. of Commerce, Washington, D. C. 
3. Meteorologist, Weather Bureau, U. S. Dept. of Commerce, Washington, 


iD: C. 


a 


. 
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duration series 2-year, 1-hour values; the 1.09 factor was based on a rela- 
tively small sample of data and, of course, even at the 2-year level there is 
some subjectivity in most methods of fitting curves to partial-duration data. 

In deriving the diagram that Sammons discussed, the writers were not 
concerned with the distribution of any of the parameters; they were dealing 
with a statistic, essentially the mean. The relationship of the 2-year value 
to the mean varies trivially among the several methods for analyzing annual 
series. In figure 1, no particular distributions were explicitly assumed. The 
scales used were largely a matter of convenience, with their inter-relation- 
ships fixed, of course, by the relationships of the data themselves, The em- 
pirical relationships of the figure are not dependent on subsequent discoveries 
as to the nature of the distributions of the variables used. 

It is not clear how Sammons obtained “Weather Bureau ratios” for the last 
table of his discussion. In the paper the writers explicitly suggested 5-year 
and 10-year ratios to the 2-year value of 1.4 and 1.7 as provisional average 
values in lieu of data for computing them. Weather Bureau Technical Paper 
No. 251 shows these ratios for about 200 United States stations. There is 
considerable variation among stations, but it is uncertain how much of this 
variation can be ascribed to physical and climatic causes, and how much is 
sampling error. 

Sammons raised a question about the writers’ stated record length averag- 
ing 40 years with respect to a base period. To avoid secular trends it is a 
common practice to restrict the data to a common period. This means elimi- 
nating data from outside this period, and stations whose records do not occupy 
all of the period. This period is sometimes made rather short to avoid ex- 
cluding too many stations. The writers felt that most of their records were 
long enough and for nearly enough the same period of time to be included. 
They felt that they would lose more by restricting the volume of data than by 
trying to get a homogeneous period. At present the writers are testing this 
issue, and at the same time are experimenting with the effect of duration on 
the frequency distribution. The writers might possibly find that the frequency 
distribution of daily amounts represents a different population than of hourly 
amounts, in some parts of the country. 

The writers might expect to find for some localities that hurricane rainfall 
defines the regime of daily extremes, while thunderstorm rainfall defines the 
hourly regime. In other localities they may find that the same storms fre- 
quently provide both the hourly and daily extremes. The writers want to find 
out under what conditions the relationship of daily to hourly extremes is a 
climatic index, and under what conditions it is a physical relationship. 


1. U. S. Weather Bureau, Department of Commerce, Weather Bureau Techni- 
cal Paper No. 25, Rainfall Intensity-Duration-Frequency Curves for Select- 
ed Stations in the United States, Alaska, Hawaiian Islands and Puerto Rico 
Washington, December 1955. : 
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Discussion of 
“GROUND WATER PHENOMENA RELATED TO BASIN RECHARGE” 


by Paul Baumann 
(Proc. Paper 806) 


PAUL BAUMANN, ! M. ASCE.—Mr. Lee correctly evaluates the writer’s 
equation 1, although expressing the acceptance rate in c.f.s. per square foot 
instead of per wetted acre would have been preferable for the following rea- 
sons: 

In the study of two-dimensional ground water movements a strip one foot 
in width is normally considered. Hence, the spreading grounds and the 
stream channels, the width of which rarely is constant, may be divided up into 
as many such strips as are required to reasonably reflect the effect of the 
variable width on the total acceptance rate. 

In constructing the flow net Mr. Lee apparently was unduly influenced by 
the classical textbook example of percolation from a reservoir per unit length 
of a floating dam of infinite length on a saturated, porous foundation, with the 
less featured case of flow from a spreading ground under a river levee. The 
two cases are fundamentally different. The get-away of percolation from a 
saturated stream channel is approximately at right angle thereto, namely, in 
the direction of residual pore storage and of a control, as shown in the writ- 
er’s Figure 3-a. In the longitudinal direction of the saturated stream channel 

: there is no residual pore storage and no control and therefore no get-away for 
additional sub-surface flow. Hence, in the absence of lateral escape, all per- 

Colation under a floating dam from a reservoir must rise to the surface, 

whereas, percolation under a levee from an adjacent spreading ground must 

hot so long as there is a get-away, by means of residual pore storage and a 
control. However, return flow to the river from the spreading ground may 
take place, depending upon the gradient between the water surfaces of spread- 
ing basin and stream channel as stated by the writer. 

The flow net shown by Mr. Lee does not consistently satisfy the basic re- 
quirements therefor. These are: (1) the vertical drop Ah between the 
€quipotentials along any phreatic line must be constant; (2) the flow 4q 
through each of the flow paths must be constant; (3) the flow 4q must be 
equal for all of the flow paths; (4) stream lines and equipotential lines must 
form quasi squares between themselves except near the boundaries where 
distortion is unavoidable. Not one of these conditions is reasonably satisfied 
by Mr. Lee’s flow net. It therefore is unbalanced and cannot be accepted. 

_ Fig. 11 shows a balanced flow net for the rather extreme condition similar 
to that chosen by Mr. Lee. The gradient between basin and stream channel 
water surfaces used in Fig. 11 is 2:1, as compared to the 5:1 gradient re- 
ferred to by the writer as the limiting value for insignificant return flow to 

pe stream channel. By assuming a number of phreatic lines through the levee it was 


rao Asst. Chf. Engr., Los Angeles County Flood Control Dist., Los Angeles, 
Calif. 
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possible to find the one for which the return flow is a maximum. For this : 
condition the return flow was found to be 7.5% of the total flow toward the : 
stream channel of percolation from the spreading basin. Hence, for a gradi- 
ent of 5:1 the return flow would indeed be insignificant as stated by the writer. 
In arriving at the return flow to the stream channel, Mr. Lee computes the 
flow through two vertical sections on the premise that the respective gradi- 
ents are equal to the difference in water surface elevation divided by the 
horizontal distance to the section. This is assuming horizontal flow through- 
out the entire reach in accordance with the Dupuit theorem. Under “Introduc- 
tion” the writer pointed out that... “the Dupuit theorem is applicable except 
in special cases” ... Subsequently, under the heading “Strip Basin” he stat- 
ed, “It is apparent that the Dupuit theorem is applicable with fair approxima- 
tion except near the spreading basin and near the control.” There the flow 
conditions cannot be ascertained without a flow net. The total flow computed 
from the flow net checks that in accordance with the writer’s eq. (1) within 
15%, which is a close check considering the extremity of the example. 

For the water year 1940-41 (October 1 to September 30) the rainfall for 
Los Angeles County was in excess of 200% of normal. Sustained runoff from 
the San Gabriel Mountains took place throughout the winter months and re- 
charge of ground water basins through percolation from stream channels and 
spreading grounds caused mounds to rise to the surface in the “forebay” of 
the Central Coastal Basin by the Spring of 1941. Fig. 12 shows the ground 
water contours and ground surface contours in the vicinity of the San Gabriel 
and the Rio Hondo Spreading Grounds of the Los Angeles County Flood Con- 
trol District as of May 6, 1941. These spreading grounds are located in the 
forebay. At stream channels the ground surface contours show average ele- 
vations between banks. Low water channel beds would be two to three feet 
lower. Of special interest is the merging of the mound beneath the San 
Gabriel Spreading Grounds with the channel mound, particularly in view of 
the fact that the get-away of percolation the east of the San Gabriel stream 
channel was more restricted than the get-away to the west of the Rio Hondo 
channel. Throughout this period there was no evidence of rising water in the 
saturated San Gabriel stream channel adjacent to the levee separating it from 
the spreading grounds, although the gradient from the spreading mound to- 
ward the stream channel is conspicuous in Fig. 12. Certainly there is no 
evidence of blocking by the channel mound of percolation from the spreading 
grounds. Of further interest is the curvature of the streamline trajectories 
normal to the ground water contours as referred to by the writer in connec- 
tion with the channel flow net. A typical trajectory is indicated by a dash line 
in Fig. 12. It is hoped that this will also serve to clarify the question in re- 
gard to curvature of flow nets raised by Prof. Brooks. 

Mr. Lee questions the similitude between the model test results and those 
of a prototype expressed by the writer’s eq. (1). The model test which was 
relatively free from capillary flow!3 showed a maximum discharge parame- 
ter q/K = 0.097 whereas eq. (1) for W = 2.4" (not 1.2" as used by Mr. Lee) and 
D = 12" gives q/K = 0.092, a difference of 5-1/2%. | 

Prof. Brooks undertakes to correlate the evaluation of the writer’s eq. (1) 
and his eq. (2) for the basic parabola. This is in keeping with fundamentally 
sound reasoning. However, in developing his eq. (2) so as to appear in terms 
of D/W Prof. Brooks apparently overlooked the fact that the depth D 


13. Ibid, p. 1059. 
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in his Fig. A can never be equal to D in the writer’s Fig. 3a. As previous- 


. iy14 shown the phreatic line between the spreading ground and the control 


can be approximated by a double parabola. D in Prof. Brooks’ Fig. A could 


| therefore be considered as the depth below the point (P) of counter flexure, 
_ common to both, the lower and upper parabolas. Hence, for identification, D 


and L in Prof. Brooks’ Fig. A should be changed to d and 1. Prof. Brooks’ 
example for arbitrary values of D/W = 2 and L/D = 10 and the apparent dis- 
crepancy between the results based on the writer’s eq. (1) and on the basic 
parabola are therefore misleading because of the fact that for potential flow 
between the spreading ground and the control the ratio q/K must be common 
to both the spreading ground and the lower basic parabola. As explained by 
the writer, eq. (1) was developed from flow nets which necessarily terminat- 
ed in basic parabolas. Hence, all corresponding values are found as the in- 
tersection of the 1/d lines with the (1 - e-4/W) curve in the writer’s Fig. 13. 
It is interesting to note that the 1/d = 3/4 line becomes tangent to the (1 - 
_e-d/W) curve as q/K approaches zero. In developing e-4/W into a series it 
_ also follows that as 


d/ weeds: G/Kee a/W 


Actually, as may be seen from Fig. 13 the q/K ratio closely checks the 
| d/W ratio if the latter is less than 0.1. 

In designing a spreading ground it is necessary to ascertain the economic 
limitations. Ordinarily, the depth D = 2d to the impervious stratum is known. 
Hence, by means of the writer’s eq. (1) and Prof. Brooks’ revised eg. (2) cor- 


responding values of q/K may be determined as shown in Fig. 14. 

| For the case quoted by Mr. Lee, namely, D ~ 100 feet, it follows from 

Fig. 14 that for D/W = 0.10 and therefore W ~ 1000 feet, the infiltration 

| parameter per square foot q/K = 0.05 or q = 3.35 c.f.s. per wetted acre for a 

-Meinzer Unit = 1000 as used by Mr. Lee. 

: Since Figs. 13 and 14 are correlated, the absolute value of 1 (which is half 

the distance between the edge of the spreading ground or stream channel and 

the control) may be ascertained from Fig. 13. 

As previously stated, extraction of ground water through pumping is to be 

expected where recharge through spreading is performed, Coordination of 

the two operations should include the locating of water wells as closely as 
possible to the theoretical control. However, the curvature in plan of the 
flow paths, previously discussed and indicated in Fig. 12, mitigates the need 
for precise location of the control, based on two-dimensional flow. For, 
prototype ground water movement, contrary to that in models, is seldom 
strictly two-dimensional but can, within limits, be approximated as such. 

Pumping in the spreading area proper, although effective in regard to in- 
filtration, may not be desirable because of a short path of percolation and 
possible contamination especially in shallow ground water basins. 

The dash line near the division of flow or axis of symmetry in the writer’s 
Figs. 3a and 3b as questioned by Prof. Brooks is not intended to be a bound- 
ary. If it were, the equipotential lines would not extend across it. Hence, 
flow takes place along the vertical division line and the dash line merely 
‘meant to indicate the effective limit of the flow Aq in the flow path next to 
the boundary. Admittedly, the effective limit of flow could be modified to 


T4, bid, p. 1042-1045. 
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advantage as shown in Fig. 15, the symmetrical flow net for the free water 
surface of stream channel and spreading ground at the same potential. 

Prof. Knapp, in presenting his interesting and constructive discussion has, 
indeed, found a flaw in the crystal ball in that it prophesied ultimate condi- 
dions when, in fact, they were just “prior to ultimate.” Of particular signifi- 
cance is Prof. Knapp’s reference to interface slopes under the influence of 
pumping. As may be judged from the writer’s Fig. (4) the draw-down gradi- 
ents adjacent to the test site are of the order of 1/2 of 1%. Hence, in using 
the 40 to 1 ratio as correctly suggested by Prof. Knapp, the slope of the in- 
terface would be of the order of 20%. This phenomenon is somewhat startling 
at first glance. However, normally sea water spread along the bottom of an 
aquifer is “sucked” into wells at moderate rates as compared to fresh water 
because of the lowest well perforations being at such a distance above the 
bottom of the aquifer as to require the interface to assume its maximum 
slope in rising to the well. 

This clearly brings out the difference between the saline wave and sea 
water intrusion proper, questioned by Prof. Hotes and Dr. Harder, in that 
during the passing of the saline wave the lowest perforations of affected welis 
may be temporarily blocked, whereas no measure of protection will save a 
well from abandonment if sea water intrusion is permitted to continue un- 
checked. Also, those wells adversely affected by the saline wave might be 
temporarily used to extract sea water cut off from the ocean by the barrier. 
In more timely terms the dissipation of the cut off part of sea water intrusion 
may perhaps be compared with containing the cut off part of an invasion army 
as compared to containing the whole invasion army. 

It was the writer’s privilege to observe the ingeniously conceived model, 
shown in Prof. Hotes’ Fig. 1, at the Richmond Field Station of the University 
of California. As is clearly brought out in Prof. Hotes’ discussion, this 
model lent itself to simulate well-nigh all probable prototype conditions and 
to evaluate results quantitatively within reasonable limits of accuracy. As in 
all cases where ground water studies are concerned, the time scale would 
have to be used on the usual assumption of homogeneity of the aquifer as well 
as uniformity in magnitude. 

Prof, Hotes’ eq. (1) which was developed by Dr. Harder and which is an 
adaptation of an example of flow through an earth prism analyzed by Dr. 
Muskat, is based on the assumption of a constant potential or equipotential 
along a vertical line through the injection well; a constant potential of equal 
magnitude along the horizontal sole and a variable potential along the verti- — 
cal line through the seaward control. Hence, the fresh water flow passes 
through the control between the bottom of the clay cap and the sole of the 
aquifer with no flow along the latter. Exactly the same result as obtained by 
Dr. Harder can be obtained by applying the Dupuit theorem to flow between the 
clay cap and the interface. The shape of the interface is such as to have bal- — 
anced pressure between fresh water and sea water at all points. Its shape is 
expressed by 


z=M+\/ MorciQextic Lemtaeae (2) in which, as 


shown in the writer’s Fig. 16 . 
z = Ordinate of interface above sole (x - Axis) 


M = Magnitude of aquifer 

x = Abscissa relative to center of recharge well 
q = Rate of flow per unit width of aquifer 

K = Permeability coefficient 

S = Specific Gravity of sea water 


Rs DISCUSSION 1041-11 
: 
| 
; 1 = Specific Gravity of fresh water. 


For any stable condition of the interface, the rate of flow q must be constant 
and a maximum. Hence, because the radical in the writer’s eq. (2) must be 
i the two terms under the radical for qmax must be equal. Thus for 
L 

| 

: 

| 

) 


M’=2qL ---------- (3) 
K(s-1) 
and q=1 (s-1) K = 1 (s-1) M T --- (4) 
2 L i i 


which checks Prof. Hotes’ eq. (1). 

now follows that for the above condition to obtain, the product qL of the two 
ariables must remain constant. Hence, as q increases, L decreases propor- 

FB otely and vice versa as is evident from Prof. Hotes’ Fig. 3. From the 

vriter’s eq. (2) it would follow that for any stable condition of the interface 

nd therefore any corresponding product qL, the height z of the interface at 
@ control is equal to M. As this would leave no aperture for the discharge 
fresh water, it is evident that eq. (2), based on the Dupuit theorem, is an 

\pproximation of the true solution. To satisfy z = o for x = o the sign before 
@ radical must be negative. 

From the writer’s eq. (2) it follows that the slope of the interface 


fiz ate 2 De GET aw i Nee nie (5) 


dx K(s-1) V MZ - 2qx 


K(s-1) 


ind that for x = L the slope for any corresponding q and L values is 90°, 
Vhile this answer again reflects the afore-mentioned approximation, it ap- 
ears to be reasonable if we remember Prof. Knapp’s 40 to 1 ratio between 
he fresh water gradient and the slope of the interface. In accordance with 
his ratio, any fresh water gradient in excess of 2-1/4 degrees would call for 
‘slope of the interface at the control of 90° which, of course, cannot be ex- 
eeded. 

If the interface were treated as an inverted, basic parabola and if the 
oundary condition z = M - q/K at x = L were imposed to bring the Dupuit 
olution in accord with the Kozeny!5 solution, then the writer’s eq. (2) would 


ead to 


q/K =- _L V_ 1? + M2 ---- (6) 
6-1) (6-1)? 


5. Ibid, p. 1040. 
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and q =e KL- +) K2L40) eh ine eet (7) 
(s-1) (s-1) 2 


While no authenticity is claimed for Ha (6) and (7), it is encouraging to 
find that the boundary conditions L=o’.q=TandL > ™©.°. q > 0 are 
satisfied. Also that q is of the same order as that from eq. (4). 


Putting eq. (6) inthe formq/K= L ( ) 1/4 (s- 1) ca) 
* 


_ L 
and expanding the radical into a binomial results in 


Mesa} iy = isiacai She Bieta (8) 


q/K = M’(s-1) - 
PANE L? 


i 
8 


from which it follows that for any finite discharge face the rate of flow q 
must be smaller than in accordance with Prof. Hotes’ eq. (1). Hence, the 
latter equation should be corrected to read 


in which n > 1 and increases as the ratio M increases. 
L 

Obviously, if fresh water from the recharge well were discharged from a 
pressure level, h feet above M.S.L. at the recharge well into a fresh water 
lake with the water surface at M.S.L. and if the interface were fixed, yo at 
the control would be close to zero and the pressure gradient or piezometric 
line would be as indicated in Fig. 16. If now it were possible to suddenly co 
vert the back pressure due to fresh water to that due to sea water, the dis- 
charge would stop just as suddenly. But the pressure head h would now exe! 
itself and would cause fresh water pressure to exceed sea water pressure e 
cept at the interface where the two pressures would be balanced. Hence, fo! 
a liquid interface the blocking of flow at the base of the control is offset, at 
least in part, by the pressure head y, on top of the control, 

This moves into focus the danger of “blow-out” at the control which groy 
with q and diminishes with the increase in the distance L between the toe of 
the salt water wedge and the control. Furthermore, this danger is greatest 
during periods of extreme low tide, a feature particularly significant at loca 
tions subject to large tidal ranges. Fortunately, the extreme tidal range 
along the coast of California is moderate, : 

It is possible that offset in the model, shown in Prof. Hotes’ Fig. 1, be- 
tween the back of the sea water weir box and the control face gave rise to tt 
trapping of air split out of the fresh water and that the q measurements wes 
thereby affected to some extent. 

Of considerable interest would have been the determination of the equigi 
tential lines and therewith the flow distribution between the interface and th 


“clay cap” of the model. This is just one of many facets meriting additiona 
research. 
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| The discharge through the control (16) to the ocean of fresh water could be 

more readily comprehended if the interface, as a basic parabola, were ex- 

tended so as to terminate at the horizontal discharge face, q/2K in width and 
lush with the bottom of the clay cap. The discharge of fresh water would 
hen be vertically upward and would thereby satisfy the physical premise of 
uoyancy of fresh water relative to sea water. 

Of special interest are the 10 stages shown in Prof. Hotes’ Fig. 4 and par- 
icularly those pertaining to high injection rates. It must be assumed that the 
ondition as shown at the bottom of the right-hand pictures would represent 
lose to the ultimate condition of the sea water wedge cut off by the barrier. 

Dr. Harder computes the contamination of 240,000 acre feet of water toa 
oncentration of 250 parts per million of chloride ions (Cl-) because of dif- 
usion of that part of the sea water wedge which had intruded landward from 

Ce injection wells. If sea water in intruding a fresh water aquifer mixed as 
eadily with the fresh water through diffusion as Dr. Harder seems to as- 
sume, the correction of sea water intrusion would indeed be a next to hope- 
less case unless it were always possible to locate the injection wells land- 
ward from the extreme tip of a saline wedge as suggested by both Prof. Hotes 
pnd Dr. Harder. At Manhattan Beach this procedure was carefully considered 
put had to be abandoned because of three principal features. First, the land- 
ward pressure gradient of the aquifer would have made it necessary to create 
barrier above mean sea level of approximately twice the height of that at 
he location where the barrier was finally located; second, the right of way 
necessary for the barrier including pipe lines, pressure regulators, chlori- 
tors, etc., would have been far too expensive to be justified in the light of 
experimental period of operation; and third, facilities for the disposal of 
ter produced during development of the recharge wells were completely 
acking. However, the writer still maintains that the effect of the saline 
ave, under similar conditions as those at Manhattan Beach, is negligible as 
compared to that of the sea water intrusion proper. It is only necessary to 
isualize the irregularities normally existing along the bottom of an aquifer 
the form of mounds and depressions to realize the facilities for ready ab- 
sorption of the sea water cut off from the original wedge by the injection and 
jorced landward. The average depth of sea water so absorbed may readily 
oe such as to require relatively short landward distance of spread. 
| In referring to the writer’s Fig. 8, Prof. Hotes concludes that the saline 
wave is worse than indicated by the writer. Apparently, Prof. Hotes over- 
j00ked the significance of the writer’s Fig. 7 which shows more tellingly than 
Fig. 8 what happened when the cut off sea water wedge was forced landward 
y the barrier wells. For example, it shows that the salinity in well G-8, 
180 feet from the line of recharge wells dropped from about 5000 p.p.m. of 
‘C1l- concentration to approximately 1000 p,p.m. in 16 months of barrier 
peration. Since October of 1955 this particular well has been immersed in 
resh water of less than 250 p.p.m. of Cl-. Prior to commencement of the 
echarge test, sea water displaced fresh water throughout the entire depth of 
he aquifer at a rate of some 500 feet a year. Hence, in the three years and 
vo months since the test was started unadulterated sea water would have 
»vertaken well G-8 and would have advanced some 400 feet beyond it. The 
Writer’s “optimism,” therefore, does not seem to be entirely unjustified. 
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Dr. Harder and Prof. Hotes’ points are however quite well taken in that 
sea water that had intruded into a shallow aquifer for a considerable distanc 
landward and had completely displaced the fresh water for such a distance 
could not be checked by means of injection wells without seriously affecting 
close-by inland wells by the displacement of the landward part of the saline 
wedge. On the other hand, wells close to the barrier would, under conditions 
similar to those at Manhattan Beach, be close to the ocean as well and there 
fore subject to permanent and irreparable instead of temporary damage. 

The model tests conducted by the University of California at Richmond 
have certainly paid dividends far beyond the return that could reasonably 
have been expected. 

The writer is deeply appreciative of the spirited discussions through 
which great value was added to his paper. He only regrets his inability to 
agree with all of the thoughts therein presented. 
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Discussion! of 
“EXTENDING STREAMFLOW DATA” 


by W. B. Langbein and C. H. Hardison 
(Proc. Paper 826) 


F. A. JOHNSON, 2 A.M. ASCE.—A graphical method is presented for esti- 

_ mating the relationship between monthly discharges of nearby streams having 

similar characteristics. The relation curve then is to be used to extend the 

Shorter record on basis of the other, or long-term record; and statistical 

| measures are suggested as an indication of the expected accuracy of the 

| estimates. 

_ Several experiments made by the writer with pairs of long-term records 

i led to fairly reliable extensions in some cases, but not in others. The statis- 

|tics, such as standard error of the relation curve and the coefficient of cor- 

relation, did not invariably give a dependable indication of the accuracy. 

| It was apparent from these tests that in some regions the areal pattern of 

!runoff was relatively stable from year to year, and over long periods. In 
other regions, however, the distribution varied greatly from year to year, or 
between irregular periods of several years. These instabilities or changes 
in pattern of runoff were noticed in a number of regions, both of the Pacific 
and Atlantic Slope basins. They apparently were caused by variations in the 
areal distribution of precipitation. 

| An average relationship between the monthly discharges of two nearby 

streams, say for 5 years of concurrent records, thus may differ considerably 

{oom that during another period, and lead to systematic errors in the esti- 
mates. An obvious way to minimize such errors, as well as the random 

| 4 is to establish a relation between the discharge in the subject or 

{satellite basin and the average discharge fromtwo or more nearby basins on 

jdifferent sides of the one under consideration. For correlation of monthly 

idischarge, however, this possibility often is limited by lack of sufficient long- 

iterm records having similar characteristics as to monthly distribution. 

| One purpose of this discussion is to describe a method for estimating an- 

mual or long-term runoff of the satellite stream from two or more long-term 
base records, regardless of the monthly distribution of runoff. In many parts 

jof the Atlantic and Pacific Slope basins this leads to fairly accurate estimates, 

jwhich may be compared with extensions of monthly records on an annual ba- 
is. Such comparisons are, in effect, independent tests of the accuracy of 

ithe extensions, if the annual estimates are based on different index records; 

d they sometimes may be used for purposes of adjustment. 

The need for such a test is demonstrated by an example of the authors’ 
orrelative estimates which is summarized in Table 1 and 2, and Figure 4. 

{in a number of similar experiments it was found that substantial errors of 


. Publication authorized by the Director, U. S. Geological Survey. 
2, Regional Hydr. Engr., Geological Survey, U. S. Dept. of the meee 
ee aconia, Wash. 
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the kind illustrated here were not uncommon, although generally of somewhat 
lesser magnitude. 

In this example, the monthly runoff of the Wading River near Norton, Mass. 
was estimated from that of the Taunton River at State Farm, Mass. The 
drainage areas are 42.4 and 260 square miles, respectively. The center of 
the Wading River basin is about 15 miles east of that of the Taunton River 
basin. Both basins are near the Atlantic Coast in the southeastern part of 
Massachusetts. A relation curve was defined on logarithmic paper by 60 
pairs of observations recorded from October 1945 to September 1950. (Run- 
off figures in inches, as adjusted for diversions and storage, were readily 
available and were used without allowance for differences in length of months). 
The graphical method illustrated by the authors in Figure 1 was used for fit- 
ting the curve and estimating the statistical measures. The standard error 
of estimate was determined as 0.10 log unit, and the coefficient of correla- 
tion as 0.975. 

Estimates made for the period 1930 to 1945 from this curve are summa- 
rized by annual totals in Column 2 of Table 1 for comparison with the re- 
corded figures of Column 3. The runoff for two 4-year periods was under- 
estimated by about 20 percent, which is approximately the standard error as 
determined for monthly estimates. The relationship during the sample peri- 
od, 1946 to 1950, evidently was not representative of other periods between 
1930 and 1945. 

The authors pointed out that the statistical measures estimated from the 
sample are applicable in other periods only if the mean line is correctly lo- 
cated. However, it is difficult to see how this can be judged from internal 
evidence. It was pointed out also that these measures have the usual signifi- 
cance only when deviations from the mean line have a normal distribution. 
Although adjustments such as the logarithmic transformation often appear to 
bring the deviations into a normal distribution, the deviations may resemble 
random variations only if regarded without relation to time. The writer has — 
noticed that large positive or negative deviations of monthly values some- 
times persist for periods of from 4 to 6 months. These consistent deviations 
sometimes result in an accumulation of error in the annual totals almost as 
large as the percentage errors of the monthly estimates. Because of such 
errors it is evident that the computed standard error of the relationship may 
be an underestimate even for the sample period. It may be grossly in error 
during other periods if there is a bias of the kind illustrated in Table 1 and 
Figure 4. 

It was noticed also that the seasonal distribution of discharge of some 
streams of the Pacific Slope is such that the monthly values tend to cluster 
in zones above and below the mean and median. Under these circumstances 
the coefficient of correlation does not have the same significance as for a 
normal distribution, and is apt to be fictitiously high. Coefficients of more 
than 0.99 in fact did not appear to be uncommon. 

The recorded annual runoff for the example of Table 1 is shown in Figure 
4. As would be expected, the annual values of runoff computed from monthly _ 
estimates fall approximately on an average curve through the recorded an- 
nual figures for the years used in definition of the monthly relationship. The 
degree of correlation between annual discharges of long-term records ina _ 
given region thus furnishes an indication of the reliability to be expected of 
monthly estimates for other streams in that region. Such comparisons 


would seem to be a useful preliminary step in extensions of streamflow data 
such as proposed by the authors. 
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It has been the writer’s experience that the possibility of large systema- 
tic and random errors can be minimized most effectively by establishing a 
relationship between the subject stream and the average of several nearby 
streams, rather than between two streams alone. Along the Pacific Coast 
correlations seem to be the most reliable where the index basins are on op- 
posite sides of the subject basin, and on a line roughly normal to the west- 

t, or average direction of the storm paths. A similar condition may exist 

other regions. In connection with their discussion of snow-fed streams, 

e authors mentioned the possible need, at times, for use of more than one 
ex of discharge to improve the correlation. It would seem that such pre- 
aution should be taken as a matter of course, since it can seldom be known 

t a group of observations in one period is closely representative of condi- 

ions in another period. 

If the monthly distribution of runoff in the several basins is similar, the 
erage of several records may be used in the same way as was illustrated 
the authors with a single index record. However, where there are sub- 

stantial differences between the altitude characteristics of the basins, and 
tunoff results both from rainfall and snow melt, the monthly distribution may 
fer considerably. In that event, average indices of annual runoff may be 
ed for defining relationships of annual runoff. Estimates of annual runoff 
jhen will serve to verify the accuracy of annual sums of monthly figures es- 
jimated in the way proposed by the authors. 
| Although the error due to variations in the areal distribution may be mini- 
ized by using averages, a relation line defined by only a few pairs of annual 
alues obviously may not be dependable for purposes of extension. Fortunate- 
, in humid regions where the precipitation is substantially more than the 
ormal evapotranspiration, fairly reliable estimates usually can be made by 
. computed relationship described in the following paragraph. This has been 
ised by Knox and Nordenson(1) for studies in the northeastern part of the 
outs, but its wider applicability may not be generally recognized. 
Ignoring systematic and random variations in water loss, and changes in 
tural storage, the relationship between the annual precipitation on nearby 
reas can be expressed as follows: 


| : Ry + L = b(R, + L,)s 
| 
: 


iyhere Ry and Rx are values of annual runoff, Ly and Ly are estimated values 
lf mean annual loss, and b is the ratio of precipitation on basin “y” to that on 
isin “x.” Both terms of the equation are subject to error equal to the devia- 
ion of the actual losses from the assumed average losses, plus the changes 
n natural storage. However these errors tend to be the same for nearby 
reas, and thus tend to cancel out. The ratio, b, often can be estimated 
osely from the average of only a few pairs of annual observations, even if 
meir range of magnitude is small or negligible, so that a relation curve could 
be defined graphically. The relation, Ry = b(Rx + Lx) - Ly then can be - 
ed to estimate values of Ry from Rx during periods of missing record. 
he results usually are affected only slightly by fairly large differences in 
ie assumed values of average loss, Lx and Ly, from actual values, or by 
riations in the amount of annual loss with variations in precipitation. In- 
mation regarding natural water losses can be obtained from a number of 
urces. Some data were compiled by Rafter(2) as early as 1903, later by 
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Williams and others(3) in 1940, and by Lee(4) in 1942. 

This method of estimation is illustrated in Table 2, as a supplement to 
Table 1. (In this example the ratio, b, happens to be 1.00. The method 
seems to work just as well where there are considerable differences in the 
values of areal precipitation). The annual estimates of Table 1, as well as 
for a number of unlisted experiments, are in general substantially more ac- 
curate and consistent than totals derived from monthly correlations between 
two streams alone. The improvements were due to the use of more than one 
index record for the annual estimates. The same improvement in monthly 
estimates could be obtained by this means, where the monthly distribution is 
closely similar in all of the basins. 

In an actual extension of records, the reliability of the estimates can be 
judged only on basis of the errors for the years of definition, or by compari- 
son with estimates made from independent correlations with other groups of 
records. The standard error of annual estimates computed for only a few 
years of record obviously may be subject to considerable error, but it at 
least affords a crude measure of the probable magnitude of error if bias has 
been substantially minimized. Since relationships between annual figures 
can be established rapidly, this provides a convenient means for selection of 
groups of index records that probably will be the most reliable for extension 
of monthly records in the way proposed by the authors. 

This method of analysis was applied to the example given by the authors 
for estimation of the discharge of Murder Creek near Evergreen, Alabama 
from that of Sepulga River. As shown in Table 3, there was a better correla- 
tion between the annual runoff of Murder Creek and the average of two 
streams, Sepulga River and Escambia Creek, than between that of Murder 
Creek and Sepulga River alone. The Escambia Creek basin is on the other 
side of Murder Creek, about 30 miles to the southwest. (The indices for the 
Escambia Creek basin were multiplied by a constant factor, 0.88, to avoid 
giving greater weight to the variations in that basin than to the variations in 
the Sepulga River basin, because of the relative magnitude of the figures). 
Since the monthly distribution of the runoff in the three basins is similar, 
the comparison of the standard errors for the period of definition, 1941 to 
1945, suggests that a substantial improvement in the monthly correlation als« 
could be obtained by use of average indices. The relationship between the 
monthly discharge of Murder Creek and the average of the other two streams 
(in amounts per square mile) in fact is nearly a linear one for the entire 
range, and the standard error for the sample of 60 values is only 0.05 log 
unit, or about 60 percent of the authors’ value. This improvement is roughly 
the same as between the two sets of annual estimates listed in Table 3, for 
the years 1941 to 1945. 

Statistical tests, whether of monthly or annual values, will not show 
whether the relationship for the period of definition is biased in relation to 
that for another period, as in the example of Figure 4. For that reason it 
would seem desirable to compare the results of monthly estimates with an-— 

nual estimates based on the average of several index records whenever the 
monthly indices themselves cannot be averaged to minimize the possibility ° 
bias. 

If it appears that estimates of annual runoff for a given stream are apt to 
be more reliable than estimates of monthly runoff, the latter could be multi- 
plied each year by the ratios of the annual estimates to the sum of the month 
ly estimates. It would have been desirable to make such adjustments in the 
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| estimates of Table 1, for example. This would have brought the estimates of 
mean runoff for the 4-year periods, 1935-1938 and 1942-1945, from within 
roughly 20 percent to within a few percent of the correct values. If the peri- 
| ods happen to be critical ones from the standpoint of dependable water sup- 
plies or regulation by storage, such improvement would be of considerable 
| importance. 
: The graphical method proposed by the authors for estimating the relation- 
ship between monthly values of runoff is simple and can be applied rapidly. 
It is much more practical than numerical analyses, and considering the pos- 
sible magnitude of sampling error in short-term correlations, more elabo- 
rate methods of analysis usually would not seem appropriate. The writer 
agrees that extension of short-term records is of growing importance and 
that the suggested methods constitute a useful approach to a sound procedure 
for extending records of monthly discharge. It is believed, however, that the 
| statistical measures should be viewed with caution as indices to the reliabili- 
ity of such estimates and that there is a need for independent estimates and 


)controls of the kind proposed in this discussion. 
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Table 1 


Estimated runoff of Wading River near 
Norton, Mass., compared with recorded runoff. 


(i) (2) (3) (4) (5) (6) 
Runoff, 
Runoff, inches Error inches Error 
a in @ of b in %@ of 
Water Sum of monthly estimated Annual estimated 
Year estimates Recorded values estimates values 
1930 13.4 12.6 + 6.3 1365 + 6.7 
31 20.4 Corel -16.0 24.5 + 3.3 
32 15.5 15.1 + 3.1 L738 +15.2 
55 31.0 29.65 + 5.6 31.4 + 6.7 
34 2lee ome = 9.6 Coe ag 0.4 
1935 tec 25.00 -36.6 24.2 + 2.9 
36 20.6 PASI S -22.3 24.4 > Katey | 
37 24.1 24.8 - 3.2 26.0 + 4.6 
38 etel Shirts! -24.7 Sen - 0.3 
39 24.2 Pa, +1.5 26.2 + 8.8 
1940 20.7 22.0 - 6.3 22.4 +1.7 
41 15.4 14.8 + 3.9 LS eit + 5.7 
42 Jip 16.1 -39.9 15.3 - 5.2 
43 17.7 ale -19.8 21.5 +1.4 
44 Sel 11.8 -21.8 ASiet +929 
1945 22.0 eit -25.9 25.2 -. 9.9 
Mean 19.5 21.8 -11.8 22.4 + 2.7 


a/ Monthly estimates were based on the runoff of Taunton River. 


b/ Annual estimates were based on arithmetical averages for 
Taunton River at State Farm, Mass., Charles River at Waltham, 
Mass., and Blackstone River at Woonsocket, R. I. 

The relationship of the indices to the runoff of Wading 
River was established for the period, 1946-50, by the method 
illustrated in Table 2. The Charles basin (227 sq. mi.) is 
about 20 miles north of Wading River, and the Blackstone 
basin (416 sq. mi.) is about 25 miles northeast. 


Figures of recorded runoff were taken from Water-Supply Paper 
1301 and rounded to the nearest tenth in this table. 


+SCE 


n Ol. 
ater R, Ry R, R pilus times minus 
ear * 25.0°. 1,00 - 25.0 

1946 55 ck bee6<o o 2Oe0 1120.5 55.6. 53.6 =. 28.6 

47 «19.9 17.1 18.8 18.6 43.6 43.6 18.6 

eemeOesO6 C8al 24.9 . 27-1 52.1 52.1 27,1 

fee 2225 = 15.4 16.6 -18.2-4,43.1- 43.1. 18.1 
38.9 38.9 


Annual runoff (adjusted for diversions and storage ) 


Taunton R. at State Farm, Mass. 


R_: Wading R. nr. Norton, Mass. 
Ratio, b: 


DISCUSSION 


Annual runoff of Wading River as 
estimated from relationship to 
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Table 2 


indices of areal precipitation; 1946 to 1950 


All data in inches 


Charles R. at Waltham, Mass. 


Blackstone R. nr. Woonsocket, R.I. 


(Ry + Ry + Rz)/3 


231.1/251.5 = 1.00 


R Error 
dy 
29.6 -1.0 
16.7 +1.9 
29.6 -2.5 
18.6 -0.5 


Standard error of estimate: V Z ay* /(n - 1) = 2.02 inches 


Estimates made for the period 1930 to 1945, as in column 8, 


are listed in Table l. 


period was 1.40 inches. 


The standard error during this 


ae 


Columns 6, 7 and 8 are unnecessary in this example, 
but are shown to illustrate the usual procedure. 


Note: 
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Table 3 


Annual runoff of Murder Creek near Evergreen, Ala., 
estimated from indices of areal precipitation. 


Year 1 aby 2 Ry R dy : oy, 
1941 46.4 42.6 43.0 42.8 15.5 17.0 = PSs 15.6 =Let 
42 57.82 -51).0r 49:59) 6 50.42 oD all Ome ewe cool eles 
43 56.9 «50.12-449.3-49.7 122.922.640.631 922.9. = O5k 
44 6824 60.35 ~65.5., 62.9 936.8 35.4 “+154 39.9" 44.5 
1945 53.6 47230. 42.5. 44.8 9 27.6 619.9 © -2.S me lO el see 
46 66.2 58.4-- 68.0 63625" Siu eo Gera 10. 96425 Oia 
47 62:55) 5439.) 55 .4.254.2 2 2760 8 2064) O65 2 COs ees 
48 6055) *53:..4.5 5056 al in Coes a coe Ose mee eae 
49 69.9, 161.5. 62.0 mrOece a OOst | Ooet . toes Obes aetoalf 
1950 49.35" 245.0) 6 olet 40,4. 12.98 4 Peele 9.8 -4.3 
Mean 25.5 20.4 40.1 25.7 +065 


Average loss: 30 in./yr. (estimated) 
b= tere Pe oerl 062 
bt= (R+150)/ P, = 1.065, 
Where values of P and R are totals for period, 1941-45. 
Ry = 1.062 P - 30 
Ry' = 1.065 P, - 50 
Py = Annual runoff of Escambia Crk. + 30 inches 
' = 
PL 0.88P, 
PS = Annual runoff of Sepulga River + 30 inches 
E = Average of P,' and Po 
R = Annual runoff of Murder Creek 
dy and dy' = Errors of computed values, Ry and Ry'. 


Standard error of estimate 


Estimated from deviations d 
1.90 inches, 1941-45 
1.69 inches, 1946-50 

Estimated from deviations, dy! 
3.46 inches, 1941-45 


4.35 inches, 1946-50 
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LEGEND 


Circled points: Runoff for 
indicated years ending 
September 30. 

Solid circles: For period 
used in definition of 
monthly relation. 

Dashed line: The average 
of annual totals (shown 
by crosses) as derived 
from monthly estimates. 
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Discussion of 
“FLOOD PLAIN ASPECTS OF RIVER PLANNING” 


by Anthony M. Lunetta 
(Proc. Paper 1040) 


ROBERT L., SMITH,* A.M. ASCE.—The role of the flood plain in river 
planning is seldom limited to a specific type of water problem, nor is it al- 
ways limited to consideration of just the water resource. Under certain con- 
ditions problems of land, wildlife, or mineral resource management can be 
expected to influence materially the planning decisions governing water re- 


. source development. By contrast the remarks in the foregoing paper have 


been limited to a specific water problem—floods. Likewise the remarks on 
the possible ways to minimize the damages caused by those floods were 
limited to the proposal that the use of land within the flood plain be subject 
to regulation. Further, those remarks were summarized in a very logical 
sequence and they set forth quite adequately the factors to be considered. 
Notwithstanding those various limitations, the most significant feature de- 
veloped in those remarks was the acknowledgment of the breadth and, in 
essence, the complexity of the problem discussed. Of particular importance 
is the concept that this problem of wise utilization of the flood plain is inter- 
disciplinary in nature. In fact, engineering consideration of the problem can- 
not dictate the answer, but can only serve to crystallize the solution to 
fundamental economic, legal, and social questions. 

The subject presented is important. In every area of the country, in areas 


both developed and underdeveloped, day by day expansion is continually in- 


creasing the potential flood damage. The author suggests that this expansion 
is twofold. First, there is encroachment which is motivated by economic de- 
sire to utilize the resources of the flood plain. Secondly there is encroach- 
“ment which did not seek the flood plain for a specific economic purpose, but 
which, because of ignorance of the risk involved, blissfully located in the 
green valley. It must be acknowledged, however, that all investments under 
our free enterprise system acknowledge and allow successfully for various 
risks, and it is submitted that in the final analysis either inability or refusal 
to evaluate the risk is responsible for all unwise encroachment. Thus, there 
is presented the need for educating John Q. Public on the virtues of determin- 
ing his calculated risk. 

Engineers and laymen confronted with flood plain losses can agree on a 
common objective, namely—“The reduction of flood damages by more intelli- 
gent use of the flood plains.” In the paper just presented it is implied that 
one way of attaining this objective is through the utilization of land use regu- 
lations. It remains to be seen whether engineers, the courts, and the public 
will be willing to adopt a regulatory type approach. Four reasons were pre- 
sented why widescale adoption of regulatory procedures has not been realized. 
Those reasons were briefly (1) lack of hydraulic data, (2) lack of land use 


¥ Director, Iowa Natural Resources Council, Des Moines, Iowa. 
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data, (3) objections by property owners, and (4) lack of legal precedent. It is 
submitted that these four points again illustrate the present inability of indi- 
viduals, courts, developers, and others to evaluate properly the risk involved 
in certain flood plain development. 

A secondary primary problem to be considered is tightly interwoven 
among the various disciplines. Simply stated the problem becomes one of 
defining what is the right of the public and what is the right of the individual. 
Primary arguments for public control hinge on protection of the general wel- 
fare from the direct adverse physical effect of unwise development, and the 
secondary effect of modern floods on the community and the regional economy 
and public health. Countering these public needs is the desire for preserving 
greatest possible individual freedoms. The only adverse effects sustained 
from many encroachments is the damage to the facility itself. If that same 
facility can be shown to represent a long term economic gain can it be clas- 
sified as a non-conforming use? Further, floods come in many shapes and 
sizes—there is the extraordinary and the routine. Does the same public 
right hold for all possible areas of inundation? The answer to these and 
other questions go far beyond the scope of this discussion. They will not be 
answered succinctly or quickly. Their answers will, at best, evolve with 
time for they are simply a portion of a more diverse question of local, state, 
and national policy. A question which the Engineers Joint Council recognized 
when they stated “The application of desirable water policies can only be car- 
ried out ... . when an overall policy of land utilization is available.” (1) 
Paradoxically, development of such a policy for flood plain areas must await 
a reasonable attainment in public understanding of the problem. 

The foregoing comments implying the need for an educational approach to 
the problem of flood plain use have not been asserted in opposition to the 
author’s implied plea for positive zoning procedures. However, the educa- 
tional approach has been grossly overlooked at the local, state, and national 
level, Presumably the responsibility for development of adequate education- 
al program could be a proper function of government. 

The writer once suggested the desirability of interpreting and presenting 
available hydraulic data in the form of river plan and profile sheets.(2) Ade- 
quate strip maps can be developed in many parts of the country from avail- 
able aerial photos. Sufficient information such as river mileages, section 
centers, and major highways would have to be identified to allow correlation 
of plan and profile. The profile sheets might well include a low water profile. 
profiles for floods of various percent chance of occurrence, and a profile for 
a so-called provisional or project design flood. Seldom could a stream be 
found where the available data seemed sufficient in every detail. However, 
treated as a design problem the task would reduce to certain standardizing 
procedures. For example, recent utilization of regional analyses of stream | 
flow data could be a useful tool for such work.(3) It seems apparent that such 
data would be in much demand once investment houses and others realized the 
significance of the message it contains. Particularly would this be true | 
where those who seek flood plain development are required to pay their pro- 
portionate share, based on benefits received, of the cost of future protection 
works. 

Finally one cannot lose sight of the importance of the author’s last state- 
ment. Bigger floods are yet to come. Zoning or less formal approaches to 
land use regulation can only be effective up to a certain limit. Thus land use 
regulations cannot be advocated as a panacea or complete cure, they simply 
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represent another tool for minimizing flood damage. If the thoughts ex- 
pressed in the paper are sufficient to reawaken in the minds of only a few the 
possibility of utilizing that tool, the paper will have served its purpose. 
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COAGULATION AND SEDIMENTATION@ 
J. M. DallaValle> 
(Proc. Paper 1052) 
SUMMARY 


By assuming initially a monodisperse system of particles or nuclei, the 
chance of collisions leading to 1, 2, 3, etc. aggregated particles is computed 


_ according to the procedure used by Smoluchowski. The effect of settling is 
_ ignored by assuming the process is very fast. A discussion is included of the 


general law of settling of dilute suspensions as well as highly concentrated 
suspensions. 


The breakup of matter into finer and finer particles can proceed only to a 
certain point before the phenomena of coagulation or coalesence begin to 
bring them together again. It may well be that the forces tending to disinte- 
grate matter are eventually counterbalanced by others which act to bring par- 
ticles together. In industry, fine particles are made to join together by a 
form of coagulation called polymerization and this is the basis of plastics 


_ manufacture, rubber synthesis and many other products. Often it is a transi- 


_.tory phenomena on the way to making a final product and thus receives little 


attention. Or, as in the case of water purification, it is a carefully controlled 
phenomena. Coagulation is not a phenomena characteristic of liquids alone. 
The manufacture of metallic oxides by burning pure metals exhibits this 
phenomena. In fact, with many metallic oxides, particle size control is 
achieved by letting the particles “grow”—in other words, to aggregate or 
coagulate while still in suspension. ** 

The tendency of particles to coalesce or coagulate is far more widespread 
pan supposed. There can be no question but they act in all fine particle 


Note: Discussion open until January 1, 1957. Paper 1052 is part of the copyrighted 


- Journal of the Hydraulics Division of the American Society of Civil Engineers, 
Vol. 82, No. HY4, January, 1956. 


a. Presented at the Dallas, Texas, Meeting of the ASCE, Feb. 13-17, 1956. 


‘b. Georgia Inst. of Technology, School of Chemical Eng., Atlanta, Ga. 


** We shall use the terms coagulation, coalesence, aggregation, cohesion and 
accretion synonymously and understand that they apply to the growth of 
__ particles by accumulation to their surfaces of other particles of the same 


species. 
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suspensions and are essential to the formation of stable sedimentary de- 
posits. The constant stirring up of tidal basins, channels and shorelines by 
wave motion, of seasonal upheavals of lake bottoms by thermal inversions 
and even atmospheric dusts as they are swept about by winds would make 
this a-difficult world in which to live in the absence of coagulation. 

This paper proposes to outline the mechanism of coagulation and the bind- 
ing forces between particles once they stick together, and then to show how 
particles settle to form the tremendous sedimentary deposits which consti- 
tute so large a part of the earth’s crust. 


Smoluchowski’s Theory of Coagulation 


It is noteworthy that since the fundamental investigations of Hardy (1900), 
Wiegner (1911), Smoluchowski (1918) and more recently, Muller (1926, 1928), 
Tuorila (1927) and Fuchs (1934), the basic theory of coagulation has been only 
slightly modified. It is especially to Smoluchowski that we owe our present 
appraisal of the mechanism of coagulation. Muller and Fuchs have extended 
certain aspects of Smoluchowski’s theory, but for all its simplifications rec- 
tified in many respects by later investigators, Smoluchowski’s basic concepts 
remain unchanged. His theory is applicable to all fluids—liquids and gases— 
and his equations have more universal applications with regard to coagulation 
phenomena than any of those derived with greater attention to detail. Consider- 
ing the many extraneous complications affecting coagulations suchas (in liquids 
the electrolyte atmosphere and accompanying double-layer effects, and (in 
gases) the existance of electrostatic charges, it is indeed remarkable that 
Smoluchowski was able to develop so successful a theory; for these factors 
affect the mathematical treatment so as to render a solution of the problem 
impossible. 

Smoluchowski used a very simple model for deriving his coagulation equa- 
tions. He assumed that the particles were initially uniform in size and ran- 
domly distributed. He regarded them as subject to Brownian motion during 
the whole period of coagulation. He assumed that there were collisions be- 
tween the particles as a result of their random motions and that on collision 
the particles formed a single particle. Rapid coagulation was regarded as a 
process in which particles on every contact adhered to each other on every 
contact; slow coagulation was attributed to the fact that only a fraction of the 
particles adhered on collision. We shall be concerned only with the first of — 
these two types of coagulation in this paper. Why particles should adhere to 
each other on first contact in some instances and not in others will be dis- 
cussed briefly later. 

As indicated or implied above, coagulation is that process which leads to 
the growth of particles in a fluid atmosphere by accretion of other particles. — 
If initially there are no particles per unit volume and no settling occurs, at 
atime t after coagulation begins, there will be no change in the mass of 
particles suspended, but there will be some kind of a distribution consisting | 
of single particles, double particles, triple particles, etc. If the particles — 
remain small even after the coagulation process has proceeded for some 
time, or if the process is very rapid, we may ignore the effect of gravity. 
However, with time, all particles come under the influence of the gravitation- 
al field and the suspension will no longer consist of randomly distributed 
particles. In general, if the particle diameter is a function of time, D(t), the 
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density of the suspension at any point z (from the bottom of the container) 
and time t is, after coagulation has ceased 


| De 
B(z,t) = po +6, [(e, - p,)/5] j F(D,) d(Dg) (1) 


° 


where Po and ps are the densities of the fluid and particle, respectively, 
Co is the concentration of the suspension at the moment coagulation ceases, 
De is the effective Stokes diameter and F(De) its size distribution. In using 
Eq. (1), we may take t = 0 at the instant the rapid coagulation process is 
completed and the particles begin to settle. Both ¢ andC can, of course, be 
readily expressed on a number basis when it is convenient to do so. 

We shall let n denote the number of particles in a unit volume at time t 
after coagulation has started. The reciprocal of this number, 1/n, can be re- 
garded as the average space occupied by a particle. If we have some method 
| of measuring the change of n with time, we shall find that a plot of 1/n vs 
t yields a straight line with a positive slope. Experimental methods are 
available for determining n from instant to instant in a coagulating system 
and for both sols and aerosols, provided gravity has not affected the random- 
ness of the distribution, data have been obtained which yield the simple but 
important equation 


(1fn) - (I1fnj) = xt (2) 


in which K is called the “coagulation constant.” Since the n’s are usually 

| quite large, it follows that the constant K for a given system is very small. 
| The constant K will, of course, depend on temperature and whether the sys- 
tem is quiescent or turbulent. 

Actually, Smoluchowski derived Eq. (2) using the following considerations: 
| First he regarded the system to be initially monodisperse and randomly dis- 
| tributed and consisting of an average of ng particles per unit volume. The 
| volume of all the solids together was assumed negligible in comparison with 
| the volume of the liquid. Next, with time increasing and with Brownian mo- 
tion, all or a stated fraction of particles on collision adhere to each other. In 
order to calculate the number of particles colliding with a particle (assumed 
| fixed), we may consider the particle to be surrounded by an imaginary sphere 
) of radius ro and apply the Fick diffusion equation giving the number of par- 
| ticles n' diffusing through a surface. This is a procedure for a steady state, 
| but it can be shown to be sufficiently accurate. The Fick equation is 


nis= Peal (dn/or ) 


where D is the diffusion coefficient. We can integrate this equation subject 
to the conditions: 


~~ 


= for r = 00 


n= Osfor r= 2, 
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The last condition states that the particle concentration at ro is kept at 
zero, that is, the particle adheres. Integrating, subject to the conditions 
stated, the number of particles intercepted by the sphere is 


nt = UnronD (3) 


By way of simplification, we assume ro and D to be the same for all par- 
ticles. If each particle acts as a center at which coagulation takes place 
within a sphere of influence of radius ro, then 


(1/2) har jn@D = 2nrgn°D 


since two particles are involved to make a single particle. It follows from 
Eq. (3) that 


nto = — = = 2mronD (4) 


an equation which is integrable if ro and D can be shown independent of 
time. At this point, Smoluchowski proceeds along lines well reasoned in the 
kinetic theory of gases. Since diffusivity is an additive property, and since, 
because of size differences as coagulation proceeds, we have in the case of 
two particles, diffusivities Dj, and Dg, then we may write 


5 hs Dy + Lo 


Also, on the basis of kinetic theory of gases, the new sphere of influence of 
two coagulated particles will be 


Po = (1/2) (ro, + P69) 


in which roj and rog are the radii of influence of the two particles before 
coagulation. We note that if ro, and rog are equal to twice the radii of the 
particles, the particles touch each other sufficiently close to pull together and 
coagulate, If the radii of the sphere of influence are greater than twice the 
radii of the particles’ field of influence, then obviously the particles will al- 
ways coagulate. With these assumptions 


gn 
at 


n(D, 4 Dy ) ry * Fo jn (5) 
Einstein has shown that for particles in Brownian motion, the diffusion 
coefficient is 


D = 2/N . 


number and u the mobility defined by Stokes law of resistance 


where R is the gas constant, T the absolute temperature, N is Avogadro's 
u = 1/3mD, 
. 
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being the viscosity of the fluid. Thus, 


: D = (RD/X) (1/3mD,) 


ind from Eq. (5), if Dey and Deg are the actual effective diameters of the 
oagulating particles, we get 


= -1 pn cL 2 
-—= (RT/3LN) (D,3 + Doo ) ro4 + Too )n 


low again, as a further simplification, we may say that the radii of influence 
re as their respective diameters—that is, ro, = aDey and rog = aDe2 , 4 
ing a constant of proportionality. The last equation then becomes 


wane 
= x = (RTa /3uN ) (roy < Poo) n ry rO2 (6) 


t 


_ the sol or suspension is initially regarded as monodisperse and if as time 
roceeds the differences in sizes is not very great, then we may put roj = 
Io2 = etc., and we obtain an expression which is independent of particle size 


yay = {i /.) 


K 


(2a/3) (RT/uN) +t = Kt (7) 
(2/3) (RTa/uN) (8) 


ihere the limits of integration are obviously n=ng when t=0 and n=n 
then t=t. Equation (7) is identical to Eq. (2). However, it includes the 
ation of state and an unknown constant a which must be evaluated by ex- 
riment. It is indeed curious that Eq. (7) with all the simplifications and 

Sumptions made in deriving it should hold even for the early stages of 

agulation. But there is much evidence that it applies very well to many 
pagulating systems until gravity begins to take over. 

Muller (1926, 1928) has re-examined the general theory and has shown 
ow it can be applied to a system consisting of spherical particles of two dif- 
rent sizes. He has also indicated the modifications needed for dealing with 
pagulation of needle-like particles. Muller’s treatment of these special 
ies are quite involved. 

Fuchs (1934) appears to have been the last investigator to tackle the the- 
etical aspects of coagulation. However, he dealt with aerosols and by in- 
oducing concepts well-known in quantum theory advanced our knowledge of 
agulation a great deal. His results, however, while more precise, remain 
be tested. 

Equation (7) can also be applied to aerosols and this has been done with 
nsiderable thoroughness by Whytlaw-Gray (1932). For such applications, 
is necessary to introduce a correction in mobility 0 as defined by Stokes 
uation, which takes account of particle slippage between gaseous molecules. 
‘this case 


a= (1 + (24./0,)]/34D, 


5 ae 
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where A is the mean free path of the gas molecules, A is a parameter de- 
pending on the ratio A /De and other terms are as previously defined. The 
correction applied above is known as Cunningham’s correction and a discus- 
sion of its significance as well as values of A are contained in most texts 
on physical or colloid chemistry. 


Forces of Attraction 


Any discussion of the forces tending to hold particles together leads to the 
very nature of matter itself. The following discussion is intended to elucidate 
on mechanism. It is necessarily brief and intended only to give an idea of the 
magnitude of the binding forces. 

Let us assume that a particle in a coagulating system is momentarily fixec 
and is struck by another particle which has a kinetic energy of the order of 
10-8 ergs. This is a very small amount of energy easily acquired by a par- 
ticle 5 microns in size moving at a speed of about 10 cm per sec. But this 
amount of energy, small as it is, when given up in an inelastic collision with 
a cube of nickle 100 Angstroms on edge is sufficient to heat this cube to a 
temperature of about 100 deg C. (in air). Thus, if this energy were possessec 
by a particle of 100 Angstroms per edge and it struck a larger particle in- 
elastically, then it would appear that there is sufficient energy to cause the 
two particles to coalesce. Collisions in such instances “fuse” the particles 
together. 

But there are other forces of cohesion when particles come in contact that 
are referred to as van der Waal forces. These are forces exerted by unsatis 
fied atoms (as might occur at surface discontinuities of a freshly made crys- 
tal) and which on contact with similar atoms in another particle tend to bind 
each other together. These binding forces are of great intensity and are of 
the order of several tons per square inch. Thus, the probability of very smal 
particles joining together on contact is indeed very great. 


Particle Distribution Resulting from the Coagulation Process 


As coagulation proceeds, single particles will collide to form doublets. 
Single particles will collide with doublets to produce triplets, etc. Again 
doublets will be colliding with doublets to form quadruplets, or with triplets 
to form quintuplets, etc. Since Eq. (7) gives the law of diminution of particle 
numbers as coagulation proceeds, it is of interest to see whether it is possi- 
ble to determine the number of singletons, doublets, triplets, etc. present in 
a suspension at any time t after the process has begun. In this regard, we 
may again follow Smoluchowski. Let no as before be the number of particls 
initially present in a unit volume of the suspension, and n the number of par 
ticles present at time t and composed of 


Le 


nj particles which are singles 
n2Q particles which are doublets 
n3 particles which are triplets 


af, 


ae 
te. "ss | 


3 

= 

etc. 4 

nj particles which are i-lets 8 
a ar 3 


| 
ASCE DALLA VALLE 1052-7 


In deriving Eq. (7), we have assumed that the coagulation constant K applies 
to all the particles. If each collision of two unit particles decreases the num- 
ber of units by 2, while collisions with other sizes diminishes the number by 
one, then 


a4 BS x 2 
~dn,/at = 2Kn,“ + 2kn,n, + 2knyn, + oo. my 2knn 


| For the decrease of doublets formed by coagulation of 2 unit particles 


2 ake - 2 
edno/it = -Kn)° + 2Knon) + 2Knp + 2Knpnz +. - - = -2Kn,"+ 2Knon 


and for triplets 
ep Cone 
-dn3/dt = -2Knyny + 2Knzn, + 2Kn3no + 2kn,° + 2kn ny +. « 
= -2Kkn,n9 a ekn 3n 
Proceeding in this manner, we then obtain a series of expressions for i- 


particles. Since the rate of decrease of all the particles is controlled by 
| eq. (7), we may use this equation to eliminate time 


- dn/dt = Sic dn, /dt 


ta thus obtain differential equations with nj, ng,..., nj as variables. 
'These equations are easily integrated and we obtain the geometrical sequence 
| 


a né /n5 


No = n° (n, - n)/ng? 


ny = n° (ny = n)*/n3 


ee n°(n, - Ae ‘ fost (9) 


‘These equations, when added, form a geometric series with a common term- 
-term ratio of (no - n)/No. These are the most convenient form to work 
ith, although (following Smoluchowski) the following arrangement for the 
eneral term is also useful: 


my = n, Go) o/s @A jt * 


to = 3uN /uRTn ° 


which is of course, the reciprocal of K with a = 2, a generally accepted 
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value of this constant. Figure 1 shows how the ratio n/no changes with the 
ratio t/t with different values of nj as a parameter. We see, therefore, 
that if we have a means of determining the number concentration of particles 
at any time t and the initial concentration when t = 0, then the distribution 
of coagulated particles in terms of singletons, doublets, triplets, etc. is 
readily determined. 

We note in passing that the most frequent number of coagulated particles 
is given by the expression 


N = (2, - 2) (10) 
or, in terms of the coagulation constant K and the time t 
N = 2kKnt + 1 (11) 


Note also as shown in Figure 1 that the maximum of each grouping shifts to- 
ward the right as time increases. 


Sedimentation 


When the effect of gravity takes over, the coagulated particles will settle— 
provided there is no further accretions—in accordance with well defined laws. 
At the instant gravity may be said to take over, the density of the suspension 
at any point z measured vertically in the suspension and time t is con- 
trolled by Eq. (1).* 


Free Settling 


We shall now recognize that there are two kinds of settling: (1) that form 
in which the particles in the act of settling do not interfere with each others 
motion. This is usually the case when the volume concentration of the sus- 
pension is less than five per cent of the total volume. (2) Hindered settling 
which occurs with thick suspensions and is such that the particles settle as a 
mass with an interface between clear fluid and suspension. 

In free settling, the downward velocity of each particle is in accordance 
with Stokes law 


us = (P, - gD, ee = 2/t (12) 


where z is the distance moved by the particle without acceleration in the P 
time t and other terms have already been defined. If now in Eq. (1) z is ~ 
fixed, we obtain by differentiation of and application of Eq. (12) that | 


D 3 
oF) ee FD : 
dt 2t e) a3) Y 


in which @ is the ratio of concentration of particles at z at time t and a 


* Throughout this section time t=0 is taken at the moment coagulation 7 
ceases and gravity takes over. , 
. 
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oncentration when the particles were uniformly distributed in the suspen- 


ion (time, t = 0), that is, C(t)/Co=¢. If, on the other hand t is constant 
nd we wish to obtain the distribution at various levels in the suspension 


: ae (14) 
: 


din this case, Z refers to the ratio of concentrations at various z-levels 
the initial concentration, C(z)/Co. Equations (13) and (14) serve to give a 
omplete history of the free-settling phenomena. 


| 
indered Settling 


: 
_ In the case of hindered settling, the interface between the liquid and the 
uspension moves at a constant velocity ug. This problem may be handled 
1 much the same was as flow through a packed bed of particles with the par- 
icles moving instead of the liquid. In fact, the Kozeny-Fair-Hatch equation 
“the form proposed by L. Hatch (1943) may be used. This equation applied 
an expanded bed of fine particles with a porosity «. The form we use 

re applies to slow settling, but Hatch indicates well the general form of the 
quation for other types of flow. The equation is 


| 
| 


| (up?/ed,) [(1 -€)/E7] = 5.56 x 100, - ,) uD (45) 

| 

nere Dg is the specific surface diameter obtained by permeametric meth- 

ils (DallaValle, 1948). Other terms in the equation are defined elsewhere in 

lis paper. Note that as Eq. (15) is'written that the left hand portion is a _ 
dified Froude group and the right hand side a modified Reynolds group. 

is in general, we may expect the following equation to govern hindered 

Httling: 

| Fr, = @ (Re), (16) 


: m 


nis equation holds very well up to porosities of 0.85 and for particles rang- 
~ from 5 to more than 1000 microns in water, oil or glycol suspending 
edia. In fact, it is a convenient method for measuring the characteristic 
ameter Dg of particles. However, with coagulated or aggregated particles 
ich in settling carry along with them an adhered layer of fluid, density cor- 
ctions are necessary. In order to correct for this adhered liquid, it is 
cessary to follow along the lines indicated by Steinour (1944). 


SUMMARY 


By assuming initially a monodisperse system of particles or nuclei, the 
ince of collisions leading to 1, 2, 3, etc. aggregated particles is computed 
ording to the procedure used by Smoluchowski. The effect of gravitation 
ignored by assuming the process is very fast. A discussion is included of 
pb general law of settling of dilute suspensions as well as highly concen- 
ted suspensions. 
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HYDRAULICS DIVISION 


___ Proceedings of the American Society of Civil Engineers 
| SSS 


NEWS BULLETIN-August 1956 


) The April issue listed the committees for the Society business year 1955- 
1956 and referred to the Official Register for the names of members serving 
on Task Forces. The following additions and corrections should be made in 
your Register to make the record of those at work in the Division complete: 


Task Force On Energy Dissipators for Spillways and Outlet Works: 
Add Rex Elder 


Task Force on Hydraulic Aspects of the Design and Operation of Spillways 
| and Crest Gates: 

| Alvin J. Peterka is now chairman 

| Marvin J. Webster is now a member 

Messrs. Nalder and Brown resigned. 


Task Force on Nomenclature for Hydraulics: (Note the change in the 
name of this task force) 
C. A. Betts is also on this task force 


: Task Force on Water Measurements: 
| Add A. L. Jorissen and C. A. Wright 


Task Force on Sediment Distribution in Reservoirs: 
Samuel Shulits, chairman 
Paul C. Benedict 
Lloyd C. Fowler 
Carl R. Miller 
John B. Stall 


The Hydrology Committee has established a Task Force on Hydrologic 
Data. The Task Force is staffed by the following members: 


D. W. Van Tuyl, chairman 


C. O. Clark W. W. Wheeler Victor A. Koelzer 
L. L. Harold Robert R. Balmer Edward J. Schaefer 
H. E. Hudson Verne Alexander Clayton H. Hardison 


Walter G. Schulz 
C. C. McDonald, Coordinator 


A Task Force on Spillway Design Floods has been authorized. 


Harvey O. Banks is chairman 


Note: No. 1956-18 is part of the copyrighted Journal of the Hydraulics Division of the 
American Society of Civil Engineers, Vol. 82, HY 4, August, 1956. 
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The Committee on Research has two task forces at work. The Task Force 
on Aerated Flow in Open Channels consists of: 
Warren W. De Lapp, chairman 
William J. Bauer 
Donald Colgate 
Joseph P, Lawrence 
(Frank B. Campbell is Contact Member) 


The Task Force on Cavitation in Hydraulic Structures consists of: 
F. R. Brown, chairman 
R. A. Elder 
V. E. Johnson, Jr. 
M. B. McPherson 
A, J. Peterka 
(J. M. Robertson is Contact Member) 


Mr. Thornton J. Corwin, Jr., past chairman of the Executive Committee, 
recently resigned from the committee for reasons of health, It is regretted 
exceedingly that illness has made it necessary for him to give up his work in 
the Division. A recent note from Thor indicates that he is making a satisfac- 
tory recovery and is now able to be at his office in the Pacific Gas and Elec- 
tric headquarters. 

Dr. Arthur T. Ippen has been appointed by the ASCE Board of Direction as 
a member of the Executive Committee of the Hydraulics Division. This ap- 
pointment is effective immediately for a term ending in October of 1960, Art 
brings with him a wealth of experience in hydraulics and committee work 
from which the Division will undoubtedly profit. 


Committee and Task Force Meetings 


The following meetings have been held so far this year since October 1955 


Executive Committee Denver November 7-8, 1955 
Committee on Tidal Hydraulics Washington, January 17, 1956 

D. 
Committee on Hydrology Dallas, Tex. February 13-14, 1956 
Committee on Design Knoxville, May 3-4, 1956 

Tenn. 
Committee on Sedimentation ub June 4-5, 1956 


Task Force on Sediment Distri- Urbana, Il, May 15-17, 1956 
bution in Reservoirs 


The Committee on Research will hold its next meeting at Chicago, 
August 25, 1956 


The Executive Committee will hold its next meeting at Madison, 
Wisconsin, August 21, 1956 


HYDRAULICS DIVISION TO HAVE BIG PROGRAM 
AT PITTSBURGH CONVENTION 


asc Hydraulics Division 1956-18--3 
| October 15-19, 1956 


Place: William Penn Hotel 


Technical Sessions have been formulated by the following committees: 


Design 1 session 
| Floods 1 session on flood insurance 
Research 1 session on turbulence in fluid flow 
1 session on cavitation in hydraulic 
structures 
1 general session 
Tidal Hydraulics 2 sessions sponsored jointly by the 


Committee on Tidal Hydraulics and 
the Committee on Coastal Engineer- 
ing of the Waterways Division 


Committeemen and Task Force chairmen are invited to use the Hydrau- 
ics Division News Bulletin to: 


a) Insert outline of proposed work and invite suggestions thereon 


| b) Insert resume of work or report and invite suggestions thereon 


| ¢) Insert questions on specific problems on which information is required 


| 
) Committee on Publications 
| 


Chairman H. G. Dewey reports that since the publication of JOURNAL No. 

2, the Division has passed on 14 papers to the Manager of Technical Pub- 

ications for publication in future issues of the JOURNAL; eight (8) papers 

ire being revised by the authors; and three (3) papers are under active re- 

tp . This indicates another very active year in publication of Hydraulic 
ers. 


| 
) HYDRAULICS DIVISION MEMBERSHIP 

Questionnaire 

Repo 

| It has been several years since the membership of the Hydraulics Division 
been canvassed to determine the interest and availability of the individual 
embers for service in Division activities, All who are interested in par- 
cipation in Division activities are urged to fill out the following question- 
ire and return to the Secretary of the Division. 
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QUESTIONNAIRE OF PARTICIPATION OF HYDRAULICS DIVISION 


Name 
Address 
Position Division Committee Service 
Consultant a Design —s 
Manufact’r pens Hydrology _—_— 
Government 2 ae Research Ss 
School Barat Publications — 
Private Soe s Floods ——— 
Other abate toed Tidal Hydraulics — 
Retired SS Sedimentation —s 
Stevens or 
Field of Interest Hilgard Prize _ 
Other Member —= 
Design et 
Hydrology aes Review Papers 
Research =a se ae 
Ind. Hydr. a Yes —- 
Other ees No —- 
Floods ss 
Preparation of Technical Paper Sedimentation = 
Gr’d Water —- 
By self San Snow, Ice and 
With other eats Permafrost 3 
Design a 
Edit Papers Hydrology —s 
Yes dene Research General = 
No Toes Research Special- 


ized (Indicate 
specialty) 


| 


Comments and suggestions 


Subjects recommended for discussion at meetings 


—_— 
7 


——_— s 
Suggested Division meeting place a 


Please return to Harold M. Martin, Secretary, Hydraulics Division, ASCE, i 
Bureau of Reclamation, Denver Federal Center, Building 56, Denver 2; 1 
Colorado. 4 
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k | PROCEEDINGS PAPERS 


© *@chnical papers published in the past year are identified by number below. Technical- 
Sponsorship is indicated by an abbreviation at the end of each Paper Number, the 
referring to: Air Transport (AT), City Planning (CP), Construction (CO), Engineering 
(EM), Highway (HW), Hydraulics (HY), Irrigation and Drainage (IR), Power (PO), 
by Engineering (SA), Soil Mechanics and Foundations (SM), Structural (ST), Surveying and 
(SU), and Waterways and Harbors (WW) divisions. Papers sponsored by the Board of 
i are identified by the symbols (BD). For titles and order coupons, refer to the appro- 
“assue of “Civil Engineering.” Beginning with Volume 82 (January. 1956) papers were 
d in Journals of the various Technical Divisions. To locate papers in the Journals, the 
etter the paper numbers are followed bya numeral designating the issue of a particular 
z 2 |in which the paper appeared. For example, Paper 861 is identified as 861 (SM1) which 
es that the paper is contained in issue 1 of the Journal ofthe Soil Mechanics and Founda- 
Wision. 
‘ 2 VOLUME 81 (1955) 


; 
i 


be 


usr: 761(BD), 762(ST), 763(ST), 764(ST), 765(ST)°, 766(CP), 767(CP), 768(CP), 769(CP), 
"O(CP), 771(EM), 772(EM), 773(SA), 774(EM), 775(EM), 776(EM)°, 777(AT), 778(AT), 
f (SA), 780(SA), 781(SA), 782(SA)°, 783(HW), 784(HW), 785(CP), 786(ST). 


i: [BER: 787(PO), 788(IR), 789(HY), 790(HY), 791(HY), 792(HY), 793(HY), 794(HY)°, 
(EM), 796(EM), 797(EM), 798(EM), 799(EM)©, 800(WW), 801(WW), 802(WW), 803(WwW), 
J4(WW), 805(WW), 806(HY), 807(PO), 808(IR)°. 


: 809(ST), 810(HW)©, 811(ST), 812 (ST)°, 813 (ST)°, 814(EM), 815(EM), 816(EM), 
IM), 818(EM), 819(EM)°, 820(SA), 821(SA), 822(SA)°, 823(HW), 824(HW). 


BMBER: 825(ST), 826(HY), 827(ST), 828(ST), 829(ST), 830(ST), 831(ST)°, 832(CP), 
33 CP), 834(CP), 835(CP)S, ’836(HY), 837(HY), 838(HY), 839(HY), 840(HY), 841(HY)°. 


IER: 842(SM), 843(SM)°, 844(SU), 845(SU)°, 846(SA), 847(SA), 848(SA)°, 849(ST)°, 
sr), 851(ST), 852(ST), 853(ST), 854(CO), 855(CO), 856(CO)°, 857(SU), 858(BD), 859(BD), 
D). 
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Y: 861(SM1), 862(SM1), 863(EM1), 864(SM1), 865(SM1), 866(SM1), 867(SM1), 868(HW1), 
1), 870(EM1), 871(HW1), 872(HW1), 873(HW1), 874(HW1), 875(HW1), 876(EM1)°, 877 
, 878(ST1)°. 


¥: 879(CP1), 880(HY1), 881(HY1)°, 882(HY1), 883(HY1), 884(IR1), 885(SA1), 886(CP1), 
Ai), 888(SA1), 889(SA1). 890(SA1), 891(SA1), 892(SA1), 893(CP1), 894(CP1), 895(PO1) 
P01), 897(PO1), 898(PO1), 899(PO1), 900(PO1), 902(PO1), 902(AT1)°, 903(DR1)°, 904 
a 905(SA1 ne 


¥ 906(WW1), 907(WW1), 908(WW1), 909(WW1), 910(WW1), 911(wW1), 912(WW1), 913 
, 914(ST2), 915(ST2), 916(ST2), 917(ST2), ie eed eee 920(ST2), 921(SU1), 
1), 923(SU1), 924(ST2)°. 


926(WW2), 927(WW2), 928(SA2), 929(SA2), 930(SA2), 931(SA2), 932(SA2)°, 
sr ouscwi2), 936(WW2), 937(WW2), 938(Ww2), 939(WW2), 940(SM2), 941 
mene)” 943(EM2), 944(EM2), 945(EM2), 946(EM2)°, 947(PO2), 948(PO2), 949(PO2), 
» 951(PO2), 952(PO2)¢, 953(H¥2), 954(HY2), 955(HY2)°, 956(HY2), 957(HY¥2), 958 
9(PO2), 960(PO2). 


5 , 962(IR2), 963(CP2), 964(C'P2), 965(WW3), 966(WW3), 967(WW3), 968(WW3), 969 
} ers, 971(ST3), ome pues ona(ars, 975(WWw3), 976(WW3), 977(IR2), 
°2) aco. en ea aie ca 983(HW2), 984(HW2), 985(HW2)°, 986(ST3), 


3), a , 993(PO3), 994(PO3), 995(PO3), 996(PO3), 997(PO3), 998 

9) Be ‘ na tere eran, 1003(SA3)°, 1004(H¥3), 1005(HY3), 1006 

(HY3), 1008 (HY3), 1009 (HY3), 1010 (HY3)¢, 1011 (B03)°, 1012 (SA3), 1013 (SA3), 

15(HY3), 1 1016(SA3), 1017(PO3), 1018(PO3). _ 

ie: 1021(ST4), 1022(ST4), 10231874), 1024(ST4)°, 1025(SM3), 1026 
oe “1030(E3), 1031 (8813), 1032(EM3), 1033(EM3)°. 


F 1039 1040(HY4) 
1037(HY4), ios’, ast utisas, oan 
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